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"Über das Restglied der Mittelwertformeln für angenährte 
ro Ra  Quadratur | | ER 
ER r Sr None Tollmien in Göttingen 


Für die Quadraturformeln nach dem Mittelwertwerfahren werden die exakten, Restglieder in Integral- 

form auf eine möglichst einfache Weise hergeleitet. Diese Restglieder enthalten von der nüherungsweise zu 

RI NEN Funktion Ableitungen beliebig wählbarer Ordnung sowie die periodischen Bernoullischen 
ı Funktionen. ee IE 


Ri The remainder terms of quadrature formulas ewpressed in. the exact form of integrals are ded uerd ina 
most simple manner. Derivatives of Ihe approwimately integrated function of an arbilrary order and the 
periodic Bernoulli’s functions are contained in Ihese vemainder terms. 2 


B 5Toi pa6ore BEIBONAITCA HAMÖOoNTee IPOCTEIM O6PA30M TOYHBIE BEIPAKEHHSA, HMEIOIIME 
POPMy HHTETpana, ANA OCTATOYHBX WIeHOB B POPMYıax AI NPUÖNNKeHNOrO BbLyncneHnA 
. HHTeTPaJI0B. OTU OCTATOYHBIE YAIEHBI CONEP:KAT IPOH3BONHEIE IMWOOTO HOPAAKA OT NOAHHTC- 
TpansHoN# AyKılmu N KpoMe Toro nepmonuueckue hynkumu Bepnyann. 


Die Mittelwertverfahren zur angenäherten Quadratur approximieren das Integral durch 

ein lineares Aggregat von Werten des Integranden für einige ausgezeichnele Argumente, so daß 

_ eine solche Formel bei einem über das Intervall 1 erstrecklen Integral in der folgenden Weise 
geschrieben werden kann: Be 


1 N f 4 
[ Flu)du = 23.4, F(u,) , DEE U EN EN TRETEN NER 
0) \ v=1 


Die Gewichte A, werden nun so bestimmt, daß die Quadraturformel Polynome bis zum Grade 
m—-1 exakt integriert, woraus m Beziehungen zwischen den Gewichten A, und den ausge- 
zeichneten Argumenten u, folgen: 7 


Am = 7 DEINER EN Be 


Bei den meisten Quadraturformeln, etwa denen vom Cotesschen Typ, ist m=n oder n-+1, bei 

den Gaußschen Formeln ist m =2n. 

Es ist sehr bekannt, daß man das Restglied der Euler-Maclaurinschen Quadraturformel 
exakt mit Hilfe der Bernoullischen Polynome!) ausdrücken kann. W. Wirtinger?) hat seine 
Herleitung dieses Zusammenhanges sinngemäß auf die Feststellung eines exakten Restgliedes 
für beliebige Quadraturformeln, die nach dem Mittelwertverfahren arbeiten, übertragen. Das 
2 . Wirtingersche Restglied enthält die m-te Ableitung des Integranden F und die periodischen 
“ Bernoullischen Funktionen?), diesich aus den Bernoullischen Polynomen im Grundintervall O1 

durch periodische Fortsetzung ergeben. 

SE Die nachstehenden Bemerkungen, die ich im Sommersemester 1948 in Göttingen vortrug, 
verfolgten das einfache Ziel, das genannte Restglied möglichst direkt herzuleiten. Insbesondere 
wird der üblicherweise eingeschlagene Weg über die. Fourierenlwicklungen der periodischen 
Bernoullischen Funktionen vermieden. Weiter wird das Restglied mittels einer Ableitung von # 
von niedrigerer als der m-ten Ordnung ausgedrückt. Da die höheren Ableitungen der Integranden F 
oft schwer zugänglich sind, besteht eine dringende Notwendigkeit, derartige Restglieder mit 
beliebig niedrigen Ableitungen von F aufzustellen, worauf R. v. Mises*) besonders hingewiesen 
hat. Die Lösung dieser Aufgabe wurde wesentlich durch das anregende Buch von G. Kowa- 


} 
. 


2) Vgl.etwaK.Knopp, Theorie und Anwendung der unendlichen Reihen. 3. Aufl. Berlin 1931, 8. 541. 

2) W. Wirtinger, Über das Fehlerglied bei numerischer Integration. Z. angew. Math. Mech. 13 (1933), 
S. 166—168. 

3) In der Bezeichnungsweise schließe ich mich dem Werk von J. F. Steffensen, Interpolation, 
Kopenhagen 1925, Baltimore 1927, $ 14, an. 

4) R. v. Mises. Zur mechanischen Quadratur. Z. angew. Math. Mech. 13 (1933), 8. 53—56. 
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lewski®) gefördert. In den folgenden Ausführungen ist wenig Gewicht auf Neuartigkeit 
sondern a Hauptanliegen ist die Klärung des Zusammenhanges verschiedener estgliedfo 
und die möglichst naturgemäße Einführung der periodischen Bernoullischen Fur 1 


% 
= Ye 
>) 
Bi 


1. Aufstellung eines ersten Restgliedes 


Wir gehen bei unseren Betrachtungen von einem Restglied aus, das die erste Ableitung 
von F enthält. Da u 


ie 1 1 n I; 

| F(u,) = F(0) + ALHOLEE / F(u)du= F(0) + fr (1—2)dz, . Z4u=1 
ist, so wird dies Restglied 
4 N U, s 
R= / Fe) (1-2) — & 4, / F'(e) dz 


nn 


n . u, 
= 34, | [Ft — 2) de— f F'‘(«) ae) 
v=1 d 
Um ein einziges Integral zu erhalten, wird die Funktion o eingeführt: 
5 o(Ww—2)=1 für 0 <sz Sf (4) 
et ee Pe ‘ 
Damit wird 
LE R 
= [F@24 (1-2) w—Dldz ... 2.2 22.. (5). 


Dies Restglied soll nun mittels der periodischen Bernoullischen Funktionen ausgedrückt und 
durch partielle Integration umgeformt werden. . 


2. Einführung und Eigenschaften der periodischen Bernoullischen Funktionen | 


Die Bernoullischen Polynome B,(v) werden durch drei Festsetzungen definiert: | 
a Be, ER (6) | 

(240) dv=(, BET Be (7) 

Belt. ee (8) 


Die nächsten vier Bernoullischen Polynome lauten danach: 


Bor B,b) = W—v BT 


>) ’ 
De B RR RL Be el.) 
2 2.9 a 30 ' 
Die Werte der Bernoullischen Polynome bei »=0 sind die Bernoullischen Zahlen Bz: 
i 1 1 1 
il, B=—7>: B=7: B=%V, B=—-5 B,=0*.-- 


Einige sich fast unmittelbar aus den Definitionsgleichungen ergebende Eigenschaften der 
DL N Polynome im Grundintervall 0 <v <1 seien hier zusammengestellt. Da für 
> 
1 


| B-1()do=0=-[B,1)— 3,0], 


so ist 
3A )=B,0)=B für 523 


5) G.Kowalewski, Interpolation und genäherte Quadratur. Leipzig und Berlin .1932. 


Me ge Ye: R r n 
k ionen von ur Er sind, die Bernollischen Polynome ungerader Ordnung dagegen ungerade, 


Es 
x Funktionen vo von %—z: mit der Beziehung. [C)% ergibt die letztere Eigenschaft, ARTEN 
Br die. Bernoullischen Br 'olynome ungerader Ordnung außer bei » = 2 auch bei vo 0 und o—1 | 
_ verschwinden, wobei B,(v) auszunehmen ist. Daß die Bernoullischen Polynome ungerader Ar 


Ordnung keine weiteren Nullstellen außer denen bei v»—0, z und 1 haben, läßt sich leicht aus. 
E dem Widerspruch zeigen, den die Annahme einer weiteren Nullstelle nach dem Rclleschen Satz 
_ ergeben würde, wenn man gemäß (6) von Bernoullischen Polynomen höherer Ordnung zu solchen 


‚niedrigerer Ordnung herabsteigt. Daraus folgt, daß die Bernoullischen Polynome gerader Ordnung 
im Inneren des Grundintervalls ein einziges Extremum und zwar bei v =7 haben. Hieraus“ | 
kann man den für Abschätzungen wichtigen Schluß ziehen, daß | 


Ba (0) — Ba; 


im Grundintervall nekche Vorzeichen hat. Beiläufig sei Be, daß die Bernoullischen 
Polynome gerader Ordnung zwei Nullstellen im Grundintervall haben. Das Aussehen der ersten eA 
Bernoullischen Polynome kann für das Intervall O<vsi1 der Abb. 1 entnommen werden. Rn 


“Die periodischen Bernoullischen Funktionen Br (v) stimmen im Grundintervall mit den = 
Bernoullischen Polynomen überein und ‚haben die Periode 1: ‚2 A 


B4(v) = B.(v) nur < 1 
Bw +1) = By(v) für alle v. 


B;(v) hat Sietige Differentialquotienten bis zur 
Ordnung k—2, Unstetigkeiten im k—1-ten Differen- 
tialquotienten bei v=0, +1, +2, +3,... Die perio- 
dischen Bernoullischen Funktionen gerader Ordnung 
sind gerade Funktionen von v, die ungerader Ordnung 
ungerade Funktionen. Bild 1 zeigt die ersten vier Cr 
periodischen Bernoullischen Funktionen. 


z zZ 


1 
7-2-O(u,-Z) a; ES u} 
Brti uy-2) 


Bild ı and - 


3. Aufstellung der allgemeinen Restgliederung in Integralform 


Nach diesen Vorbereitungen, welche die periodischen Bernoullischen Funktionen der An- 
schauung möglichst nahebringen sollten, ist es leicht möglich, die Restglieder der Quadratur- 
formeln mittels periodischer Bernoullischer Funktionen auszudrücken. Für die in dem Kern 
des Restgliedintegrals R, nach (5) auftretende Funktion gilt: 


120m) + Bm) Re See (10), 


2 ZB he en a N Aa un aa za DD ma a sin ne au 


' Nun ist wegen der Periodizität 


und wegen der exakten Gültigkeit der 'Quadraturformel für Polynome 2. Grades: 


ist, so erhält man X et En RE 
| =/ro Een A: 


Durch partielle Integration erhält man die weiteren nern, welche ar 
leitungen von F enthalten. Zunächst bekommt man % 


R=—zF) 2 >>} A,B, (u,—1) +50 2 2 ArBata) u Fr Be 


ee @) EuBw— 2 de. 


B,(w—1)= B, (w,), 


1 


24 B,(u,) = | B,(u)du=0, 


wenn m—1 <2 ist. Danach wird. 
2 
Klar NE RE DE 
2-4 [P% 2 AyB, (up — 2) de DW ET Tr Par Sa a, © 3 rn" (12). 
ö x 


Durch Fortsetzung dieser Schlußweise erhält man allgemein 
1 f a - 
Rs Lu pi I» (2) 2 4, B, (U, — 2) dz TE a ar BAT En (13), 
y- 


falls p<&m— 1. Damit ist das Restglied in Integralform mittels einer beliebig niedrigen 
Ableitung von F ausgedrückt, was erstrebt war. 


Will man das Restglied mittels der m-ten. Ableitung von F ausdrücken, so hat man 
das letzte Restglied der Form (13), nämlich 


1 
1 we 
Bu PO E Ann (u, — 2) dz Car Pa Te. .,; (14) 


(m 


durch partielle Integration umzuformen. Man erhält: 
1 
Auen 1 m— m— S 1 5 ir: 
ee a We Bm (4) + Si Jr ZrBn (w—2)dz (15). 


‘Das Polynom B„(w), das mit w” als höchster Potenz beginnt, wird durch die Quadraturformel 
nur mehr angenähert integriert, und zwar ist der Fehler: 


1 \ 
[rau Au = 
11) vl 7 


N. 
FRE ya insb ©. m 
nr Pa ? 


Also ist 


1 
n 1 n 
2, Ar Bmw) anaiTz Ara | Batadu=0 A EN naar (16). 


ER TTE 


ae u 


Folglich ist 


Rn = ee] ih 


Be vi 2 
Hamm (2) P> A, Bu (u, — 2) dz 
wofür man auch schreiben a 


1 
An [Fo e,, en wa Se. A,ur +34 Bu —a)} 2: ES 


4. Einige Anwendungen 
Einige Anwendungen der hergeleiteten Restgliedformeln auf die Simpsonsche Formel und 
die Trapezregel mögen zur Erläuterung zusammengestellt werden. 
Für die Simpsonsche Regel ist 


u,=0, w=4, Tr, 
1 2 1 
4A,= % = en A, Ya 


und m=4. Nach u wird das Restglied mit der zweiten Ableitung von F: 


R,= - [md 148 B,(— 2) +58 B B,(2 2) nor I nualar 


kn 
ee 2 —- {1 
== [r (z) Sn B,(z) +7 B, (4—:)}ae 

Bezeichnen wir den Kern dieses Restgliedintegrals R 

v(2) = 2 B,(z) en PD, Me )\ NR PER AN CH (20), 
so wird 
(dw ee 
Ya(2z) = 32273 ür 2<Sz 

a RL N (21). 


1 2 
)=ze@—1) («5 für 225 


Auf anderem Wege war diese spezielle Formel von G. Kowalewski l.c. S. 136 hergeleitet. 
ys(z) ist in Bild 3 dargestellt. 


Will man das Restglied der Simpsonschen Formel mittels der vierten Ableitung von F 
ausdrücken, so erhält man nach (17): 


1 
1 — {1 
= 2 vn RW vH 2B we 2) d END 
Das Integral kann man durch partielle Integration noch weiter umformen, wonach man schließlich 
unter Beachtung der Nullstellen von B,(z) auf die Restgliedform: 


1 
Pl 
ER R0 vo Bu —2B,(5 —2)}dz. ...(@83) 


er 


“Formel noch Polynome 5. Grades exakt integriert werden, da 


i (& - 2 N nt ” ’ e ö 
« £ } w N a dur i % fer u Pe 
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kommt. Der erste Teil ist das Chevillietsche Restglied, mit dessen 


£ 


nz 
x 


r 2 "ii, 

ee R k 
a “ 

vr de er 


| [240 _13, (U-)}ae=o0 


Bei der Trapezregel ist = 0,» =1, 4, = 3’ A= DE m=2. Nach; dnETIANNE 
gliedformel (15) erhält man | 


wird. 


1 . 1 
| Fi) du OFF _ Br + / F'(e) Bz()dz . . . (M). 


Indem man das letzte Integral durch partielle Integration entwickelt und das Verschwinden 
der Bernoullischen Zahlen ungerader Ordnung beachtet, erhält man die Euler-Maclaurinsche 
Quadraturformel: 


Fee | 
[rw a0 _ Brno 
{0} r 
un var Dr: 3 22-1) (1) — FA—N) (0 En 
— Zr (1) — F"”(0)] an?“ (1) — F@-2) (0)] (25) 
1 


1 
Sa ee Baizı (2) dz 


Für das exakt in Integralform angegebene Restglied kann man leicht eine Abschätzung geben. 
Es wird durch partielle Integration: 
1 . 
1 1 | 
a Fertl) = _ ——___ Bars [FRA KT) — FRIHDIO 
am? (2) Barzı ) de = — ggg; Bass [F@3H00N) 0) 


e 
| 
1 


ı 1 
1 1 
+33 er Bai+2(2) de= EEE aaa [Berrel2) —Barreldz. 


Aue 


Da-nach einem der aufgeführten Sätze über Bernoullischen Polynome gerader Ordnung 


Bar+2(2) — Barr2 


einheitliches Vorzeichen hat, so wird nach dem Mittelwertsatz der Integralrechnung das Restglied 
gleich 


1 
1 B 
PETE | [Barrat) — Bar lde= — a7 291 Fett)... (26), 


wo0 s’sıt ist. 


5. Schlußbemerkung 


Es sei bemerkt, daß U. Wegner‘) auf anderem Wege auch das Restglied der Lagrangeschen. 
Interpolationsformel mit Hilfe der periodischen Bernoullischen Funktionen ausgedrückt hat, 
woraus natürlich die Restglieder der Quadraturformeln ebenfalls gewonnen werden können. 
Nachdrücklich sei auch auf den Vortrag von E. Weinel?) auf der Göttinger GAMM-Tagung 1948 
über das Restglied der Lagrangeschen Interpolationsformel und der Formeln zur numerischen 
Differentiation und Integration hingewiesen, wobei aber in den Integraldarstellungen der Rest- 
glieder die periodischen Bernoullischen Funktionen nicht zur Anwendung kamen. 


& ) “2 Wegner, Zur Theorie der Interpolation. Dtsch. Mathematik, Jahrg. 1 (1936), S. 103—107 und 
3. 352— 858. 

?) E.Weinel, Über das Restglied der Lagrangeschen Interpolationsformel und der Formeln zur nume- 
rischen Differentiation und Integration. Vortrag auf der Göttinger GAMM-Tagung, September 1948. Z. angew. 
Math. Mech. 29 (1949), S. 32/33. 


Eingegangen: 25. 10. 1948, 
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Differenzenverfahren zur numerischen Integration von 
gewöhnlichen Differentialgleichungen n-ter Ordnung 

| | Von L. Collatz in Hannover s; % | 

Es werden zur numerischen Haan von Anfangswertaufgaben bei gewöhnlichen Differentialglei- 


Ordnung als Spezialfälle enthalten. Unter den verschiedenen Verfahren zeichnet sich auch bei Differential. 
gleichungen höherer Ordnung das Verfahren der zentralen Differenzen durch einfache Koeffizienten. und. 
günstigere Konvergenzbedingungen aus. i 


For the numerical compulation of initial value problems of ordinary differential equalions of Ihe n-Ih 


order, two methods of extrapolation and two others of interpolation are staled, which conlain most of the 
difference methods hitherto known for differential equations of the first and the second order as special cases. 
4A ' the several methods, the central difference method is distinguished by simpler coefficients and more 
favourable criteria of convergence also in the case of differential equalions of @ higher order. 


' 


na vucnekHoro pemeHuna IIPOCTEx ndepenmmassesrx ypaBHennf n-HOTO NOPAAKa C 
HAYAJIbHBEIMH TPAHHYHBIMH YCHOBHAMH ANAMTCH ABA DKCTPONONANHOHHHX MH ABA HHTEP- 
HONANHOHHBIX METONA, KOTOPHEe Conep:KaT B cebe B BHNe WACTHHX CIIyyaeB IIOYTH Bce 
cymecTByiolnue 10 CHX HOP METOANkI PeImeHusd B KOHEYHBIX PASHOCTAX AuhepeHNHanbHEIX 
ypaBHeHnf mMepBoro H BTOPOrO TopAnka. M B cıyyae andepeHnHuanbHkIX ypaBHeHunt Öorree 
BBICOKOTO HOPANKA METOAN HeHTPANbHEIX KOHCYHBIX PA3HOCTEH OTIHYAETCAOT APYTUX METONOB 
'Haubonpımeä IPOCTOTOH KoshunneHroB u Öoree YNOÖHBIMAH IIPUSHAKAMH CXONHMOCTH. 


1. Einführung | 
Für die numerische Integration gewöhnlicher Differentialgleichungen sind die Verfahren, 


. welche Differenzenschemata benutzen, in steigendem Maße verwendet worden, und zwar haupt- 


sächlich die Extra- und Interpolationsverfahren von Adams, Störmer, Nyström u.a. für 
Differentialgleichungen erster und zweiter Ordnung. Im folgenden sollen verschiedene Arten 
der Durchführung für Differentialgleichungen n-ter Ordnung 


DIÄTEN YES Ye DD (1) 
besprochen werden: 
Es liegt an sich nahe, eine Differentialgleichung n-ter Ordnung auf ein System von n Glei- 
chungen erster Ordnung zurückzuführen: 


Y=Yo» %=Y %=lo -:» I-ı= Ya Yn=fl® Yo Yo Ya --» Yo)» 


und dieses mit den für Differentialgleichungen erster Ordnung bekannten Verfahren zu behandeln. 
Jedoch wird dadurch nicht nur die Rechenarbeit größer (z. B. hat man dann n Differenzschemata), 
sondern es leidet dabei auch die Genauigkeit; beim Differenzenverfahren sinkt die Genauigkeit 
infolge des n-maligen Einsetzens der Funktionsverläufe durch Polynome um ganze h-Potenzen 
der Taylor-Entwicklung!). So bietet die Behandlung einer Differentialgleichung »-ter Ordnung 
grundsätzliche Vorteile gegenüber der eines Systems von Differentialgleichungen erster Ordnung, 
undes ist zuempfehlen, ein System von Differentialgleichungen niedriger Ordnung zur numerischen 


"Behandlung in ein System von weniger Differentialgleichungen höherer Ordnung zu verwandeln, 


sofern diese Umwandlung analytisch bequem durchführbar ist, und dabei die Funktionsausdrücke 
in den Differentialgleichungen nicht etwa zu kompliziert werden. 
Zur Festlegung einer bestimmten Lösung der Differentialgleichung (1) seien für =, 
Anfangswerte vorgegeben: 
vayed. Yan. Frey d.: ine.) 
Die Funktion f der Differentialgleichung (1) genüge in einem konvexen Bereich ® des x, %, Yy, 
..., Y®—D-Raumes, in dem die Lösungsfunktion und alle Näherungen enthalten sind, einer 
Lipschitzbedingung 
s * %* * 
a en ER TEE Eos) EEE ER (8) 


n—1 * 
< z K, |y» — yo 
v= 
für alle Punktpaare 9, Y, ..., y*=U und y, y, EN ya—d aus ®. Für praktische Zwecke 


2) L.Collatz— R. Zurmühl, Zur Genauigkeit verschiedener Integrationsverfahren bei gewöhn- 
lichen Differentialgleichungen. Ing.-Arch. 13 (1942), S. 34-36. D 


chungen n-ter Ordnung zwei Extrapolationsverfahren und zwei Interpolationgverfahren angegeben, welche i 
. die meisten bisher bekannten Arten der Differenzenverfahren für Differentialgleichungen erster und zweiter 
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Amy: f u 0 a" uaN se 
wird man oft setzen: | ee: ER yo 
Be En 


Überhaupt verlaufen die Überlegungen bei Differentialgleichungen höherer Ordnung ähnlich = | 


denen bei Gleichungen erster und zweiter Ordnung?), so daß wir uns hier kürzer fassen können. 

Wie bei.den Gleichungen erster und zweiter Ordnung besteht die Rechnung aus zwei ganz 
getrennten Schritten: 1. Man muß sich zunächst mehrere Näherungswerte für y und seine Ab- 
leitungen und damit für die Funktion f, die in der Differentialgleichung (1) auftritt, verschaffen 


(Berechnung des „‚Anfangsstückes‘‘). 2. Bei der ‚fortlaufenden Rechnung‘ werden die Werte 


von y und seinen Ableitungen an der Stelle z,+ ermittell, wenn die Werte bei =, &-ı1, 
&—2, ... bekannt sind. Es werden mehrere Rechnungsarten angegeben 


2. Hilfstormeln aus der Differenzenrechnung « : 
Um später die Gedankengänge nicht unterbrechen zu müssen, stellen wir hier einige nach- 


"her gebrauchte Hilfsformeln zusammen. Es seien fz=(#;) die Funklionswerte einer Funktion f(z) 


an äquidistanten Stellen z, = a, + kh (für k= 0, 1,2,..., N, manchmal auch nicht ganzzahlig) 
mit h als Maschenweite. Mit den üblichen Differenzenoperatoren 


ee (5) 
vh=hh-ı 


hat man die bekannten Interpolationspolynome (mit der Abkürzung u = =, ): Das ‚‚New- 


tonsche Interpolationspolynom mit absteigenden Differenzen“: 


i 1 1 be 
tert: mini 


mit dem Restglied?®) 
Ba 1 (x) =f(&) Fr N,(8) FR 


PPfo - - (6) 


u(u +1)---(u+p) 
(p+1)! 

(mit & als geeigneter Zwischenstelle in einem Intervall der x-Achse, welches die Punkte z,, 2_,, 

-., £_pund zenthält), das ,‚Newtonsche Interpolationspolynom mit aufsteigenden Differenzen“: 


PHlfo+DE 2.2.2... 


u (u — 


R 1 1 
af + HA + Ap+e.- +) et )erfene. . @) 


und das Stirlingsche Interpolationspolynom, welches wir nur für gerade p aufschreiben: 


UNE RE rh+rfhi, w uw—1)pehh +r°f, 
Sl = UT Tata ee 
+ ————prpn 


init dem Bestlglied 


FERN a 

R,+1(2) =f (2) — St,(a) | Frlaryerng ti). 
Br iH 

Aus dem Stirlingschen Interpolation»polynom gewinnt man die symmetrische Darstellung®), bei 


ai It. Y. M in 3 Zur numerischen Integration von Differentialgleichungen. Z. angew. Math. Mech. 10 
(1930), 8. 81--92. — a Schulz, Interpolationsverfahren zur numerischen Integration gewöhnlichtr 
Differentialgleichungen. 2. ungew. Math. Mech. 12 (1932), S. 44—-59. — G. Schulz Fehlerabschätzun 
für das Störmersche Integrationsverfahren. Z. angew. Math. Mech. 14(1934), 8. 224—234.— W. To Ilmie x 
sr die I NET umeialn Verfahren zur Integration gewöhnlicher Differentialgleichungen. 
4. angew. Math. Mech. 938), 9. 83—90. — E. Lindelöf, Bemarques sur l’interrati eri N 
equations differentielles ordinaires. Acta soc. sei. Fennicae A 2 (1938), Nr 13 21S _- Le HTKe EZ 
2 .. ” In er h * u ” R . - ß ee ‚ x Bi: x x 
mühl, Beiträge zu den Interpolationsverfahren der numerischen Integration von Differentialgleichungen 
erster ur „weiter Ordnung. Z. angew. Matlı. Mech. 22 (1942), 8. 42—-55. —_ L.Collatz Natürliche Schritt 
weite bei numerischer Integrati r ilferentk: ichungssystemen. Zu a ch 2 ai j 
F re terischer Integration von Differentialgleichungssystemen. Z. angew. Math, Mech. 22 (1942), 

Ev r 
») Vgl. z.B. Fr. A. Willers, Methoden der praktischen Analysis Berli Leipzi 
DB. F Anı s. Berlin und Leipzig 1928, hier i 
besondere 8. 776, und G. Schulz, Formelsanmlung zur : Hr i er. 
£ var, li. \ AR els: zur praktischen Mathe u HÖSE 
Bd. 1110. Leipzig und Berlin 1937. 5 Pe u MEN ER 
JE. Steffensen, Interpolation, $. 29. Baltimore 1927. 
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der das Restglied von höherer Ordnung ist: 


3. Integrationen der Interpolationsformeln 


Integriert man bei der Newtonschen Interpolationsformel (6), (7) die Funktion f(x) über 
das Intervall von z, bis & (esist = x, + uh gesetzt), so a 


[1a de=hf N, (x(w)) du +h [Any ıda= fx) — BR 


| (13) 
38 +2u3 
ln kkens Track Pet] + Ru>+1 
} und bei n-facher Integration 
3 x © el 
| [He22= "ro — I emp) EZ 
„nm . nfäch RUE FE (14). 
(he Da Pn,e(u)v®efo+ Rn,» +1(8) 
o=u 
Dabei bedeutet f ein v-faches unbestimmtes Integral der Funktion f(x) = "f 
wf= [endfde (Bl ZEN A a TE (15) 
und es ist (mit Ro» +ı = R,;+ ı) gesetzt 
Ban mlBesıoride vH >> 5 u ee ee ee (16) 


und P,„,,(w) erhält man durch »-fache Integration 


u 


en) du--.du für o z1 


Ö N) n-tach (nl, 2, 8.) will 
1 
Pr, 10] (u) = nl au" 


I 
—— 
ARE 
——ı 


Pu,»(%) 


Einige P,,. für kleine Werte von n und o sind in der Tabelle angegeben. 


Polynome P, „(uw 


we 1 2 1 2 3 TR 2 3 4 EN 2 3 5 4 5 
n=]1l u. | gu 758w +2u°) 54 (Au +4wW+ u) nn (90 u + 110 u + 45 u +6 uw) 

RR ' N 2u: i -(20 W415 ul -3 u) { 60 u -+ 55 uw + 18 u +2 u) 
u=2 zu Tre u? + u!) | 360 u Ta | u 5 ut 
LE ! 3 AR ae DrbN | a 4 512] u 2 u? 
n=-3|;W |57 u a +2 uP) = (Wu +6w + uf) oogg (105 w +77 wW+21u° 2u) 
n=4 Er L — (3 u° + u$) 1: (14u® +7 u® -- u’) | 
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fix, + uh) an rap er 4) va, + 
N )) cher (11)N 
PATIENT 
(u —1) are ”) 
mit Ry= ae (£) ae, ' rege (12). 


ET FEN 


a 
ann 


ge are < rt vo n 


Collatz, Bir u nun Tg 


Mit Einsetzen Ber oberen Grenze x= x, entsteht N Be ER z = 


[f- fdz---da= "fln) — ER v®fo Hase) u 
Core Li n-fach , 


Dabei bedeuten die Koeffizienten Pn, q “ Zahlen 


me = - ft De ') du: du . wm. ” . . “ a9. 
fach 

Die ersten Zahlen fn,, sind 
Zahlenß,, 


In (18) läßt sich nach (7), (17) das Restglied abschätzen durch 
IBn,o+1l)| SPn,otı MP Hr HF Alma nenne (20). 
Setzt man als obere Grenze = z_,, so entsteht 


© mfach m talk 


[ 5 f = [ fs: de da = fie a Zenrn 4 


= (tree Diyu are fo + Rnatı (a) | 
Dabei ist ) 
; —1u U 
nem Se CHEN) ea [F-[da 
VEREIN n-fach Aue (22). 


= Fe An) 
Die ersten Zahlen yn, , sind 


Zahleny,, 


Natürlich kann man diese Zahlen y„,. auch durch Integration der Newtonschen Formel (8) 
mit aufsteigenden Differenzen A®f, gewinnen. 


& ® © 


p 
[rem 0 Mah-r Dmdriots 
e=0 


Bo fe % n-fach 


Er R Seas Fe an Stelle von &% in (13) bis 19) verwendet. Man erhält bei m-facher Integration 


REN. ndach | a 
a ER 2 y e nr WrhRinn Le | 
’ PO O, EUR). 
Dabei entstehen die ee P#,, durch n-fache Integration: art Bu: 
S . \ Y # 
> ER ek ka du--- du füro>i 3 { * j ve 
€ SE N RER etach % ha =ld..).-..B 
E \ - I . . } . e 
a ; P,oW)= u tl? | 2 ’ 
Die Den Tafel gibt die Polynome PP, . für einige Werte von n und e: i RUHE 5 
ar Pr, «m. re AR 
| 0=0 ; | oe=1 | R | o=2 BR: 
| EEE DEN N ER | RE Ur ER 2 
| | ’ 2 en he 
n=2 DR 1)? ige _9 ee) u—]) . | 
® 2 6 24 Ei 
n=3 ur Wu) We + u 0w—10u—3) 
n=4 Zu+ ı)* WR — 0159 u + 35 —10#— 1502 —9u—2) 
ger + 4 +4W—]) En 45 u! + 110 u? + 90 u — 19) 
n=2 Zu Be + u + 20#—15u—7) aa @ut + 18u5 + Su + 000° — 38 u— 1m) 


A 


sl ee ru . 
Außer ee Koeffizienten werden noch andere Größen ß7, . gebraucht. Dazu wird als ie 


5“ 


a Ge Ua da = 62 Ya | -3 enfa „“ al! | & ' | x 


gas 
n=3 70 r6W+ 10m 15 W— 14u—4) 


E 4t 7 6 5_ IE an & 
n=4 z00 (" +7u +l4u 35 u? - 49 u? -28u-6) 


Für 2 = Ro geht (24) über in 


P 2 z n—1 
f ER j fix. -de=Mf(x,) Ran 
ER | “_ı1 n-fach- Dat Dune (26). 
; » 
DR, 07° fo + Ru,» + ı(%) 
e=0 


Dabei bedeutet mit Ru,,+1ı= Rp» +1 
E7 
Bar a Me re re. 


lf- Be )o u 
—ı -i n-fach 


"Durch Einnetken der Definitionen beweist man | ER NEE 
Pn,e + Pi o+1 = Pa, o+1 KeulLand we ET 4 9. ns 


Wegen ßn,o= Pi,o = erhält man BERSBEN ION der Gleichungen (29) für =, 1,..„p—1 x 


2 Pe> —— Pa,» ” ” ” ” ” ” ” ” * a We Zu ” (30). 
a das Restglied Rt %,0+1 (20) folgt Sn aus (7), (27), (28) | 


IRkorılao)l SIAY,o lH H Er De nn Bi. 


Schließlich Klare wir bei mehrfacher Integration der Stirlingschen Interpolationsformel (10) 
mit 


Ba= 1 +uh) +f{m —uhl]- ...... = .".(82) 


mit = z, als oberer Grenze 
u 2% % 


| pr [| Ba2. a2 de vr Sm £p2ef, + R%,3a,) 


#o 0 To n-Tach e=u 


1 BEE 
[m f(@o +uh) + fa — uh)] — 2 Ba" oz) für gerades n 


n—3 


(33). 


"2 2A2e+1 

[fa +uh) — f(x, — uh)] a Bermi 1-20) f(x.) für ungerades n 
Dabei bedeutet 

Rt sea tre 2.) 


Ry= — Ri nach (11) 
Die ersten Zahlen 87*, sind 


N a 


Zahlen} 5 


PR EN Abi 
Fu 


i Es stehen im wesentlichen die gleichen Hilfsmit 
Ordnung zur Verfügung: 


I, Verwendung eines anderen Näherungsverfahrens. Auch hier wird man in erster Linie 


Die ersten Zahlen rf%, sind 


7“ Tl x# 
Zahlenry*, 


4. Berechnung des Anfangsstückes 


tel wie für Differentialgleichungen erster 5 


das Runge-Kutta-Verfahren heranziehen®), und zwar wegen der beim Anfangsstück besonders 


- wichtigen Genauigkeit am besten mit halber Schritiweise oder Viertelschrittweite. 


II. Benutzung der Taylorschen Reihe für y(x) und die Ableitungen. 


III. Benutzung von Quadraturformeln. Man kann hierzu die Formeln (25), (33) verwenden ; 


p Er 
TO +HAr—n I ya-mertfo (37). 


eu 


| h 
am = ym) don yantD 


Für eine Rechnung mit zweiten bzw. dritten Differenzen benötigt man noch Formeln für y(m) 


bzw. auch für yi®. 


hr —m—1 


h2 
BE 0 A Ee NEE NT TEN 
ae en 


n—m—3 
2 
h2e+l 
m 2 1) 
52 m +22, Bann 
+17 n—m—2 > 


au hr, 


Diese Formeln verwendet man als Integrationsformeln, indem man sich zunächst Rohwerte 
für y{m, y(m, y(m) (etwa grob mit der Taylorschen Reihe bis zu drei Gliedern mit y ein- 
schließlich) verschafft, dann f,, f, fs und die Differenzen ausrechnet und diese Werte auf den 
rechten Seiten einsetzt, um die auf der linken Seite stehenden Werte dadurch in nächster 
Näherung zu erhalten. Wir beschreiben nun mehrere Verfahren der fortlaufenden Rechnung. 


Es seien bis für die Stellen x, (füroe=1,...., r) Näherungen Y, Yo - , 
exakten Werte y(zo), Y'(Xe) - - »» y®=2 (x,) bekannt, und es sollen Näherungen %r--1, Yr+1, 


5) Formeln für Differentialgleichungen n-ter Ordnung wurden aufgestellt von R.Zurmühl, Runge- 
Kutta-Verfahren zur numerischen Integration von Differentialgleichungen »-ter Ordnung. Z. angew. 


Math. Mech. 28 (1948), 8. 


hre 


2 
m 


e=u 


+2 m (Bir, + Barı 7? +ı) 


\ n ungerade 
| n gerade 


5. Extrapolationsverfahren 


173—182. 


TR 


9°» für die 


206 Collatz, Differenzenverfahren zur numerischen Integration BR En ri 


ou, Y%jD ermittelt werden. Nach (18) gilt für die exakte Lösungsfunktion, wenn man dort n 


durch n— m, ferner x, durch x, ersetzt und y®) = f benutzt, 


: n—m—1 m ; - ir 
ya Ft Hin! BumoPtRlz) + Ra-matı : 89) 
} „= gi Be 
wo a ya), ya)...» ydla)=Fla) ......: Pe, 
- gesetzt ist, und für das Restglied die Abschätzung (20) gilt. Läßt man das Restglied R fort, 
so erhält man die Falknersche Extrapolationsformel®) für Näherungswerte y(, : ne 
n—m—1 r $ 
ym, = u yirtn mon ee, N a a 
v=0 vl 0e=0 


Für n=1 ist die Adamssche Extrapolationsformel für Differentialgleichungen erster Ordnung 
als Spezialfall enthalten, für n=2 ergibt sich 


Yyr+1ı= Yy+hyr 4m an+4 7% +z 7’ rs] 
’ ’ | 1 5 3 (42). 
Yr+1= Yr + nttrntgrn tern + 

Allgemein sind die Koeffizienten ß„,, durch (19) gegeben. 


Es lassen sich noch auf mannigfache andere Arten Extrapolationsformeln aufstellen. Man 
kann z. B. durch Kombinationen Formeln aufstellen, in denen nur y- und f-Werte vorkommen, 


in denen also alle Zwischenableitungen eliminiert sind. (Solche Formeln sind besonders geeignet, t 
wenn die Differentialgleichung die Form hat y”=f(z, y).) Zur Herleitung schreiben wir die 
Gl. (14) für die Argumente & = 2%, 2_,, 24... £—_n+ı an und multiplizieren diese Gleichungen a 


einzeln mit den Faktoren ) (2) E= Pal () (8). N (jr () multiplizieren also die 


N 


Gleichung für <= x_, mit (—1)e+1 Pe 7) Addition dieser Gleichungen ergibt mit "f= y 


vi) = (1) ytr) +) van HN Yen = Prye)= Ar Din Pe +R 3) 
mit 


n—1 
N = n zT +1 n n SE rn R PERT 
En. Pu.) HZ, (-Ntt( 7) Pan (44) 
und 
< n—1 
R= Rus+ı la) +2, ET, 4) Be a Te (45). 


Die ersten Zahlen &,,, sind: 
Zahlen,, 


le= ole=1le=2]e=3]e=2]o=5 


? 1 5 3 | 3ı| 
"II I Nm |:T 170 (88 
1 1 Pig ls 
—=2 Et kn er ee 
9 L 253 93 
1 1 
n=38 1 re nn 
g [08 ° | 0 
1 1 


n=4 N ee Pr et Dee 
.6 720 

Die Formel (43) für Differentialgleichungen n-ter Ordnung enthält speziell bei Fortlassen 
der Restglieder für n=1 die Adamssche Extrapolationsformel und für n=2 die Störmersche 
Extrapolationsformel. 

Bei Differentialgleichungen dritter Ordnung fallen wegen £,,=233—=0 die Glieder mit pr 2f 
und p?f heraus: 


YHı=3y— 3 Yyr—ı + Yr—2 + (rt +). n. (46). 


©) V. M. Falkner, A method of numerical soluti f di i i i 
Sl Weiner solution of differential eunfione: Philos. Mag. 7. ser., 21 


NERS 
Eh 


De a) Tann Dr de Ad a 1 N am 2 


a 


‚höhere Differenzen benötigt). 


u Mi 


man hier mach‘ ehe mit v ab, so hat man eine Kar: ur Kind a 


ui Em. Man braucht nur die ersten Differenzen mitzuführen (sofern man nicht wegen Auf- ER 


tretens von y’ oder y’’ in der De und damit megeh, der Formeln für y oder 6. er 


Auch bei Differentialgleichungen vierter Ordnung ist die De Abbrechen ve; den Gliede N 
4 Waren 

ae‘ Hay ya +h er u 720 ve Er (a7) 
Enistchinde Formel sehr einfach. j ET 
6. .Tnterpolationsverahren R 

Wieder seien für die Stellen ,(e=1, ..., r) Näherungen yo für die a! Werte 
ym(z.) fürv=0,1,..., n—1) bekannt, DM Näherungen y) , "sollen ermittelt werden. 


Nach (26) gilt für die ke Lösungsfunktion, wenn man ni N durch a %, durch Gr +1 


ersetzt und y'9=f benutzt: | ER 5; = ” 
m 1 mv u g i BE ad 
yım) (& +1) =: sh yinın (a) mn en ve F(&r+ 1) + Ri _mp+ 1 . (48), N e: 


wobei ı wieder die Abkürzungen (40) für F, (28) für Pä—m,e benutzt sind und für das Restgied 
die Abschätzung (31) gilt. Bei Fortlassen des Restgliedes erhält man die Interpolationformel 


für Sr ae YD:- ER BR: 
Dan Be; 
ur, ee Bunt v) > Pr-m, eve fı+ ee BON 0 


Für Ber. m=0 ist die Adamssche Aue ER für Differentialgleichungen erster 
Ordnung als Spezialfall enthalten. 


In (49) treten die Unbekannten y(”), auf beiden Seiten auf, man verwendet daher die au Be; 
Gleichung zur iterativen Berechnung, : man die »-te Näherung y{mP1 auf der rechten 5 Tu 
Seite einsetzt und aarane die (» +1)-te Näherung berechnet: 

un w+1] Au . yım+ ») > rm A ER ER RN (50). 
v0 
Die Näherungsstufe » ist dabei zur Unterscheidung von ee Ableitungen durch eckige Klammern BR 
[v] gekennzeichnet. F 

Gelegentlich verwendet man die Gl. (49) auch in der Form, bei der die Differenzen Fef Bi. 

durch die Funktionswerte ausgedrückt sind. Die Gleichung lautet dann | k 
n—m—ı1 hv , Br 
\ ym,= = ge — yim+r » + hm P Piom an Srrie Er  AHIETN E: 


mit den neuen Koeffizienten 
r n | 
BVZ ße (%) RESTE ..(82). 


Wie bei den Differentialgleichungen erster Grasung kann man auch hier noch andere Inter- 
polationsformeln aufstellen. So gilt nach (33) 


T=-1 120 
ym) (41) = — ym(—) +2% Boll 1 y(m +20) («,) 
e=0 
“ £ . 
L2m—m ZpEn er?? F(%,+.) +2 Rx m,» für gerades n — m (53) 
oe=0 j ® . . . 
5 er 
ym) (241) = yM (0 —ı) Bir: ee er) 


Hamm} BR-m, or? F(&,+.) +2 Ri m,» für ungerades n— m 


Für die Restglieder gilt En Abschätzung (35); läßt man sie fort, so erhält man genau wie oben 
Gleichungen für die Näherungswerte y/”),. Das so aufstellbare Interpolationsverfahren ist das 


+amon2 Pme BEN 


' Gewöhnlich wird man die letzte Se beim Gliede p=2, also 2 1, abbrechen, da für -2 
bereits außer der Unbekannten +1 noch die Unbekannte en ine na en 
Durchführung des Interpolationsverfahrens umständlich 


Es seien noch einmal die Formeln (54) ausführlicher: map. (Ur enden © DZ 


chungen reichen aus für Differentialgleichungen bis zur vierten Ordnung 


yazd = ya +1(24, Bars, ke % = 
| Da yera gene ans +? (r+ a pr) I 2 (55). x = 
ya 3) zn +2 hyn-d +5, *m FR EN 
ee: Fr ger eg yr-d — UJAzENe + A: yr2 + (2, 15h tr: 3607 fı = 
N | %. Konvergenz der Iteration bei der fortlaufenden Rechnung 


Zur Konvergenzuntersuchung werden die Interpolationsformeln in der auf die Funklions- 
werte umgeschriebenen Form entsprechend (51) verwendet. Schreibt man die Gleichung für » 
und »—1 an, so ergibt Subtraktion für den Unterschied die Beziehung 

öl) Lv] — ylm) HN] — lm) 1 


om) WI] — hr —m Ban ll, En Ze 2] 
Dabei sind die Zahlen ßf,o,„ mach (52) und nach (30) 
Bao,» = e= Ba,» 
Die rechte Seite wird mit Hilfe der Lipschitzbedingung (3) abgeschätzt 
n—1 - 
N N 1 ol, 2 k,[sab-uU (m=0,1,..,n—1) ... . 157). 


Multipliziertt man diese n- RR ai Ko Ky ..., An, so erhält man bei Addition 
für die Größe. 


n—1 ö 

N= I K,|lsam| , .... SER (58) 

e=0 i 7 

die Ungleichung | 

VS OR NE ee (59). | 

mit Ä 

C= A" |Bm,p| Kr + ro |Pa-ı, »] K, ee. +h|Pßı.»| Ko Be a gern (60). \ 

oo 4 

Nach bekannten Sätzen ist O<1 hinreichend für die Konvergenz von I «{) ei da wir die 
v=0 

K, als positiv vorausselzen können, auch für die absolute Konvergenz von 3 swb) beim 

Iterationsverf ahren, v0 


Die hinreichende Konvergenzbedingung lautet somit 


A 1 N 1 n 1 . 
 fürp=1: zh Ku-ı RE h?2 Kung +.% Ks +’ +2. <1, 


für p=2: Kn—ı + Kuh zii Ans pres + in.) <1 


”E. Lind elöf n.a. 0. gibt bei Differentialgleichungen erster it i - 
kichtigung dee GH ER g gen erster und zweiter Ordnung Arten der Mitberück 


ii ‚pnehme 
‚legung für das Terlanren der ) bei A je mit 
pr? er E Führt Reheahen, auf vr, | SM ER ee 

) . j 5 a 2, ho us 1 ent ER N | = 
oder Ausführlich 


| VE me | n | | wi W 
N ER 4m a im NS A +2M 7: HRa<ı. LE RPe „00. Me 


Re ei ‚obigen Fällen für py=1 Ai ?=2 sind die Koeffizienten abermals verkleinert i 
a darf h noch etwas Ba als oben gewählt werden. 17 
1 8 Prinzip der Fehlerabsehätzung für die fortlaufende ee | 
ce Die Fehlerabschätzungen können für die verschiedenen Arten von Extra- und Inier- 
polalionsverfahren in derselben Weise durchgeführt werden. Wir beschränken uns auf die Ab- 
‚schätzung bei der Interpolationsformel (49). Durch Subtraktion der Gleichungen (48) und 49, 
. wobei die Differenzen 7 © f entsprechend (51) durch die Funklionswerte ausgedrückt sind, erhält 
man für die Fehler 


13 En a ee 
die Gleichungen Dun AR 
Nn—m—1 hv ? ? » R a i Y h g R a 
EN E } ih, en ramn2 Bu—m, oo» Ar sh a3 1 —.o) ))— Ra—m,p+1 (63). 
a Dabei ist die Abkörkung (40) verwendet. Die Differenz der Funktionswerle kann mit Hilfe‘ der, 9a 
2 N ee 8) abgeschätzt werden: 
2 Er. $ jem , nz K, WI ler} DEN mit ART 
x BR : 2. (64). 
(mn= 2, ar R jem+n] > So K, | Br m,.o, ‚leer, +lBr_mp+1]l 
y= 


In (64) haben wir ein System von » ce für n Fehler En „ aus denen man sie ab- 


schälzen kann, wenn man die rechten Seiten em, d.h. Schranken für die früheren Fehler zu, 
em) ,... kennt. Hierbei wird ein Hilfssatz von G. Schulz benutzt!), da in (64) 1 — hı—m K, 
ee mr| > 0 (nach der bereits in der vorigen Nummer nötigen Voraussetzung C <1 in (60)) 
und alle anderen Koeffizienten auf der linken Seite negativ sind. 
Für obere Schranken Y) der Fehler (Y® > |e(®]) erhält man somit ein Gleichungs- 
system, indem man in (64) Y{® statt e(”) und anstatt des Zeichens < ein Gleichheitszeichen 
schreibt. 
Damit ist eine rekursive Fehlerabschätzung durchführbar. 
. Man kann auch wie bei den Differentialgleichungen erster und zweiter Ordnung inde- 
i pendente Fehlerabschätzungen durchführen und Schranken 


ii) SYymt Cm? 


aufstellen mit y„ und c„ als Konstanten. Die für, das Wachstum dieser Schranken entscheidende 
Größe z ist dabei Wurzel der Gleichung D= 0 mit 


2PI—WK Bl) 7 —a! n Be si eu Nele ei ei m MN Kr 
—h"1 K,D,(2) FE a N 2 ih, Bl) — m Fe Tale Ka-ı P; 
S E h"2 K,d;(2) — N K,D;(2) 2 — zP1 — MR,D,(2).- — zn ee — N"? Kaı Ds; 
—h Kuda lt) ER Due) SER,dnil).... Det Kenn 

wobei. zur Abkürzung © (2 —. | Br m 5iel 2? 7@ Besetzt ist. 


8) G. Schulz, Z. angew. Math. Mech. 12 (1932), S. 553—554. 
Eingegangen: 30. 10. 1948. 
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Numerische RER einer Differentialgleichung der 
Wellenmechanik | 


Von Th. Schade in Aachen 


Die vorbiegende Arbeit behandelt die numerische Auswertung einer speziellen Schrödingergleichung des 
Mehrkörperproblems der Alomphysik unter Zugrundelegung gewisser analylischer Aussagen vor allem im 
an Verhalten der Sa der Different ichung. Damit kann das bisherige Auswertungs- 

ahren erheblich gekürzt werden. \ | 2 A, 
ie m I. und II. Heil wird die DET AERO. RISSE ee rs Eigenwert m ar. 

obei auch an Hand von Störungsrechnungen quamlilative Aussagen über Eigen derungen. ' 
Ei des Potentials und anderer Funktionen der Differentialgleichung gemacht werden, die für die Iterations- 
rechnung wesentlich sind. 


Der III. Teil zeigt dann an einem Beispiel die numerischen Auswertung und die Genauigkeit des Ver- . 


jahrens: 


The present r deals with Ihe numerical treatment of a special8chrödinger equation ooncerning 
the Eat ofa Ale of particles in atomic physics that is based on cerlain ya stalemenis primarily 
on the asymptolic Pia of the solution of the diffential equalion. Thus the method of Ihe Ireatmenti 
hitherlo used is shortened. t 

Im part I and II, the differential equation of the second order is treated as eigenvalue problem. Quanti- 
lative siatements of Ihe alteration of Ihe eigenvalues with the variation of Ihe potential and other functions 
of ihe differential equalion are given that are imporlant to Ihe process of ileralion. . 

Part III illustrates the numerical treatment and Ihe accuracy of ihe melhod by an example. 

B 9Tof pabore N3NATaeTCA METON YHCHEHHOTO PeIIEeHHA ONHOTO CHENHANBHOTO YPABHeHHA 
Ilpenusrepa, oTHocameroca K IpoÖneMe CHCTEeMEI CBODOAHEIX WACTHL 'aTOMHOH PuauKH. 
ITo pemeHne UCXOAHT U3 HEKOTOPEIX AHAJIHTUYECKHX YTBEPKNeHHÄ, KACAMIMHXCH B IEPByW 
oyepenb ACHMITOTUYECKUX CBOÄCTB pemeHHuä nudepeHuualbHbIX YpaBHeHHl. ITEM NyTeM 
SHAYUTENBHO YOPOMAETCH IPe3KHHH METON PeIueHngn. > 

B I u Il yacru pacematrpusaerca IPpoÖTeMa COÖCTBEHHEHX PemeHHuä AHdepeHlHansHsX 
ypasHeun& BTOporo nopanka. ‚KonuyecTBeHHbä AHaJIUs H3MeHeHHA CODCTBEHHHX 3HAYEHHÄ 
IPuU WSMeHeHHU NOTEHNHAAIA U APyTux hyHEuuH nmpepenHmmansHorO VPABHeHHRA B BHAYHTENL- 
HOHU Mepe yupomaer mpolecc HTEepanun. \ 

B Ill yacru unmoceTpHupyEmTcA Ha IpHMepe xOA YHCHEHHOTO PeIIeHHA H-TOYHOCTb MEeTONA. 


Vorbemerkung 

Die in der Arbeit durchgeführte Behandlung der Differentialgleichung 2. Ordnung ist auf 
Grund des vorgelegten physikalischen Problems entstanden. Es ist aber hier weitgehendst in 
den Hintergrund getreten und der mathematischen Durchführung das Hauptaugenmerk ge- 
schenkt. Die abgeleiteten Formeln sind möglichst allgemein gehalten, um gegebenenfalls ähnliche 
mathematische Probleme in analoger Weise der numerischen Auswertung zugänglich zu machen. 
Die numerischen Rechnungen im Teil III der Arbeit sind mit üblichen elektrischen Rechen- 
maschinen im Max-Planck-Institut für Physik in Göttingen durchgeführt worden, wo infolge 
der erfolgreichen Versuchsrechnungen auch weitere laufende Rechnungen des gleichen physika- 
lischen Problems mit der hier beschriebenen neuen Methode bearbeitet werden. 


I. Einleitung 

Das quantenmechanische Mehrkörperproblem der Elektronen einer Atomhülle wird mathe- 
matisch durch die bekannte Schrödingergleichung für die sog. Wellenfunktion formuliert. Unter 
der Voraussetzung eines kugelsymmetrischen Feldes läßt sich eine Trennung der Variablen in 
räumlichen Polarkoordinaten durchführen, wobei dann noch das einzige Problem in der Be- 
stimmung des nur vom Abstand r vom Mittelpunkt abhängigen Teilfaktors der Wellenfunktion 
besteht (vgl. Hartree: Proc. Cambr. Phil. Soc. 24, S. 89). Die Schrödingersche Wellengleichung 
reduziert sich hierdurch auf 


P"(r) +|20- .— dan P(r) +F(r)=0 


r? 


mit V= —v als Potential des Bohrschen Atommodells und E= 5 als Gesamtenergie. / wird 


als Nebenquantenzahl bezeichnet. P ist der noch mit r multiplizierte Radialfaktor. Im Falle 
des Wasserstoffatoms mit einer s-Schale (2=0) erhält man eine einzige Differentialgleichung, 
die exakt mit Besselschen Funktionen gelöst werden kann. Bei schwereren Atomen erhält man 
für jede Schale n, 2 (Haupt- und Nebenquantenzahl) je eine Differentialgleichung, die aber durch 
das Potential miteinander gekoppelt sind und die somit ein System von Differentialgleichungen 
darstellen. Über die Aufstellung der Differentialgleichungssysteme in angeregten und nicht 
angeregten Zuständen der Elektronen soll im Rahmen dieser Arbeit nicht gesprochen werden 
(vgl. dazu die entsprechenden Arbeiten von Hartree, Fock u. a, Z. Physik 61 [1930] 126). 


ee 
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en Le, Ditferentialgleichung der Wellenmechanil 
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und Mehrelektronenspektren der Atomphysik lautet im Prinzip | 


ER | Birmi) + 7 HD 


‚wobei bei Mehrelektronenspektren ein entsprechendes System von Ditferentialgleichungen von 


diesem Typus zu befriedigen ist. 


T(r) und F(r) selbst enthalten bei Mehrelektronenspektren die Wellenfunktion Py(n In), | £ 


als Lösungen des Differentialgleichungssystems in Form von Integralausdrücken dieser Funktion, 
so daß die richtige Wellenfunktion iteriert werden muß. Die Haupt- und Nebenquantenzahlen n 


und / gelten als Parameter der Wellenfunktionen. Sie bestimmen sich aus dem Zustand des 
betrachteten Elements. P(nl/r) soll’ im weiteren Verlauf der Rechnung durch P ersetzt werden. 


‚Py soll normiert sein, d.h. 


VE STEINER EN BE 
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Für die homogene Differentialgleichung (1,1) ist diese Bedingung unproblematisch, da jedes Ba 


‚Vielfache einer beliebigen Lösung auch eine Lösung der Differentialgleichung ist. 
‘ eist der Eigenwert der Differentialgleichung, dessen Bestimmung das eigentliche Haupt- 
problem bedeutet. Für Py gelten die Randbedingungen \ 


fürr=0 IN, 
Pr=0| fürr=o@. 


Die erste Bedingung für r—0 läßt sich dadurch erfüllen, daß für Px in der Umgebung des Null- 


punktes eine Potenzreihe angesetzt wird, die mit a, r!! beginnt. Da 7(r) und F(r) mit r— oo 


nach Null gehen, würde die dann übrig bleibende Differentialgleichung die Lösungen e 
arten, zulassen, wovon die erste die gesuchte wäre. 


°’ und 


Zur Integration der Differentialgleichung kann man also mit r=(0, Innenintegration ge- 


nannt, und r—, Außenintegration genannt, beginnen und beide numerischen Lösungen unter 
richtiger Wahl des Eigenwertes und Are 

unter Beachtung der Regularität an 
der Flickstelle zusammensetzen, was 
einen erheblichen Rechenaufwand er- 
fordert. Diese gebräuchliche Methode 
soll nicht im einzelnen beschrieben 
werden. 

Zur Ermittlung des Py bedient 
man sich eines Iterationsverfahrens. 
Aus einem geschätzten P, das in T(r) 
und F(r) eingeht und einem geschätz- 
ten e berechnet man die Wellenfunk- 
tion .P aus der Differentialgleichung 
unter Beachtung von Regularität (Ei- 
genwertvariation) und Normierung 
(a,-Variation), die zunächst nicht 
mit der in 7 und F eingesetzten über- 
einstimmen werden. Durch Mittlung 
dieser Funktionen und Änderung von e 
lassen sich neue Näherungsfunktionen bestimmen und damit auch iterativ die richtige Lösung 
und das dazugehörige efinden. Die Regularitätsforderung in jedem der Iterationsschritte wird 
unter Variation von e dadurch erreicht, daß die logarithmischen Ableitungen an der Flickstelle 
für Außen- und Innenintegration in Übereinstimmung gebracht werden. Da die Flickstelle in 
den Bereich des Wendepunktes gelegt wird, kann man statt der Ableitungen den Differenzen- 


quotienten setzen und erhält: 
B2, > LlaR).. 
h\P en h\P /außen ; 


rn ren 


Bild 1 


14* 


-e) Py (ni) +Fn=0 3 2 Bye Ra 


2 

Em SES sei hier nur noch erwähnt, daß die Bestimmung der Wellenfunktionen für die Berechnung N 2 R 

gewissen Größen der Atomphysik wie Gesamtenergie des Systems oder Oszillatorenstärken not- Br 
wendig sind, soweit sie nicht mehr der Beobachtung entnommen werden können. u 

_ Die vorgelegte Differentialgleichung für die Wellenfunktion P zur Berechnung von Ein- 


a 
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bei r, das Argument eines Funktionswertes in der Umgebung des letzten Wendepunkter 
ehentet (die eat Abszisse wird mit r=r, bezeichnet) und A die Intervallänge von r. Der 
letzte Wendepunkt ist in der Differentialgleichung durch 


AU 
1m -e=0 
gekennzeichnet, wenn man die dort schon kleine Größe von F nicht berücksichligt. Alle weiteren 
Wendepunkte für r<r, sind durch P=0 ausgezeichnet. 


Eine Innenintegration über r=r, hinaus würde in die e/er-T.ösung der asymptolischen 
Differentialgleichung führen, da der geschätzte oder jeweils iterierte Wert von & nichl genau 
dem richtigen Eigenwert entspricht und die Lösung für große r damit sehr instabil wird. Die 
vorliegende Arbeit soll nun einen Weg zeigen, die Funktion für die Außenintegralion von der 
Umgebung des Wendepunktes r=r, an zu bestimmen, so daß zumindestens die Ordinaten- 
übereinstimmung der Wellenfunktion mit der Innenintegration erreicht werden kann, da man 
über die Größenordnung der Anfangswerte für eine Außenintegration nur sehr wenig aussagen 
kann, wenn man nicht zu sehr großen Argumentwerten übergeht und damit mit sehr kleinen 
Funktionswerten rechnen muß. Der Zusammensetzungsprozeß wird dadurch sehr mühsam. Im 
numerischen Teil der Arbeit wird außerdem noch gezeigt, daß es durchaus möglich ist, daß die 
Wellenfunktion mit großen Argumentwerten ihr Vorzeichen noch ändern kann. Das bedeutet 
eine weitere Schwierigkeit in der Schätzung solcher Anfangswerte zur numerischen Integration, , 

Es soll unter Benutzung von asymptotischen Darstellungen für F und 7, die vom letzten 
Wendepunkt an gelten, die numerische Außenintegration durch eine exakte Lösung der asymp- 
totischen Differentialgleichung ersetzt werden und aus ihr eine zumindestens für Näherungs- 
lösungen brauchbare Funktion ermittelt werden. Wie man später sehen wird, können für F 
auch die numerischen Werte benutzt werden, was auch zu genaueren Ergebnissen führt; der 
Vollständigkeit halber sollen aber beide Möglichkeiten diskutiert werden. 

Im numerischen Teil wird im einzelnen die Praxis der numerischen Auswertung erläutert 
und an Beispielen die Genauigkeit der entwickelten Methode gezeigt. 

Zum Abschluß sei hier noch als Beispiel das Differentialgleichungssystem des Ca*+ an- 
gegeben und das Resultat in der folgenden Abbildung skizziert. (Die Anzahl der Nullstellen der j 
Wellenfunktion ist a—1.) (D.R. Hartree und W.Hartree, Proc. of Royal Soc. Vol. 164 re 
(1938) 167.) 
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Y£(o,'ß/r) = j Pyta/r,) PyxtPlr,) (*) dr,. 


Es.sei noch bemerkt, daß das positive Vorzeichen im jeweils letzten Term der Differential- 
gleichung den sog. Triplettzustand und das negative Vorzeichen den Singulettzustand bezeichnet. 

Die Wellenfunktionen des Rumpfes P(1s/r), P(2s/xz), P(3s/r), P(2pjr), P(3p/r) sind 
numerisch vorgegeben. 

Bild 1 zeigt .den Verlauf der Wellenfunklionen P(4s/r) und P(4p/r), wobei der- 
Singulettzustand durch die ausgezogene und der Triplettzustand durch die gebrochene Linie 
dargestellt wird. 

Die Eigenwerle sind: 


Triplett: &4,,4s = 0,518 Singulett: &4,, 4, = 0,505 
&4p,4p9 = 0,306 F4p,4p = 0,172. 
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Form: 


 undetwar=3r,. Essoll später im numerischen Teil noch gezeigt werden, welches der Bedeu Eu 
Wert für 5 ist. 


an Diffretilgleichung ek: EN | 


2 x 2 +(r — Il a ” 1 P Re . pr a np HE = 


Wie schon in der a erwähnt, wird die Innenintegration mit r=0 unter. ad N 
% 


ech 


Ks: einer beliebigen (-+1)-ten Ableitung der F bektion P begonnen und bis zum Meder 1 
BE durchgeführt. 


Für die anschließende Lösung wird eine approsimalive Darstellung f für H benizt in der 


T=2 +5 ER ENLERLE Bo 


Die Konstahfeh a und b bestimmt. man aus (2,2) durch Be der Fiıktiönswertz für BR, 


Se man weiterhin be &, so lautet a asymptotische Wert für T, 
re SLIE4-1) “ 


$ w KIN y2 


. Der Ansatz für die homogene Lösung der Ditferentialgleichung (2,1) lautet dann 


Pearl tar BAHR) 22222222. 2A). 
Der Koeffizientenvergleich liefert nach Einsetzen in 2 1) FE: 


K=—_ rl 
2 Ye RR 


Kik—ı) + +1) © (K4 ei 


= K+D)(K—1—1)]-{[(K +4—1)(K—1—2 
a. ap EINS iR +1) ( ) 


on -h! 


Für die a, läßt sich weiterhin eine el lonsurmel ableiten 


RE KTNKHK-T N Rn HH KH) 


n-1 
een En K— KR ENE D e  A ‚6). 
a RR (2,6) 


Damit lautet die homogene Lösung 


Pair (1 SAN: = I ER 2,7. 


; \ u=1 
Diese Reihe ist semikonvergent. Zur Darstellung einer numerisch brauchbaren Lösung kann 
man aber den Teil der Reihe benutzen, für den die a, mit wachsendem » kleiner werden und den 
Fehler berechnen, der entsteht, wenn man die Reihe bei n=N abbricht. Man wird die Reihe 
bei dem Glied abbrechen, bei welchem |a„| ein Minimum ist, d.h. 


Mila u ee er ER AR T WM 
Damit würde die » Differentialgleichung lauten 
P" +T,(r)P=Ki(r) 
mit 
Kir) =ay ee sK—2(K— N+H)(K— N—-I—-N)...... (2,9). 
Zur Berechnung des Fehlers X(r) der rechten Seite der ursprünglichen Differentialgleichung 


genügt es den maximalen Fehler für r=r, festzustellen, da die Funktion X (r) mit wachsendem r 


sehr rasch abnimmt. 

Ist K (r) von der gleichen Größenordnung wie die numerische Berechnung der Funktion F(r), 
so kann man ohne weiteres X e) gleich Null setzen. Dieses wird zumeist der Fall sein. Anderen- 
falls kann man u in einer im Abschn. 5 durchgeführten Störungsrechnung die Änderung des 


} E \ der semikonvergenten Reihe mit (2,7) Um s 
gleichung nicht mehr derart möglich ist, ee 
{A f 
Damit folgt für die asymptotische Lösung der homoge 
| nA. ET 
Pi=Gelrrk(i+ Dan) VE 0 

3 n=1 mr - nr ne Bez 
wobei C, so gewählt werden kann, daß für r=r, der Ordinatenanschluß an die Integration von i 
innen gewährleistet ist. Die Lösung der inhomogenen Differentialgleichung P;, läßt sich dann 
bekannterweise in der Form darstellen rn 


Pa= mr] Pr| Zi Prrar . * en) Be: 


Die Integrationskonstanten bestimmen sich aus den Randbedingungen. Da die. Innenintegration 
nach dem üblichen Schema am besten numerisch von r=0 durchgeführt werden kann, interessiert 
auch hier nur die Bestimmung der Funktion für r Zr. | 


Da P, für ro mit eV’ und auch P;,für ro nach Null gehen soll, bestimmt sich Cs 
zu Null. Wird die weitere Integration von r=r,, an durchgeführt, so lautet die Lösung für P;, 


N | 
Er } 


wobei für die Integrationskonstante (, gilt 


Gr [rar NS ee 
Fi. 


e- Schreibt man noch für — [ FP,dr=(,, so lautet (2,12) 

> | "o $ je | 
4 0, + FPıdr | 
F Pi = 0,Pı vo Pr SA; RE, dr ee (2,14).. 
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Dt Es ist hierbei zunächst nicht notwendig, die Innenintegration für die homogene Gleichung durch- 


zuführen, da ein Faktor für den Ordinatenanschluß sich aus dem Integralausdruck herauskürzt 
und die Integrationskonstante CO, den Ordinatenanschluß an die inhomogene Lösung von innen 
gewährleistet. Aussagen über Regularität und Normierung sind damit.noch nicht gemacht. 


3. Näherungsweise Erfüllung der Regularität 


Ist für die Innenlösung P die Änderung von P bei einer Änderung von g um As, so heißt 
das für die Differentialgleichung 


P'+AP" +HT—e—A)(P+AP)+F=0......,..681) 
der mit P’+Nn—e)P+F=0 | 
APULTD EN AP—L- PA A PER 


Wird Ae: A,P als eine kleine Größe zweiter Ordnung vernachlässigt, so lautet die Differential- | 
gleichung für A, P 


AP AP- Ar Pet 


Bi (A.P= const P) so daß man erhält” 


"Integriert man nur ber E08 gemeinsamen Gültigkeitsbereich der Außen- und Innenlösung 
E NE nr are so en Er Er re A 


‘ Damit ergibt sich für ar Ordinatendifferenz der BRERIOSEnBE im gemeinsamen Gültigkeits- 
bereich W 


Dr DeBueriergeräng lautet damit - 
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ersten. beiden Summanden von 2 “ entsprechen der "homogenen Ditferentialgleichung | a 


& 2 A Pan con: P #0 [% fa Padr : Ei Eu Kaas A ne 0 


a x ee 


Ty 


4 Pr) = 4, ein Her, in [ran 
8 ET: a Papa & Een 2 a ES 


A, Pin) — AeP(n) = Ae:0g: Pr | 


KAP on = Pin) N Pin — Pen) 2er Pr). = Pfr) Pe) Zt de: "60° Pit (1) ie 


A Pausen) = A Pu cinnen) = P(r,)— Der Cg: Er (r) IR u (3, A). 


Rechnet man sich für die ersten beiden ae des Polynoms zwei AP mit Verschiedeneit & 
aus, so kann man durch lineare Interpolation die Größenordnung der Änderung von & bestimmen, 
ohne die Innenintegration selbst durchzuführen. Eine numerische Nachprüfung dieses Verfahrens 
folgt im numerischen Teil ger Arbeit. 


‘4. Näherungsweise on lag der Normierungsbedingung 
Pin ist der Normierungsbedingung unterworfen, die eine Variation des a, zur Folge hat. 
Für eine Iterationsrechnung ist es nicht nötig, diesen Wert sehr genau zu bestimmen. Eine gute 
Schätzung für die Bestimmung der Anfangswerte bei r=0 ergibt aber folgende Überlegung. Yfe 
Setzt man in (2,12) für das Doppelintegral D(r), so erhält man durch Quadrieren und a 
Integrieren von 0 bis oo unter der Voraussetzung, daß 'P;, normiert sein soll 4 


1 -/ P} (Dir) — O,)?dr. 


Beer rapea ri BD: | 4 
[ Par [Par Y 
0 v ö : 


Bedenkt man, daß die D(r) enthaltenden Integrale sehr klein sind, da D(r) im wesentlichen die 
Korrektur der homogenen Lösung durch einen inhomogenen Term bedeutet, so ist, wenn die 
homogene Lösung P, normiert ist, CO, ungefähr eins. Das bedeutet aber, daß die Anfangswerte 
für r=0 der normierten homogenen Lösung auch für die inhomogene Lösung einen Normierungs- 
faktor==1 ergeben. 


Das heißt 


5. Störungsrechnung 


Es ist oft erwünscht zu wissen, wie groß etwa die Eigenwertänderung sein wird, wenn das 
Potential T oder die Austauschglieder # um einen von r abhängigen Betrag geändert werden. 


‚Da die Wellenfunktion selbst keine singulären Stellen besitzt, ist diese Rechnung ohne weiteres 


durchzuführen. Ändert man T um AT, F um AF, so ändert eum Ae und Pum AP. Man erhält 
durch Einsetzen in die Differentialgleichung unter Berücksichtigung derselben und Vernach- 
lässigung von kleinen Größen zweiter Ordnung. 

d.,F, mPin + (Ts $) Je; F, rP;r - (A T— de) Pin 4- Az (51). 


216 Schade, Numerische Integration einerDifferentialgleichungder Wellenmechanik 74.20 
Ist &, der Eigenwert der zugehörigen homogenen Differentialgleichung 
i Ph +(T,—o) Pr=9. a Te A W 


so kann man auch für (5,1) schreiben 

Ar. Pin +1, — 6) Pn+la— 9 Aur,r Pa +AT— A) Pu +AF=V 
oder wenn man T,—e, aus (5,2) einführt | 
(Au. n Pin Pa— PiAee,n Pin’ + (&ı—e) Au, rn PinPa HAT— 4e) PinPn + AFPi=0 6,3). 
Integriert man (5,3). von 0 bis @, so folgt die Beziehung 


1 = 7 f | 
de= ——— Ik T P;„Pıdr +[ AFPıdr + (4-2) JAur, nPinPudr]|. s 
v 


[Pin Pıdr y 
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Da &,--& und A, r,rP;n im allgemeinen kleine Größen erster Ordnung sind, kann man den 
letzten Summanden vernachlässigen und erhält das endgültige Resultat 


Ae& = HEN. [/4 TP;,Pıdr + [ArPıdr| an En Tara Re“ ) (3,4). 
[PinPr dr n bi 
0 


Diese Rechnung setzt voraus, daß man die homogene Lösung der Differentialgleichung mit ihrem 
Eigenwert kennt, die auch verhältnismäßig einfach zu bestimmen-sind. Man kann jedoch ohne 
diese Voraussetzung mit der gleichen Methode, die zur Abschätzung des Eigenwertes zur Erfüllung 
der Regularitätsbedingung (Abschn. 3) diente, eine einfache Beziehung erhalten, in der nur die 
Teillösungen der homogenen Differentialgleichung des gleichen Eigenwertes wie der inhomogenen 
benötigt werden. Die Rechnung ergibt, wenn man AF=0 setzt 
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T, 

1 

A Pu (außen) — A Pr (innen) P,(ro) [de(, > C,] f Fu aan a s 
far 
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0 0 


Die Funktionswerte der Außenintegration verschafft man sich aus dem Näherungspolynom und 
erhält somit eine eindeutige Beziehung zwischen A7 und As, um durch lineare Interpolation 
den richtigen Wert Ae zu gewinnen. 

In dem vorliegenden speziellen Fall der Differentialgleichung der Wellenmechanik wurde 
nach der Änderung des Eigenwertes gefragt, wenn man im Potential 7 (r) eindurch L. Biermann 
gegebenes additives Polarisationsglied der Form 


AI er f — ex» 5 ))) er Me (5,6) 


berücksichtigt, wobei a der sog. Polarisationsfaktor ist. Bevor hierfür gleich ein entsprechendes 
Beispiel zitiert werden soll, sei noch vermerkt, daß für homogene Differentialgleichungen (F=0) 
und AF=0 (5,4) für normierte Lösungen die einfache bekannte Form hat 


de= fa EPhadrs SR Sr Ua (5,7). 
ö 


Weiter unten wird außerdem gezeigt, daß eine analog gebaute Abschätzungsformel für inhomogene 
Differentialgleichungen 


de= J4 DT Pad N Dre (9,8) 


zu sehr ungenauen Ergebnissen führen kann, was die ganze Abschätzung bei der oft großen 
Empfindlichkeit des Eigenwertes illusorisch ınacht. 


Das hier gewählte Beispiel ist der Rechnung, die dem späteren numerischen Teil der Arbeit 
zugrunde liegt, entnommen. (Die Rechnung ist an der dort angegebenen Differentialgleichung II 


für Mgt+ durchgeführt worden.) 
1. Nach (5,8) ergibt sich mit AF = 0 
Ae= x [FPr,:.dr 
1) 


de = -- 0,001346 «. 
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a | | de = 0.000075 «.. DR Re 
u 3 Aus 6, 5) folgt durch Interpolation Bee A 
de = —0,000576 . AN 
Damit bestätigt sich, daß 6 8) nur für ganz grobe Schätzungen Draht ist. Die genauen 
Ergebnisse liefert die letzte Formel, die etwas mehr Rechenaufwand erfordert. Allerdings lohnt ı 
sich dieser Rechenaufwand vor allem dann, wenn man schon aus der ersten Schätzung der ersten 
Näherung auf den richtigen Eigenwert schließen will, ohne die Wellenfunktion im ganzen Bereich 
genau zu kennen. Es ist noch zu bemerken, daß diese Rechnung natürlich auch benutzt werden AR 
kann, um die umgekehrte Aufgabe zu lösen. In dem vorliegenden physikalischen Problem könnte 
man daran denken, die Größe & so zu ‚bestimmen, daß man die beobachtete Energie, dieim 
wesentlichen vom Eigenwert che erhält, um damit auch den Einfluß der PO AnyN richtig Ir 
zu erfassen. me 
Als eine weitere Anwendung I Störungsrechnung .. man die BILLEE, der Norn mnierungs- a R 
bedingung betrachten. 
: Man schreibt die Differentialgleichung in der I br 
FHRBdg tk Fer.) 822% a 
mit k als einen so zu bestimmenden Faktor, daß die Normierungsbedingung f. 


1/2 
RB k=a=|[riur| 


erfüllt ist. : 
Man hat zunächst eine Lösung für k=1 bestimmt, wobei sich das Normierungsintegral 


i | = 
: KA) — If Par) 
iv) 


im allgemeinen nicht zu Eins ergibt. Anderenfalls ist die integrierte Wellenfunktion normieıt. 
Setzt man in (5,9) wahlweise einen zweiten Wert k=Ä, ein, errechnet sich mit (5,4) die ent- 
F sprechende Eigenwertänderung, indem man AT=0 und AF=(k,—k,)F setzt, so ergibt die 
Integration von (5,7) das Normierungsintegral k®. Durch lineare Interpolation läßt sich nun- 
mehr der richtige Wert k=4 bestimmen und mit (5, 4) der zugehörige Eigenwert. Die gesuchte 
normierte Wellenfunktion ist dann 
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6. Näherungsweise Bestimmung des Eigenwertes & E: 
Um bei dem in der Einleitung beschriebenen Iterälionsverfahren möglichst schnel] zum 


- Ziel zu kommen, benötigt man eine möglichst genaue Eigenwertschätzung. | 
} Integriert man die inhomogene Differentialgleichung von Null bis oo, so erhält ınan mit 


Au  P(&)=0 nach e aufgelöst x 
_. po) ae) p + rar 


y2 


ee [ a (O 
[Pdr 
f h) 
Ist 2> 0 so ist auch P’(0)=0 und aus (6,1) wird 
ea )e+ de 
er: 120. ee 


f Pdr 


wobei es durchaus genügt, die Integrale nach Simpson zu bestimmen. 
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un 
“ * ” ” r * ben ist. Die Inte- 
Für /=0 beginnt P(r) mit einem in r linearen Glied, womit P'(0) gege ie. 
gration ist-aber deshalb ungenau, weil 7 im allgemeinen mit const/r beginnt und Fee pn 
kleine r erhebliche Ungenauigkeiten im Produkt T- P die Berechnung von & gerade im wesen 
lichen Teil des Zählerintegrals verfälschen. Um dieser Schwierigkeit zu begegnen, erg % 
man die Differentialgleichung mit r und integriert dann von Null bis ©. Es ergibt s 


für 1=0 # 
[le 9)? +0]ar & 


(rPd 
jr r 


€e> 


7. Asymptotische Lösung der inhomogenen Differentialgleichung unter 
Benutzung einer approximativen Darstellung für 


Es hat sich gezeigt, daß im vorliegenden Fall-der Differentialgleichung der Wellenmechanik 
Fr) in der Form dargestellt werden kann 


F= (0, ea "TE; DR 
mit | N 
VE ae PS ee (1). 
Die asymptotische Lösung der homogenen Differentialgleichung schrieb sich, in der Form 


’ ya 1 
Pte er rle) 


IPmiräe 


n=1 # 


| 
| 


mit 


Führt man diesen Ausdruck in die inhomogene Lösung (2,12) ein, so lautet mit 


r + Y 2 : ;; A ‘ 
fen fee), 


das Resultat 
A| o) 


Pa=0%GPr—- (UP aa 1 Bag ER INT UR TE ee. (7,2), 
a. 
j N 
w 
wobei 
Pin(ro) 
(= ——— 
u Pı(ru) 
bedeutet. 


Die Bestimmung der Koeffizienten des Polynoms P, Hißt sich auch einfach aus den schon 
bekannten Koeffizienten a„ durchführen. 
Setzt man 


so gilt für die b„ die Rekursionsformel 
ne—n + Non +) su LEN een (7,4), 
wobei b,= —1 zu setzen ist. 


Wie sich später in der numerischen Rechnung zeigen wird, sind die Werte ungenaner, 
abgesehen davon, daß die Rechnung auch nicht schneller durchgeführt werden kann. 


Schlußbemerkung 


Die in dieser Arbeit vorliegenden Ergebnisse basieren im wesentlichen auf zwei analytischen 
Erkenntnissen. 

Die eine ist die bekannte Möglichkeit der Darstellung von Lösungen inhomogener linearer 
Differentialgleichungen zweiter Ordnung durch die Lösung der homogenen Gleichung. 


B hung der Welle 


% Fr 
Se Fr = Die ee a Frkenntnis ist die Denn der Ben ana der vor- 

v gelegten Differentialgleichung der Wellenmechanik durch eine semikonvergente hypergeometrische: 2 
Reihe, wobei man einen Teil der Reihe zur Darstellung der Lösung abspaltt.. e IR 
A Die Möglichkeit ‘der Darstellung durch eine konfluente hypergeometrische Reihe nel f 

R. Hartreei in einer Arbeit ‚‚The wave mechanics of an atom‘‘ (Proc. Cambr. Phil. Soc. Vol. 24, 
(1928) erwähnt. Die hier zugrunde liegende asymptotische DE ABlaLSERUNB ist unter . 


5 der Annahme für das Potential V= T2= Chr (= Kernladungszahl) 
er RR RER EN a # DEREN 
RAN; dr +{@e Es es Kalle 1.388 
‚Auf die Praxis der Möglichkeit einer Durchtähring wird dann nich, weiter 


eingegangen, abgesehen davon, daß die angegebene Form zu ungenau: Is B- 


000. Damit ist der analytische Teil abgeschlossen und es folgt der sulföribehe) Teil, in dem Ve 
' einmal eine kurze Anleitung zur numerischen Rechnung gegeben wird und weiterhin die ‘Methode Br, 
‘an der auf andere Art durchgeführten Rechnung für das Zwei ggnenipekirum von Mt | 
einer Ber rien Prüfung unterzogen wird. ai j 
Dee Be 


III. Numerischer Teil © | | Br E. = 
5. Rechenschemata zur Durchführung der numerischen Rechnung io: 3 
Es ist zu berechnen - EEE, BZ 


Paar Helden =)" 5 E | z 


u=1 


N 


Dafür schreibt man 


P(a) = ee 4Iay a) RB 
Die Te oation der Exponentialfunktion läßt sich durch die Exponentialreihe sehr schnell, .. 2 Se 
durchführen. Es ist en ER 6% 


| wis el ras+l2e), en es 


wobei x, der Numerus des Tafelwertes ist. Im allgemeinen kommt man’bei'einer im Numerus  . 

dreistelligen Tafel mit der hier geforderten Rechengenauigkeit mit dem in, A.x linearen Glied 

aus, so daß die Rechnung für die e-Funktion in einem Rechengang. in ‚der ‚Rechenmaschine 

erledigt werden kann. Br 
Das Polynom 79 


ar a a a en 


an a = 


N 1 
Aa)=1 +2 an u 


berechnet man nach dem Hornerschen Schema. 

Für die numerische Integration zwischen festen Grenzen genügt im allgemeinen die Simp- 
‚sonsche Regel. Für die übrigen Integrationen ist es ratsam, unter Vermeidung der Aufstellung 
von Differenzenschemata folgende Formeln zu benutzen. 


| f fa)da= * Wfla— hl) +19fla) —5fta +h) Ha +24] 
h 


24 
a+h DE I | j 
fo)da= —[—fia +h) +13f(a) +13 fa +h) fa +21] f : : : (52. 
Erd N 
ade= zZLft@—h)—5fla) +19 fa +h) +9 f(a +2h)] 
a+h 


h5 j 
Der absolute Fehler ist im Mittel 3.4 T u) a—h<sEsa-+t2h. 


6. Er Schoungen 
Wie schon erwähnt, sollen die numerischen Ausgangswerte einer im wesentlichen für ‚das 
Zweielektronenspektrum des Mg+ (Singlett) durchgeführten Rechnung entnommen werden.. 


220 Schade, Numerische Integration einer Differentialgleichung der Wellenmechanik 34.29 Nr. 7/8 Juli/Aug. 1949 
PL ER RERE 2 SCHBErE SENDER EIER ENDE EB a 


Das System der Differentialgleichungen lautet: 

I. P'’(3s) +(T, — 8.) P(35) +F, = 0, 

11. P’(3p) +7, —&3,) P(3P) a Be 
wobei in T,, F,, T, und F, Integralausdrücke der Funktionen (35) und P(3p) selbst vorkommen. 
Das genaue Gleichungssystem läßt sich aus (1,6) ableiten, wenn man dort 4s durch 3s, 3s durch 
2s, 2s durch 1s, 4p durch 3p und 3p durch 2p ersetzt. P(1s), P(2s) und P(2p) sind in r nume- 
risch vorgegebene Funktionen, wobei über die Beschaffung dieser Funktionen hier nicht diskutiert 

j d 11 g ”* £ * ” 
en ER RE einer ersten Näherungsfunktion für die inhomogene Differentialgleichung I 
ird ächst die homogene Gleichung I gelöst. hi 

ie Ans ok Heckndne wird zunächst der erste Iterationsschritt zur Bestimmung des de 
herausgenommen. Anschließend wird für den letzten Iterationsschritt mit dem dort angegebenen 
Eigenwert das Näherungspolynom für r Zr, berechnet, während die Werte der Innenintegration 
nach dem bisherigen Verfahren übernommen werden. 


a) Abschätzung der Eigenwertänderung zur Erfüllung der Regularitätsbedingung für die 
homogene Differentialgleichung I 
Datenangaben. 
L’=0 
&3, = 0,463 erster geschätzter Wert 


„= = 0,03 Näherungsfunktion für die Außenintegration von homogen I 
r 


ei = e, — b= 0,493 Eigenwert zur Polynomrechnung für A e=0 
e&2 = 0,393 Eigenwert zur Polynomrechnung für de= — 0,1 


In = 


. 3,6 <r, <4, r, Wendepunktabszisse. 
r —— 3,6 D 


&: = 0,463 (de = 0) el) = 0,493 (A & = 0) | e{? = 0,393 (A & = —0,1) 


Pr) 


innen Par) außen P(r)außen 


0,644 332 0,644 332 
—0,072 557 — 0,036 382 
Damit gilt für die Interpolation nach (2,13) 
— 0,014105. + Ae - 0,561966 = — 0,072557 — 4e - 0,36175; 


das ergibt 
de = — 0,063, 
e= 0,400. 


Der in der dortigen Rechnung angegebene Interpolationswert für & unter Zulassung eines Fehlers 
des logarithmischen Differenzenquotienten für r=r, von 5% ist 0,396. Zur Ermittlung dieses 
Wertes wurden 2 Interpolationsschritte mit insgesamt 3 Außen- und 3 Innenintegrationen durch- 
geführt, während die beschriebene Methode lediglich eine Innenintegration und zwei Polynom- 
berechnungen der Außenintegration für die Anfangswerte 3,6 und 4 mit den Anschlußwerten 
bei r=4 notwendig machen. 

Selbst wenn e=0,400 eine etwas größere Abweichung der Regularitätsbedingung als 4A 
ergeben würde, was numerisch nicht nachgeprüft ist, würde die Genauigkeit dieses Wertes weitaus 
reichen, um damit als 1. Näherung in die inhomogene Differentialgleichung einzugehen. 


b) Berechnung des Polynoms für die Außenintegration der homogenen Differentialgleichung 
Datenangaben. 
a) 
&3 5 = 0,396 


„> 
T, = =— b Näherungsfunktion für r > nr. 


kei Die numerischen Rechnungen nach dem bisherigen Verfahren des Probierens und die hier angegebenen 
Versuchsrechnungen sind in der ‚‚Astrophysikalischen Abteilung‘‘ des Max-Planck-Instituts für Physik, Göt- 
tingen, durchgeführt worden, wo gegebenenfalls die Unterlagen einzusehen sind. 
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Bi A a ae Et 0 de 


imerischen ae aid: in Zahlen alel, 1 usammengestelt. 


Spalte 1: Exakte Werte der T-Funklion für 4 eiz 10. 

Spalte 2: A ae, füra= 2; b= —.0,02 mit Annäherung de T-Funktion für 

a zz 10, 

„Spalte a Polynomrechnung für a= — 2; ag — 0, 03 mit Annäherung der T-Funktion 
Bemzll;. 


Spalte 4: Numerische. DENE OL ac ‚außen mit den Ausgangswerten 0,408878 (I (r== =); Be: ; 


0,345834 (r = 5,5). 
Spalte 5: Be für Va 51730; b=- -0, 023 mit Annäherung der T-Funktion 
ürAsr s2. | 
Spalte 6: Numerische Integration nach innen mit den en k+0, 000221 In = 20), & 


k- 0,000293 (r = 19,5), wobei %k so bestimmt wird, ‚daß der Ordinatenanschluß ER 


beir=4 gewährleistet ist. 
Spalte 7: Bang Integration nach innen mit der bisherigen Methode des Probierens. | 


DEN OnCARSTT eg 

SO ( ERST ER) In ohre 

Polynomrechnung Dis DE P 

Pss P3s [a = 1,951730 N ne 

t 
0= 0,02) | 0 | 12 =— 0,08 EN 

0,480334 | 0,536830 | 0,5680 | — | 0,836830 | 0,536830 | 0,536 830 
0,416067 | 0,477526 | 0,474009 (ei 0,473185 | 0.474174 | 0.474173 


0,368480 | 0,414591 | 0,408678 | 0,408678 | 0,407576 | 0,408826 | 0,408 825 
0,331737 | 0,353138 | 0,345834 | 0,345834 | 0,344727 | 0,346087 0,346 085 
0,302345 | 0,296189 | 0,288253 | 0,288415 | 0,287362 | 0,288779 | 0,288777 
0,278147 | 0,245279 .| 0,237296 | 0,237670 | 0,236539 | 0,238153 | 0,238 152 
0,257775 | 0,200958 | 0,193309 | 0,193883 | 0,192760 | 0,194499 | 0,194497 - 
0.240319 | 0,163155 | 0,156075 | 0,156 771 0,155716 | 0,157541 0,157 539 
0.225156 | 0.131427. | 0.125046 | 0.125751 0,124847 0,126704 | 0,126 702 
0,211834 | 0,105148 | 0,099515 | 0,100098 | 0,099441 | 0,101276 | 0,101274 
0.200026 | 0,083620 | 0,078731 | 0,079057 0,078747 .| 0,080513 | 0,080511 - 
0,189477 | 0,066146 | 0,061962 -| 0,061893 | 0,062043 | 0,063700 | 0,063 698 
0,179992 | 0,052077 | 0,048538 | 0,047932 | 0,048658 | 0,050184 | 0,050 181 
Ei ee ee Br 0,038004 | 0,039383 | 0,039 380 
2 0,029574 0,030798 | 0,030796 
u 0,022935 0,024010 0,024007 
2 0,017733 | 0,018664 | 0,018661 
Ir we 0,013672 | 0,014470 | 0,014.467 
224 0,010514 | 0,011 191 0,011 188 
$ 0,008066 | 0,008636 | 0.008 633 


£ 0,006 175 | 0.006651 | 0,006 648 

#7 r= 0,004716 | 0,005 112 | 0,005 109 
Ri Bi 0,003 597 0,003 923 0,003 919 
re = f 0,002738 | 0.003004 | 0,003 001 
En ER 2 0,002080 | 0,002 297 0,002 294 
A ar I 0,001 577 0,001 754 | 0,001 751 
ar Ei re 0,001195 | 0,001 337 0,001 333 
13 Re Ws 0,000 904 | 0,001.017 0,001 013 
= 24 0,000684 | 0,000773 | 0,000 768 

si ee; 0,000516 | 0,000588 | 0,000 579 

a 2. 0,000389 | 0,000446 | 0,000.435 

I 0,000 293 0,000338 | 0,000. 323 
ne: 0,000 221 0.000255 | 0,000 236 


0,000249 | 0,000213 


Die Ergebnisse der Spalle 2 und 3 sollen zeigen, daß die maximale Dilferenz zwischen a/r 
und dem exakten Wert von 7 der geeignetste Wert von 5 ist, was sich auch in allen anderen 
Rechnungen zeigte. In Spalte 5 wird gezeigt, daß es im Durchschnitt günstig ist, die Annäherung, 
der T-Punklion über den größten betrachteten Bereich zu erstrecken. 

Eine numerische Integralion nach außen (Spalte 4) erweist sich als wesentlich schlechter 
als eine Integration nach innen (Spalte 6), wenn die Ausgangswerte der Polynomrechnung ent- 
nonmen werden. Es zeigt sich also, daß man die Polynomrechnung auch lediglich dazu benutzen 
kann, gute Ausgangswerte für eine numerische Integration zu gewinnen, ohne die Methode des 
Probierens anwenden zu müssen, wenn der Eigenwert zuvor zur Erfüllung der Regularität im 
Wendepunkt nach obigem Verfahren bestimmt wurde. 

In jedem Fall ist die Genauigkeit jedenfalls den Anforderungen einer Iterationsrechnung 


vollauf gewachsen. 


....°e) Beispielrechnungen für die inhomogenen Differentialg 
Für beide Differentialgleichungen I und II wird aus den best F 
_ der letzten. Iterationsschritte a Es yo ar " ah ze 
Doppelintegration sowie der einfachen Integration gezeigt, wie weit die Polynomrec 
ae nächst Gleichung I betrachtet. 


Genauigkeitsanforderungen genügen. Es wird zu 


Daten ae . 


Yy =4. 
Egg = 0,680 # ch . 
7, = = 9,11 Näherungsfunktion für # 210 
€ = 0,790 De 3 | | RRSERE 
FF, = 0,009308 - e Var J—K Näherungsfunktion fürs 2m. 00. 


Zahlentafel II gibt die Lösung der homogenen Differentialgleichung für r Zr, als Polynom- x 


Yr 


“ rechnüng an, sowie die Werte der inhomogenen Näherungsintegrationen aus der ‚homogenen a 
Polynomrechnung und zum Vergleich die Werte der numerisch nach Numerow integrierten 


exakten Differentialgleichung. 
Zahlentafel II 


@ 


3 
P}n(@s) = ch 
Doppelintegral ji 
4,0 0,571.039 0,413 025 — ((0,404 523)) 0,413 025 0,413 025 
4,4 0,475 308 0,331 661 En . ((0,325 468)) 0,330 747 0,332 563 
4,8 0,389 161 0,261 216 (0,261 216) ((0,256 423)) 0,259 484 0,262 237 
5,2 0,314 329 0,202 244 (0,202 244) (0,198 460)) 0,199854 0,203 157 
5,6 0,251 012 0,154 151 (0,154 108) (0,151 195)) 0,151 301 0,154946 
6,0 | ° 0,198 514 0,115 751 (0,115 679) ((0,113 512)) 0,112 649 0,116 479 
6,4 0,155 694 0,085 636 . (0,085 535) ((0,084 009)) ' 0,082 485 0,086 338 
6,8 0,121 226 .. 0,062 387 (0,062 229) ((0,061 263)) 0,059 354 0,063 077 
7,2 0,093 787 0,044 695 (0,044 424) ((0,043 961)) 0,041 893 0,045 366 
7,6 0,072 149. 0,031 411 (0,030 953) ((0,030 965)) 0,028 933 0,032 045 
8,0 0,055 226 0,021 574 (0,020 835) ((0,021 324)) 0,019 463 - 0,022 146 
8,4 0,042 083 0,014 387 (0,013 264) ((0,014 262)) 0,012649 _ 0,014881 
8,8 0,031 933 0,009 213 (0,007 586) |. ((0,009 161)) - 0,007 858 0,009 619 
9,2 0,024 145 0,005 553 (0,003 273) ((0,005 537)) 0,004 535 0,005 866 
9,6 0,018 195 0,003 013 (—0,000 103) ((0,003 013)) 0,002 298 - 0,003237 
10,0 0,013 668 0,001 300 (—0,002 899) ((0,001 300)) 0,000 847 0,001 437 
10,4 —_ 0,000 182 E= — En 0,000 238 
10,8 — —0,000 507 — — — —0,000 528 
11,2 = —0,000 891 .— . _ —0,000 989 
11,6 — —0,001 060 _— —_ — —0,001 236 
12,0 = —0,001 077 = _ = —0,001 340 
12,4 — —0,000 984 — _ _ — 0,001 348 
12,8 — —0,000 807 -—_ — _ —0,001 296 


Die homogene Rechnung ist zur möglichst genauen Bestimmung der Konstanten C, bis 
r=20 durchgeführt worden, während die innomogenen Näherungsrechnungen im wesentlichen. 
zunächst analog der Bestimmung der Näherungsfunktionen für 7, und F, bis r= 10 erfolgt Sind. 
(Die Zahlentafeln hierzu konnten aus Platzmangel nicht veröffentlicht werden.) 

Die unter dem waagerechten Strich stehenden Werte sind lediglich zum zusätzlichen 
Vergleich gerechnet worden. 

Es zeigt sich zunächst, daß die näherungsweise berechneten Funktionswerte P;»(3s) in 
dem Bereich von r=4 bis 10 mindestens für eine Iterationsrechnung vollauf genügen, wobei die 
Werte der Doppelintegration genauer sind. 

Im weiteren Verlauf über r= 10 hinaus wird die Rechnung ungenauer, da der Integrand 1/P% 
sehr groß wird und damit die Auswertung des Integrals ungenau, wenn man die Intervallängen 
des Arguments beibehält. 

Außerdem lohnt sich eine Weiterrechnung schon deshalb nicht, weil die Funktionswerte 
so klein sind, daß sie ohne Einfluß auf die gesamte Iterationsrechnung bleiben. 

Es ist nur interessant festzustellen, daß für große Werte von r noch kleine Extremwerte 


der Funktion vorhanden sind, bevor die innomogene Lösung in die positive Lösung der homo- 
genen Gleichung übergeht. 
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übereinstimmung und Regularität besonders mühsam ist, BR u Be. 
Die in Spalte 2 stehenden eingeklammerten Werte [( )]erhält man wiederum,durch nume- 


rische Integration nach außen mit den Anfangswerten P;, (4,8) = 0,261 216 und P;,(5,2)=0,202224, 


die Ergebnisse. sind schlechter als die der Doppelintegration. Die weiteren eingeklammerten _ 
Werte in Spalte 2 [(( .))] ergeben sich. aus der numerischen Außenintegration mit den Anfangs- 
werten der Polynomrechnung bei r=10 und 9,6. Der Differenzenquotient der logarithmischn 


Ableitungen der Innen- und Außenintegration und der Ordinaten am Wendepunkt differieren 
‚etwa um 2%. Man kann hier schon bemerken, daß es keine große Schwierigkeit bedeuten würde, 


die nun schon für r=10 und 9,6 ziemlich genau bestimmten Ordinaten so zu ändern, däß vor 
allem die Ordinatenübereinstimmung bis auf 10%/,, genau wird, Dasgleiche ergibtsich auch später 
in der Rechnung für Gleichung II, wobei die Übereinstimmung noch besser wird, dauch de 
T-Funktion über einen größeren Bereich angenähert wird. Die Abschätzung der Eigenwert- 3 
änderung zur Erfüllung der Regularität soll hier für den 1. Iterationsschritt der Gleichung I BR 
gerechnet werden. i \ HR AR hr 

Die 1. Eigenwertschätzung wurde mit e=0,77 angesetzt, e wurde wieder mit — 0,1 variiert, 
wobei die Ordinatenübereinstimmung für r=3,6 durchgeführt wurde. | 


% 


e=0,77(4e= 0) 


4,0 i 
P ("innen Cop @ 1/Pi dr 


- 0,421878 0,421878 Br: 


Bo 0,421 878 2 
| -0,471391 | er 
4,0 0,369 299 | | 0,338 318 0,355 478 ı 
AP —0,052579 | —0,083 560 © —0,066.400 % 
de = —0,048 RR INN 


Damit berechnet sich der Eigenwert unter Berücksichtigung der Regularität zu e=0,722. 


Ein mehrfaches Probieren ergab bei einer zulässigen Differenz der logarithmischen Ab- 
leitung am Wendepunkt von 1% einen Eigenwert von e=0,725. Zur Ermittlung wurden mehrere 
Innen- und Außenintegrationen benötigt, wobei der Eigenwert und Anfangswerte für die Außen- 
integration nach bestem Ermessen variiert werden mußten. Mit der oben erwähnten Methode 
kann man durch eine grobe Variation von e unter Benutzung der homogenen Teillösungen für 
innen und außen in einem Interpolationsschritt den richtigen Eigenwert ermitteln, wobei man 
für die inhomogene Lösung nur einige. Punkte im Bereich des Wendepunktes benötigt. Ein 
Probieren bei großen r-Werten ist damit nicht mehr nötig. Für den so ermittelten Eigenwert 
kann dann die richtige Lösung mit regulärem Verhalten im Wendepunktbereich durchgerechnet 
werden. Eine analoge Rechnung ist auch für Gleichung Il durchgeführt worden. Die Ergebnisse 
sind in Zahlentafel III aufgeführt. | 


Die einzelnen Daten sind 


il = 

ER 9 

&p= 0,2373 

= Z — 0,005436 Näherungsfunktion für r 27, 

& = 0,242736 

F, = —0,000446 - e 7%" „5°”® Näherungsfunktion für r 2r.. 


Es zeigt sich auch hier, daß in jedem Fall die Auswertung des Doppelintegrals vorzuziehen 
ist, da dieses mit einem nicht größeren Rechenaufwänd die genauesten Werte liefert. 

Die in Spalte 2 (Zahlentafel III) in Klammern stehenden Werte ergeben sich aus.der nume- 
rischen Außenintegration der Differentialgleichung II mit den Anfangswerten des Doppelintegrals 
für r=20 und r=19,6. Die Differenz zwischen dem logarithmischen Differenzenquotienten der 
Innen- und Außenintegration ergibt sich mit 1%, die Ordinatendifferenz am Wendepunkt 
mit 11/,%. Mit der schon verhältnismäßig guten Ordinatenübereinstimmung läßt sich unschwer 
eine Änderung der Ordinaten bei r=20 und 19,6 ermitteln, um auch bei r=r, eine Ordinaten- 


u 2 " 
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(l) | ( 


T Polynom lintegral 
Pı(3») B, (3) 

7,2 1,730 208 —1(),355 720 (—0,350 732) 

7,6 1,596 149 —0,331275 (—0,326 953) 

8,0 1,465 183 —0,306 863 (—0,302 736) 

8,4 1,336 627 —10,282 377 (—0,278 622) 

8,8: 1: 1213211 0,258 432 (—0,255 035) 

9,2 1,096 153 —0,235 372 (— 0,232 298) 

9,6 0,986 261 — 0,213 400 (—0,210 641) 
10,0 0,883 996 —0,192 684 (—0,190 223) 
10,4 0,789 551 . —0,173 323 (—0,171 143) 

10,8 0,702 910 —0,155 365 (—0,153 448) 

11,2 0,623 900 —0,138 821 (—0,137 149) 
11,6 0,552 230 —0,123 671 (—0,122 222) 
12,0 0,487 526 —0,109 869 (—0,108 624) 
12,4 0,429 363 —0,097357 | (0,096 295) 
12,8 0,377 284 —0,086. 062 | (—0,085 164) 
13,2 | .0,330 819 —0,075907 | (0,075 155) 
13,6 0,289 500 —0,066809 | (0,066 184) 
14,0 0,252 867 —0,058684 | (—0,058 170) 
14,4 0,220 479 —0,051 449 (—0,051 032) 
14,8 0,191 919 —0,045 027 (—0,044 692) 
15,2 . 0,166 795 —0,039 340 (—0,039 074) —0,041 539 —1,038 522 
15,6 0,144 743 —0,034315 | (--0,034 108) —0,036 386 —0,033 575 — 1,034 332 
16,0 0,125 429 —0,029 885 (—0,029 727) —0,031 828 —0,029 220 —0,029 868 
16,4 0,108 546 —0,025 988: (—0,025 869) —0,027 803 —1,025 393 —0,025 930 
16,8 0,093 815 —0,022566 | (--0,022 479) —0,024 256 —0,022 038 — 0,022 461 
17,2 0,080 985 —0,019564 | (0,019 504) —0,021 136 —D,019 102 —0,019 405 
17,6 0,069 828 —0,016 938 (—0,016 897) —0,018 396 —0,016 537 — 0,016 717 
18,0 0,060 141 —0,014 641 (—0,014 616) —0,015 993 —0,014 299 —0,014 349 
18,4 0,051 743 —0,012 635 (—0,012 621) —0,013 890 —10,012 350 | —0,012 263 
18,8 0,044 472 —0,010 885 (—0,010878) —0,012 051 -—0,010 656 —0,010 422 
19,2 0,037 805 —0,009 266 (—0,009 357) —0,010 341 —0,009 093 —0,008 790 
19,6 0,032 785 —0,008029 | (0,008 029) —0,009 045 —0,007 908 —0,007 339 
20,0 | 0,028 077 —0,006869 | (—0,006 869) —0,007826 - | 0,006 803 —0,006 400 


übereinstimmung mit einem Fehler von 1°/,, zu erreichen. Man erkennt hieraus, daß-die Aus- 
wertung des Doppelinlegrals auch für. eine numerische Außenintegration sehr gule Ausgangs- 
werle liefert. 


Zusammenfassung 


Die numerische Durchführung der Versuchsrechnungen hat gezeigt, daß sich im Falle des 
Ein- und Zweielektronenspektrums eine analylische Behandlung des Differentialgleichungs- 
systems vor allem für das asymptotische Verhalten lohnt. Die Methode des Probierens zur 
Erfüllung der Regularitätsbedingung kann durch analytische Überlegungen erheblich abgekürzt 
oder ganz fallen gelassen werden, womit sich im Rechnungsgang eine bedeutende Zeitersparnis 
ergibt. 

Es hat sich vor allem gezeigt, daß im Bereich der r-Werte, in denen die Außenintegration 
numerisch begonnen wird, im homogenen wie im inhomogenen Fall keineswegs eine reine Expo- 
nentialfunktion in » anzusetzen ist. Die damit gerechneten Anfangswerte sind weitaus zu klein. 
Im homogenen Fall läßt sich dieser Fehler hinterher durch einen Faktor korrigieren, was aber 


x 


keinesfalls im inhomogenen Fall möglich ist. 


Die aufgezeigte Integration der asymptotischen Differentialgleichung ergibt dafür recht 
gute Ausgangswerte für eine numerische Integration. 

Das Iterationsverfahren bleibt im Prinzip unberührt. Hierbei ist lediglich auf einen Weg 
verwiesen, die Anfangswerte der Integration vom Nullpunkt so zu bestimmen, daß die integrierte 
Funktion nahezu normiert ist, so daß eine nachträgliche Normierung keine wesentliche Änderung 
der rechten Seite der Differentialgleichung ergeben würde. 


Eingegangen: 5. 11. 1948. 


we Das Bildungsgesetz für die Fehlerform. | | 
von fehlerhaften Meßreihen mit Hilfe ganzer rationler 


[x FIR; Fr 


"Funktionen wachsender Ordnung 


= Von K. Stange in Karlsruhe Kt | ü > 2 “ 


Die beim Ausgleichen von Meßwerten mit Hilfe von endlichen Potenzreihen geltenden Gesetze für di 
Fortpflanzung der Meßfehler und für die Fälschung der Funktionswerte werden Ei nk 
aufgestellt. Indem man diese Determinanten nach fallenden Polenzen der Punktzahl n entwickelt, findet 

. man: äußersi einfache, für ‚‚große‘“‘ Punktzahlen n gültige ‚‚asymptotische‘“ Gesetze. BD 
The laws valid for In a of errors and the falsifisation of Ihe yalues of the function in the case. 
ai Ian Kara { wei lues by ae . ae ar stated in a „‚delerminant form“. By expanding 
e determinant in decreasing powers of Ihe number n of the measuring points, extremly simpl i 
laws are found valid for large numbers n. u ER 
3aKOHB, AMeHCTBHTeNBHLE IA PACHPOCTpaHeHuAd Morpemmocteä u AA HCKa:KeHuA 
SHayeHHÄ DyHEIIHH B CJIyyae BEPABHHBAHHN Pe3yIBTATOB IIPH IIOMOIIH KOHEYHEIX CTEMEHHEIX 
. PAAOB, IIpencTaBimmrtca B dopMe onpenennteuei. Ilpmn packpsrrun aTux onmpenennrenet 
‚NO yÖBIBAWINHM CTelleHAM YHCH& n TOYeK H8MePeHHA HONYYAITCH OYeHb TIIPOCTEIE . „ACHM- 
_ ATOTHyecKHe“ 3AKOHB, NEÄCTBUTENBHHE MIA O6ONBIIHX YHCEN n.| i 


| -1. Einleitung a 
Über die Verfahren des Ausgleichens fehlerbehafteter Meßreihen gibt es ein ausgedehntes 


Schrifttum). Trotzdem stehen zwei Fragen dabei noch verhältnismäßig im Hintergrund, nämlich 
die. Auswirkung fp des Meßfehlers u der Ausgangswerte x auf die geglätteten Werte & (und ihre 
Ableitungen %, ...) und die mit dem Glätten zwangsläufig verbundene Fälschung fu dieser. 
Werte. In welcher Weise der „physikalische“ Fehler fr und der ‚‚mathematische‘ fy des ge- 
glätteten Wertes z von dem benutzten mathematischen Aufwand, d.h. von der Ordnung R der 
Ausgleichfunktion und von der jeweils zusammengefaßten- Punktzahl N, abhängen, wurde im 
Falle ganzer rationaler Ausgleichfunktionen niedrigen Grades R kürzlich erörtert2). Da in jener 
Arbeit praktische Gesichtspunkte im Vordergrund standen, so war die dort vorgenommene 
Beschränkung auf den Fall R</5 kein Nachteil. Da jedoch das Bildungsgesetz, dem die Fehler- 
gleichungen mit wachsender Ordnung R der Ausgleichfunktion gehorchen, für fr nur andeutungs- 
weise und für fy gar nicht hervortritt, soll die Frage im folgenden wenigstens für die Funktions- 
werte und ihre erste Ableitung noch einmal in größerer Allgemeinheit aufgegriffen werden. Dabei 
wird es gelingen, für ‚‚große‘‘ Punktzahlen. n geltende ‚‚asymptotische‘ Gesetze von sehr einfacher 
Form aufzustellen. Diese Gesetze lassen sich dann mit Hilfe einer von Willers angegebenen 
Lösung der Ausgleichaufgabe®°), die für ‚„‚kleine‘‘ Punktzahlen n gültig ist, auf „‚das unbekannte 
Zwischengebiet ausdehnen‘. 


2. Das Bildungsgesetz für Le, und Ma,+9 
Glättet man eine ausgedehnte Reihe von Meßpunkten ..., %—1, &ı 41 +. welche 


in gleichen Abständen A#= aufeinanderfolgen, indem man abschnittsweise durch je N=2n-+1 
aufeinanderfolgende Punkte %—n, :.. Zp +++, Zk+n die „beste‘‘ ganze rationale Funktion 
vom Grade R, 


RR (&): 1 
Dis Tale. Re ee LE VA EUER, 
a2.) ...) 
legt, so ist der ausgeglichene Wert x; an der Stelle #; 
+n 
EN DE ee A Ne Er LA 
v=—n \ 


Die Gewichte A,, mit denen man die einzelnen x-Werte zu multiplizieren hat, sind in der Form?) 


1 Se ... Ser 
1 y2 5 a S2r+2 


kh= . . ß 
v Ds, a (3) 


var Sar+2 ae Sy 


1) Fr. A. Willers: Methoden der praktischen Analysis, S. 239 u.f. Berlin und Leipzig 1928. (Im 
folgenden Willers Igenannt.) —C. RungeundH. König: Numerisches Rechnen, S. 189 u. f. Berlin 
1924. — Blaschke: Vorlesungen über mathematische Statistik, S. 192 u.f. Berlin und Leipzig 1906. — 
Whittacker und Robinson, The calculus of observations, S. 285 u.f. London 1926. 


2) K. Stange: Die zweekmäßige Auswertung von punktweise aufgenommenen Zeit-Weglinien. Ing.- 


Arch. 1948, S. 383. 


5) Fr. A. Willers: Das Glätten empirisch gefundener Zahlenreihen. Z. angew. Math. Mech. 10 (1930), 


S. 166. (Im folgenden Willers II genannt.) 


4) Stange: a.a.0. S. 387, Gl.(10,) und (9). 
15 


eln beim Ausgleichen | Fr 5 


MM 
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Gr 
Ar: 
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= a x 
> 
Fa} 
1s2 
SEE 
4 ri 
a 
ar 
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5 > na en >. La j 
EN RMER u ea; SERIE ei . 
y n ee I e >, ry e > ev ARFE | N 41% 
ale je | a ae ee. 
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darstellbar. Mit Dar bezeichnen wir allgemein eine r-reihige 


Element die Summe 8, besitzt und deren übrige Elemente „von 8, aus“ der ch aus- 
gedehnten Matrix } par I 
S 8, S 9 
E Sa S Sg RE . . a a Te (4) 
8 De De A 
entnommen sind. Br 


Die Zahlen S, sind die Potenzsummen 
-+n n 
S=2n +1; = zm=22W; p gerade und >2 .....0). 
van y= F 


Wir setzen voraus, daß die ursprünglichen s-Werte mit dem Meßfehler u behaftet sind. 
Nimmt man für alle an der Bildung des gewogenen Mittels 7, beteiligten z-Werte &—n, . -- Ze+n 
die gleiche Meßgenauigkeit an, was bei dichter Meßfolge erlaubt ist, so ist der „physikalische 


Fehler‘ von &x 
fe(R; N)= 2» re er . (6). 


Der ‚mathematische Fehler‘ fy wird bei Verwendung einer ganzen rationalen Funktion vom 
Grade R=2r oder R=2r +1 in erster Näherung?) 


1 : L 
Ju(R; N) == Zr +2)l zer+2) err2 Ma3++2 Er (7,)» 
wobei ENKaE | 
x g Fri 
Mop=2 Mr2,=(—1f u ee (75) | 
5 y—n Do, +1 


das auf die Bereichmitte bezogene Moment der Ordnung 2r-+2 der Einflußzahlen 2, ist. 
Maßgebend für Fehlerfortpflanzung und Fälschung sind also letzten Endes die Ausdrücke 
L;, und Ma,+2. Durch (7,) ist Ma,ı+2 mit Hilfe der Determinanten D,,,+1 ausgedrückt. | 
Für L,, wollen wir die entsprechende Form noch herleiten. 
Multipliziert man (3) mit A,, summiert die entstehenden Ausdrücke über v von —n bis +n 
und berücksichtigt dabei, daß & },=1 ist und daß sämtliche Momente M, = 3 »* A, für 
v 


v 
2<08<<2r verschwinden, so vereinfacht sich die auf der rechten Seite erscheinende Determi- 
nante zu D;,,. Damit haben wir das Bildungsgesetz 


+n 
SER FE le (8) 


y=_—n Do r+1 
für die Summe der Gewichtsquadrate gefunden. 


Man sieht, daß es möglich ist, mit (8) und (7,) allgemeine Aussagen über Fehlerfortpflanzung 


und Fälschung zu machen, wenn es gelingt, die Determinanten D,, ‚+1 der Berechnung zugänglich 
zu machen. 


Aus (8) und (3) folgt im übrigen ohne jede Vernachlässigung für alle N und R die Beziehung 


+n 


Pe N (8), 


van 


d.h. die Summe der Gewichtsquadrate ist gleich dem Gewicht A, des Ausgangswertes x; in der 
Mitte des Bereichs. 


3. Umformung der Determinante Dr,r+1 


Die Berechnung von D,,r+1 ist zwar in voller Allgemeinheit (d.h. für beliebige Werte 
von nundr) nicht möglich. Wir werden jedoch für „genügend große‘ Punktzahlen n zu brauch- 


baren einfachen Ergebnissen gelangen. Für die Summen S, gilt nämlich die Entwicklung®) 
nach Potenzen von n 


N « 


’ 27 Hallo re]... 9). 


°) Stange: a.a.0. 8.392, G1.(23,) und $. 387, Gl. (11,). 
%) Willers I: 2.8.0. 8. 127, Gl. (13). — Whittacker: a.a.0. 8, 29. 


- die Bernoullischen Zahlen?) sind. Ersichtlich darf S,(n) mit wachsendem » immer besser durch 
den Ausdruck N er a Es. 


| 5 +] er - SE (%) 


ersetzt werden, der bereits eine brauchbare Näherung darstellt, wenn n>p-1 ist. Setzt man 
 8p(n) gemäß (9,) in die Determinante D,,,+ı ein, dann läßt sich aus jeder Spalte der Faktor 
2 n? herausnehmen. Außerdem liefern die Spalten 2, 4, .. .,‚2r noch die Faktoren na, nd, .. „na 

Die gleichen Faktoren lassen sich aber auch aus den Zeilen 2, 4,..., 2r absondern. Insgesamt 


kommt el TEE 
En ER De £ DH = RI E a Ay,ı+1 re Bene i re ? ... (10,) 
N N n ah: m ER r 
BE WS pa wre EN 
N id H N 1 u MN a! ; 
er enagete e gT eio 
33 . = \ 
Den RN a, An RR 
RE A RE T ER  FET ET 


In der Determinante A,,,+ı ist jedes Element in bezug auf n von erster Ordnung. Entwickelt 
man sie nach Potenzen von n,so beginnt die Entwicklung mit n’+!. Mit (10,) wird aus (8) und (7,) 


R nitnmtdr Ar 1 Ar 
Denen a on. EAN ROSA NE - (11) 
tina Ayrrı 2 Ao,rtı 
und | 
2 Ir+l 2(r+1) „2r (r+1) 
en 
2HINErUED  Aoyers Ao,r+1 - 
Da A,,r in bezug auf n vom Grade r ist, so streben die Quotienten 
n As, r SE Aeier 
—— =R,(n und —— =m RN ES 
FAR Zr (n) Sr 2r+2(N) 2 (13) 


mit wachsendem n gegen feste nur noch von r aber nicht mehr von » abhängige Zahlen, die 
wir mit 137 und m&++2 bezeichnen wollen, um anzudeuten, daß sie nur von der Ordnung R der 
Ausgleichfunktion aber nicht mehr von der jeweils zusammengefaßten Punktzahl N abhängen. 

Mit dem bisherigen Ergebnis sind wir imstande, schon ein paar wesentliche allgemeine 
Eigenschaften der Fehler fp und fu auszusprechen. Der physikalische Fehler (6) hat nämlich 
die Gestalt 


oder in erster Näherung 


und für den mathematischen Fehler (7,) gilt 
| | (Dr 


MN nr (meer timgrte(n) . (15) 
oder in erster Näherung 
| ermn 
fu(R; N) = Sera NEE AR (15,). 


Für „genügend große‘ Punktzahlen » ist demnach der physikalische Fehler fr des ge- 
glätteten Wertes 2 bei Verwendung von Ausgleichfunktionen beliebiger Ordnung R in erster 
Näherung umgekehrt verhältnisgleich zur Wurzel aus der Gesamtpunktzahl N=2n +1. 


?) K. Knopp: Theorie und Anwendung der unendlichen Reihen, 8.182. Berlin 1924. 
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Der mathematische Fehler fi ist bei Verwendung von Ausgleichfunktionen der Ordnung 


 R=2r und R=2r +1 in erster Näherung verhältnisgleich zur (2r +2)-ten Potenz der „‚Bereich- 


breite“ B=nr. fy hängt also (eigentlich nicht von der Punktzahl » sondern) nur von dem 
Produkt nr ab Ak wird infolgedessen nicht (wesentlich) beeinflußt, wenn man die Zahl N der 


in einem festen Bereich 2 B= konst enthaltenen Meßpunkte ändert. 


4. Berechnung von fp und fm in erster Näherung 
Um weiteren Aufschluß zu gewinnen, haben wir die Ausdrücke 1%, und m$,+g in Ab- 
hängigkeit von R bzw. r zu bestimmen. Gemäß (13) kommt das im wesentlichen auf die Berech- 
nung der Determinanten Ay,r+ı hinaus. Erselzen wir die Summen 8, (n) durch die Näherung (9), 
so gewinnen wir in der Entwicklung von A,,,+1 nach Potenzen von n 


Aral JRR I Earl 


die ersten zwei Glieder. 
In dieser Entwicklung ist der Faktor von n’+! die (von n unabhängige) Determinante 


1 1 1 
vv. p+1 p+3 PrF2r +1 
1 1 1 
Ar p+3 p-+5 »+2r +3|1....... (17) 
p+2r +1 p+2r +3 p+4r-+1|]: 
und der Faktor von 7 n’ wird eine Summe aus (r+1) Determinanten, 
Deus Anrana ct Anus Ara FRA er I EB (18), 


welche entstehen, indem man die Elemente der Spalten 0, 2,4, ..., 2r von A3,r+ı durch die 
Wertel ersetzt. (16) erscheint also bei Berücksichtigung der beiden höchsten Glieder in der Form 


1 
An, r+1 — n’+1 As: r+1 E= Dy N n r+1 =» 6 el) fe; dere Dr (19). 


Wir wollen zunächst zwei Rekursionsformeln zur Berechnung von A3,-+ı herleiten. Behält 
man die Spalte 2r bei und subtrahiert ihre Elemente von denen der Spalten 0, 2, ..., 2r—2, 
so läßt sich aus den Spalten 0, 2, 4,... ., 2r— 2im Zähler der Faktor 2r, 2(r—1),2(r—2),...,2-1 
und aus den Zeilen 0, 2,...,2rim Nenner der Faktor (p+2r +1), (p+2r+3),..„ (p+4r-H1) 


herausziehen, 
Läßt man dann die Zeile 2r ungeändert und subtrahiert ihre Elemente von denen der 


Zeilen 0, 2,..., 2r— 2, so erscheinen in der letzten Spalte mit Ausnahme des letzten unver- 
änderten Elements 1 lauter Nullen. Die Determinante vereinfacht sich demnach, wenn man 
noch ähnlich wie eben alle gemeinsamen Faktoren -herauszieht, zu 
A* (2’ ri)? A* 
PT (p+2r HIElp +2r +38 (ptar—lep tar Hr 
Mit der vorstehenden Rechnung haben wir in der Ausgangsdeterminante A}, ,+1ı sowohl die 
Spalte 27 als auch die Zeile 2r ‚„‚ausgelöscht‘. Man kann nun aus den Spalten 0, 2,.. ., 2rirgend- 
eine, z.B.2& und aus den Zeilen 0,2,..., 2r irgendeine andere, z. B. 28, auswählen und diese 
durch den gleichen Rechengang, wie er eben ausführlich erläutert wurde, „zum Verschwinden 
bringen‘. ‚Dann entstehen Rekursionsformeln anderer Gestalt, von denen wir die für «=0 
und ß=0 noch hinschreiben wollen. Bringt man in der Ausgangsdeterminante A} +41 die 
Spalte O0 und die Zeile 0 zum Verschwinden, so gewinnt man die Beziehung 


« (20,). 


die für p=0 im Hinblick auf (13) nützlich ist. 
Nach (13) und (19) gilt in erster Näherung für genügend große Punktzahlen, d.h. bei 
Beschränkung auf die höchsten Potenzen von n, 


und 


; Rekursionsformeln. na und 20). um, so 6 erhält man i 


di mit Hilfe 
‚schließlich die von n Be Grenzwerte. 


ee Fr 


ee: ie | = 2r +1) 
RE a0 | IREVe: er +3) TE - (dr +3) | Bi: 
eg Setzt. man diese Werte schließlich in (14,) und (15,) ein, so a ınan die Fehler N "m Re 
‘ und Ju(R; N) in erster Näherung berechnet. 


£ E ne DR Für die Gewichte N Aalet aus (Si in Verbindung mit (10,) und as) in era Näherung SEE: vr 
“ ER I Y; "S Se 'Ir 1 
Ben be) 2 ! 
anhy= n 5 2r +3 
‚r+1 . 


E 2 ae ET N, 
Be}. | ER, Ru ar, dry om 
Der Nenner Ad,r+ı läßt sich in der Form | 


en [ r—o)te! he 1 
. E Hg DarEi.lTit (3 +20) ° A +1 Ft 
ie ea ER et Be ARE 


’ ' darstellen, wobei Az, die im Zähler. und Nenner gemeinsamen Faber der ersten 
2 Spalte sind. Die erste dieser Gleichungen für e=0 stimmt mit (20,) überein, wenn man dort 
p=0 setzt. Entwickelt man nun den Zähler nach Potenzen von (v/n)2®, so findet man schließlich 

\ mit Rücksicht auf die letzten Gleichungen als asymptotisches Gesetz der Gewichtsverteilung 


ER 
ur 7 Are 1 @+20)- BEEEZTIHITEE el) 


x 5. Berechnung von fein und fu in zweiter Näherung 
B: Daß die Gleichungen des vorausgehenden Abschnitts für sehr große Punktzullen‘ n die | 
5 wirklichen Verhältnisse gut wiedergeben, steht außer Zweifel. Um sie jedoch auch bei praktisch _ er 
> vorkommenden ‚‚kleineren‘“ Werten von n nutzbar zu machen, wollen wir im folgenden die hi 
zweite Näherung noch herleiten. Dazu müssen wir in der Entwicklung (19) den Faktor Sisrtt s 


1 
von — n" berechnen. Die Determinante Ay,,+ı läßt sich umgestalten zu 


B 2 

7 = (— ne a ERST, 

; Ad,r41 = ( in (+1 +20) - -(p-+1 +20 +2r) Ap,r+1,20> 

2 a 

4 Man findet diese Gleichungen, wenn man die vor (20,) nur mit der Spalte 2r durchgeführte 
A Rechnung der Reihe nach mit den Spalten 0, 2, 4, ;.., 2r durchführt. Bildet man nun mit 


Hilfe der letzten Gleichungen die Summe (18), so kommt 


1 
Ing Aal, Ze +1 420 @+1+2r+20(,)) 
ep A 1 EEE Tg gs 2,141 
Die in der geschweiften Klammer {},,r+ı enthaltene (r-+1)-gliedrige Summe läßt sich verein- 


fachen, wenn man die Binomialfaktoren () gemäß der Beziehung (6) — iR ') 4 ES H) 


aufspaltet. Dann erscheint {},,,+1 wegen (d) = Er —=1 und (W a 1) —= 1 zunächst 
in der Gestalt 
r—1 
Un zZ (—1 17-e(7 ') P),++1, 0 +21 Malz: ern) P,,r+1,0 > 
-oder, indem ınan in der a Summe o durch o Kr a und zusammenfaßt 


Damm >= a dee Vi ') [Pp,r+1,0+1 — Pr,r+1,ol- 


Nun ist ae re a Be: 
Ports, ee = 3 +3 +29). ver zer 
Damit erhält On ER ne die ‚Form! 


Be x Wir haben Bing eine Rekursonsforml zur ‚Berechnung von onnen, 
Bi. den Be SIE | re Ri 
4 IN N v Harn N Fe 
Yetert, Die ‚gesichte Summe ba infolgedessen 0 “2 
: au=lett+l)p+2r +1) ae 4038 N 


Für die Determinante Be A| Hnden wir demnach aus (19) als Anfang der! 
Er ige von n 


Ap,r41 = + n" Ad, r+1 I?» ir +1) (p+2r +1) +04 En = on Ex: [or 

Bilden wir damit die Ausdrücke 1}, und Mar+2 gemäß 13), so un mit ei) und un 

zweiter Näherung WE 
M In tr(2r-+ 3) 
Kb reeresnres 


2n+(r+1) (2r +3) 
TETZATZI 
Setzt man l;, aus (24,) in (14,) ein, so erhält man 


A me ya Es 
ZA‘, Jr( N) YN 1+ Ps ® RM u 1 
Da sich die Glieder mit n-! aus der Entwicklung herausheben, so bleibt mit (21,) in zweiter 
Näherung unverändert gegen (14,) als Gesetz der Fehlerfortpflanzung 

3 SATT RAR ee 
Hi » (R; N ) => a wer you YN a a a | 


Aus (15,) und (24,) erhält man in Verbindung mit (21,) für die Fälschung fy in zweiter Näherung 
„mit einem von r und n Be. Verbesserungsfaktor‘‘ gegen (15,) 


ag ge u RE a 
und u = 


a = Ber . 4). 


„rt 


Mar+2(n) 7 


; : Er Aus ( j ee 1 arts Br 
Rt; Die Gl. (25) und (26) gelten für R= 2r und R=2r +1, 0. Im Falle r=0 ist das nur von r 
b FR abhängige Produkt in (25) Tr ‚d.h. es hat den Wert 1 


6. Der Gültigkeitsbereich der zweiten Näherung 


Die im vorausgehenden Abschnitt entwickelten Fehlerformeln beruhen im großen und 
ganzen auf der Näherung (9,) für die Summen Sp. Die Vernachlässigung wirkt sich am un- 
günstigsten bei der in der Determinante Ds, ,+1 höchsten Summe 


Surta = 2ntrte|,— Tr ts re: } 


aus. Wählt man bei vorgegebener fester Ordnung R= 2r>2 der Ausgleichfunktion die Punkt- 
zahl n>4r +3, so liegt der relative Fehler von S4r+2 unter 5,6%. Da alle anderen Summen 
in den maßgebenden Determinanten besser angenähert werden, so darf man für 


nz4Ar +33 „Sn Sa 27) 
bereits ein brauchbares Ergebnis erwarten. Tatsächlich ist z. B. für r=1 und n=7 der Faktor 


für die Fehlerfortpflanzung genau Z,(7) = ii 75° 1,004, während die zweite Näherung (25) 


81 
den Wert L,(7) z > TI liefert. Die THE sind also völlig belanglos. Bei dem Faktor 


En für die Fälschung. en a ad as ne Bee ee wie aus Bild 1 an- E x 


ä ‚schaulich aa ‚Die ‚Kreispunkte ee den genauen Werten®) 
Ben? | ER) 
d j ; a Fi: muln) = 5-7 7: een X } 


wahre 2% obere ‚Kurve zu 


Pr 


A SELTTeTeTTTT 
Eee Na 
Bi „ gehört. An E in (27) fest- | Näterungerh) | | | | | 


gesetzten Gültigkeitsgrenze 
' n=7 ist der relative Fehler % 
Am,/m,=1/n? etwa 2%. Da 
es sich bei der Verwendung 
der Fehlerformeln jedoch 

stetsnurumeineAbschätzunge * 

2 des physikalischen und ma- 

thematischen Fehlers handeln 
kann, so sind auch größere 0 


Bey: 


BE 
UNSNNBERBNNEE 
EI IRSCUNNEREEE 
le 

>= 


! 


E Unterschiede zwischen stten- 2 2 

E ger und angenäherter Lösung Bild1. Der beim Ausgleichen mit einer Parabel r = d ee R= P und 3) ent- 

; zulässig. Es hat vom prak- stehende mathematische Fehler /y (2, N) = —(jal)® j I (nz)imı(n) ist bei fester 

ber. P . ER Bereichlänge nr = konst verhältnisgleich zu m, (n). Die Punkte geben m, (n) nach der 
a tischen Gesichtspunkt aus ae ausnünns! a EM ist Er Sa nis! Sep ee Kraus 
Bi a : 1 { Vz r. „große‘‘ Punktzahlen n. e untere Kurve entspricht einer mit Hilfe der Grenz- - 
# keinen ‚Wert, die Genauig lösung. Ne ) berechneten Näherung für das ‚‚Zwischengebiet‘‘ II. (Vgl. den Wortlaut 

keit bei der Berechnung der zu Bild 4 und 6.) 


Zahlenfakoren Za,(n) und i 
Mgr+2(n) zu weit zu treiben, weilsowohl die Beträge der ehe u Bis auch die Ablei- 
tungen «2"+2) der Erfahrungsfunktion oft nur der Größenordnung nach bekannt oder berechen- 


E bar sind. 

4 Beim geradlinigen Ausgleichen (r=.0) hat man einfach den Mittelwert: 

e. a 1 > | 

ä | > HEIL, TR pr SEE 2, DER N ec EEE (28,), 

f zu bilden und erreicht damit den günstigsten physikalischen Fehler für % 

Er ., N AR een: 28,). 
| ON 7 (28,) 


Das Fehlerverhältnis Sr (R; N)/fr(0;N) =1%, ist in Abhängigkeit von der Ordnung R der Aus- 
3 gleichfunktion in Bild 2 dargestellt. Aa ist dabei folgender Zusammenhang. Be- 
7 . zeichnet man die zum Erreichen eines vorgeschriebenen Fehlers fr(R; N)=konst notwendige 
‘# Punktzahl N mit N,, wenn linear ausgeglichen wird, und mit Nx bei Verwendung einer Potenz- 
a 
B 


3 reihe vom Grade R=2r (oder R=2r +1), so gilt 

5 m =— x ae ae Ta Te a NEE ern (29,) 
3 Es ist z.B. N,=K?®, falls fr(R; N) = konst = (1/K)u betragen soll. Bildet man nun den 
E Differenzenquotienten der Funktion Nr/N,, so erhält man mit Ar=2 

E. AlNRIN)._ Al F Lime 5) | 

1 TE eerg — 3; —.] = aan 11. 


Mit (21,) läßt sich das in der Kan 


A(NzIN,) _ = | ) (< \ n a ne U: = ae 
Ar -[5 5 AAN E IB are 2r=2 2r 4r 
- schreiben. Da der Ausdruck in [ ] die Wallissche Produktzerlegung für 2/x ist), so strebt der 
Differenzenquotient mit wachsendem r gegen den Grenzwert 


u N) 2 m ee (29,). 


r— © Ar 


8) Stange, a.a.0., 8.394, Gl. (28). 
ER. Courant, Ditferential- und Integralrechnung, Bd.1, 8. 181. Berlin 1927. 


® ns r ” # 
5 
el F 
A Bi 
in 
58 % 


Si 
a oe a 
| 6. r} 70 72 %n E 


Die Beziehung 


ANzIN.) = 0,637 
Ar = . 
ist aber schon von Anfang an recht gut erfüllt, wie die Fäller=1.und r=2zeigen. Esist nämlich 
dar Pat = A: 625 und 0, 633. 


Daraus folgt, daß die zum Erreichen eines vorgeschrie-: 
benen Fehlers fp{R) notwendige Punktzahl N% (nahe- 
zu) linear mit der Ordnung R der Ausgleichfunktion 
anwächst. Es ist mit sehr guter Näherung anfangs 


Na=M(t+ZR)... Eh: (29,), 


was in Bild 2 anschaulich zum Ausdruck kommt, und 
später wegen (29,) 


N2= N (R+1) er: 


2 7. Die Fehlergesetze für die Ableitung £° 
Die Ableitung an der Stelle 2= z; wird nach (1) 
1 
j a = ER RFEE (30,), 
wobei der Faktor 
R yo 
0 2 & 6 Rar 8 a4,= 9.3 Ay Lei ER (30,) 
Bild 2. Die zum Erreichen eines vorgeschriebenen yv=—n 


physikalischen Fehlers fp (R, X) = konst erforder- 3 
liche, Bunktzahl N =" + 1 wächst (nahesn) ein gewogenes Mittel aus den Ausgangswerten 2%; —_n, «+» 


reise benutzten Ersatztunktion. fin ist, Welches mit Hilfe der Einflußzahlen 
1 S a Sar 
; „| De ER Der 
»=2|: { ee a EA (30,) 
2,17 


Ar Fe 
gebildet werden muß. 
Die für die Fehlerfortpflanzung 


SBEN=EBe rs. N De (31) 


%n 
maßgebende Summe ZLY,= 3 A, der Gewichtsquadrate läßt sich nach einer ähnlichen Rechnung 


van 


wie vor Gl. (8) in der Form 


N Ds; r—i + (=) 
Li, = —— =1lim I], 3 wo oe 3 
2r D,, . Sn ( 2,) 
darstellen. Das für den mathematischen Fehler der Ableitung %;, 
’ 1 
s r+1) „27 Mt 
Ju (R;:N)= @r +ıI +1i al? T M3 +1 eh (33) 
maßgebende nichtverschwindende Moment M3++ı der Einflußzahlen ist 
’ D, r 
eg ER A cn Sp. 
M3r +41 ( 1) Des ee 22 ee (34,). 


Mit Hilfe der Formeln (10,), (23) und (20,) bzw. (20,) findet man aus (32,) und (34,) zunächst 


3 RB: ... eeelk 2n +(r—1) (2r-+3) 


Ka DR 
we 2 | Zrtllr—1)l 20 Hr Ar = Sn 


ce u ee "Te er 


KR, hr h er Br 

a er +) SE, N 
a ee — h 

f Mern=t ı "ar #3) Or 75): (Ar PTR ne Done 0 f 

Die Fehlerfdrmeln (31) und (33) für = Ableitung # x werden infolgedessen, ebenso ı wie &) 3 

und wöh in zweiter Näherung, 


5° 7: .9. (dr +1) EI 1 1 rd ; i 2 Te 
tz ne Baar VN a a 
EN a 55.79. (@r +1) Vu Bi er 
3 IR: i 5 9 4- .6- -(2r—2) m+2 (nz) KR 


und Nase 


2 pr ei i : 2r + AAz, : T A ER N 

ck R; N N irei) a nr ee Prary 27, 
Er gültig für R= BE ‚und u 27, r DR “2 
; Der günstigste mit geradlinigem Ausgleichen r—=1 bzw. R=1 und R= = erreichbare b 
Er go aie Fehler der Ableitung ist k 


1;N)= En RR en FE 37). 
fr P( )= YN+2 NT 5: ( ) 
5 [In dem Falle ist das in 1 (85) orans chende von r abhängige Produkt „leer“, d. h.es hat den 
>, Wert’ 1.] Gl, (85) läßt sich demnach auch in der Gestalt - 
3 | Je(R;N) 5-:7-9.,-@r+l). 
NR 3 Fri; N) 2-4-6.--2r—2)’ 
schreiben. Mit steigender Ordnung R der Ausgleichfunktion wächst der bhyskatsche Fehler fp 
< Ableitung monoton an. Er multipliziert sich beim Übergang von R=2r—1 bzw. 2r zu 
4 =2r-+1 bzw. 2r +2 mit dem von der Punktzahl N unabhängigen Faktor (2r +3)/2r. Der 
2 a überhaupt erreichbare Fehler fr(1; N) ist verhältnisgleich zum Fehler u der Aus- 
£ gangswerte und umgekehrt verhältnisgleich zur Wurzel aus der Gesamtpunktzahl und zur 
„Bereichlänge“ B=nr. Der mathematische Fehler fr ‚hängt bei fester Ordnung der Ausgleich- 
funktion im wesentlichen nur von der ‚„Bereichlänge“ B ab. 
Die in der früheren Arbeit hergeleiteten Formeln für Fehlerfortpflanzung und Fälschung 
erkennt man leicht als Sonderfälle der allgemeinen Gleichungen (25), (26) und (35), (36). 


E 8. Eine Eigenschaft der Binomialfaktoren 

E Für die weitere Untersuchung benötigt man eine Eigenschaft der Binomialfaktoren, die im 
Be --, folgenden bereitgestellt werden soll, damit später der Gedankengang nicht unterbrochen werden 
5 muß. Es gilt für ihre auf die Symmetrielinie bezogenen Momente!) der Ordnung p 


5 für alle ungeraden Werte von p; 
4 5(2n) = > Br —1jr (er RS ",)=0 für die geraden Werte p=0; 2; (39,) 
E ER 2n—2); . | 

3 und für p=2n 

5 ee De 2n i 
z an=d van (— 1) K "= a en a VE a Er (3%). 


2 


5 Die Momente s, der mit abwechselnden Vorzeichen versehenen Binomialfaktoren (n ) ver- 


4 schwinden bei gerader Ordnung p, solange O<p<2(n— 1) ist. Das erste nicht verschwindende 


Moment ist sn (2 n). 
Es ist selbstverständlich, daß mit s, auch jede lineare Funktion der Gestalt 


ar) tree) + +, An)=0 ....... 0. (40) 
ist, solange Pp<2n—2 bleibt. Insbesondere ist aus dem Bi 


S 2n 
Be RT a Ara Ay I As = 
=, v2»? — 1] [»? — 4] [? — (n — 2) 2] (—1) [ E " 0, 
da sich jede dieser Gleichungen auf die Form (40) bringen läßt. 
10) Willers Il, a.a.O., 8.168, Gl. (13). 


(al), 


Weiter br 


auchen wir im folgenden die Summe der Quadrate der Binomialfaktoren. Dafür 
gilt”) E% 


3 ed: REN} I EA “) 


Setzt man insbesondere m = 2 n, so wird nach einer einfachen Umformung 


42 (2n\? (an) 2n +1) (2n +3) (dn— 1) 
1) 1:3 :*-(2n—]1) Er ie.) 


yv=—n 


9. Der geglättete Wert ©, als Funktion der Differenzen 47? gerader Ordnung 


Die im Abschnitt 7 hergeleiteten Fehlergesetze gelten asymptotisch für „‚große‘‘ Punkt- 
zahlen n2,4r +3. Ist nun im Gegensatz dazu die Punktzahl n nur wenig größer als r, so besitzt 
die Ausgleichaufgabe, wie Willers für die drei Fälle n=r +1, n=r-+2 und n=r-+3 gezeigt 
hat, eine einfache Lösung. Wir wollen im folgenden insbesondere die für n=r-+1 geltende 
„Grenzlösung‘ benutzen, um ‚das unbekannte Gebiet‘ zwischen n=r+1 und „großen“ 
Werten von n zu überbrücken. Dazu haben wir zunächst den Zusammenhang mit der von 
Willers gegebenen Lösung herzustellen. 


Die Differenzen 42? der nicht geglätteten Ausgangswerte 2; sind?) 


+? 2 p 
Ar=23 (-IP# Ye en RL NER (44,). 
=—y7 
Bezeichnen wir zur Abkürzung die Summe von zwei symmetrisch zu z; gelegenen z-Werten mit 
&+„ also = my +8 fürv=1;2; °'n, und setzen & = 2, so werden die s 
Differenzen 
2 $ i 2p 
42? a = (—1)ptv + y Er a BT ti a (44,). 
v=0 


Andererseits ist der geglättete Wert #, wegen ,=4_,. 


Schreibt man (44,) der Reihe nach für p=0, 1, ... n auf, so erhält man (n +1)-Gleichungen, 
welche die (n-+1)-Differenzen gerader Ordnung 4f, ..., 42” der Zeile k mit den (n-+1) 


„Funktionswerten“ &b» +++» &%4n verknüpfen. Löst man dieses System nach den „Unbe- 

PRILEN. &+» auf, so findet man die &+, als lineare Funktionen der Differenzen At, 4 
oo. Ay 

FDA FDA F RD VE BEI 


Setzt man diese Ausdrücke schließlich in (45) ein, so wird auch der ausgeglichene Wert &; linear 
durch die Differenzen gerader Ordnung dargestellt, 


N BAR Ha En en (47). 


Wir wollen nun die Faktoren b,, . .., dan dieser Entwicklung wirklich berechnen. Die Auflösung 
des Systems (44,) nach den &.+, schreiben wir in der Form 


ala arte) 
ee en de EEE DO (48). 


p Y 


Für»=1ist die Lösung offenbar richtig. Daß'sie allgemein gilt, läßt sich durch den Schluß von n 
auf (n-+-1) zeigen, wobei man den Faktor a(??) in der Entwicklung (46) zweckmäßig in der 


Gestalt 
| +9—1 2 
nl ee ee [1] ER — AT Ta SR 


schreibt. Die zum Beweise von (48) notwendige längere Rechnung soll hier jedoch unter- 
drückt werden. 


1) Willers I, a.a.0., 8.171, G1. (80). 
2) Willers I, a.2.0., 8.72, Gl. (6). 


oder mit (49) ausführlich aufgeschrieben ea Er EEE ee 
EZ = 10 +22, u + 3 P3 ) a Z. „(2 —1)2, '+ 2 
v= RE RER Br 


AR Setzen v wir nun iR Ausdrücke (as) i in (45) ein, so in der auageeNeheuen ne 


& Beat, “ 2 rar) Mla+2S, "\48, Be: = 


Re lm, 2 & ET ml) | 


Der Faktor von 10— & == 2; hat 3 Wert > = a Bei "Verwendung einer EN ai. 
ar 


‚vom Grade R=2r (bzw. R=2r +1) kenn die Momente M, = N ” A, der Gewichte A, 


u.a.für alle geraden Werte von a zwischen 2 und 2r, einschließlich der Grönzen 13), Infolgedessen 
. haben die Faktoren von A}, ..., 4%" den Wert Null, da sie als lineare Funktionen.dieser "Momente 


darstellbar sind. Die.erste in ae Entwicklung (50) on ar Ditterane ist Besnaele 
Ar. Es wird also aus (47) 


Ex % = % ar Aue 17 2. am. Bi BE Be | om 


ee ae Berge 
2 Ten le re Bu Er (52). 


EREZE 


| Durch diese Beziehung ist die am Ausgangswert x; anzubringende Verbesserung A = — 


in einfacher Weise durch die Differenzen gerader Ordnung 43%, r +1<p<n, der unge- 
glätteten Funktionswerte &—_n, . .., %%-4n dargestellt. Wir sind damit auf andere Weise zu dem 
bereits von Willers gefundenen Ergebnis gelangt. ‘Die Beziehung (51) ist einesteils enger als 
bei W., da der Satz dort nicht nur für die ‚mittlere‘ Ordinate &; sondern (in etwas allgemeinerer 
Fassung) auch noch für die „Randordinaten‘ &+, bewiesen wird. Auf der anderen Seite ist 
die Beziehung (52), welche die Faktoren ds, mit den Gewichten A, verknüpft, hier allgemeiner, 
da sich W. bei der Berechnung der Faktoren u die eingangs genannten drei Fälle n=r-+J1, 


r-+-2 und r-+3 beschränkt. 


10. Die Grenzlösung für zwei überzählige Punkte; n=r +1 bzw. N=R-+3 
Wählt man beim ee die Punktzahl n=r-+-1, so vereinfacht sich die Entwicklung 
(51) zu | 
Te = %% + Arzı AR EN NE N RATE, ee (53), 
. sie enthält also nur noch eine einzige Differenz. In diesem besonderen Falle gilt: 


für die Gewichte “ für die Binomialfaktoren gemäß (39,) und (39,) 


r+1 s r+1 N 2r +2 = 1 r+1 
Z hn=- 1) 2:1) Ka = - gr =) et! 


v 
r+1 r+1 2) >) 
SET Se 


vl 


ar RR 
Ser, =0 > var (— 17 Es =!) 


vol 
und 
r+1 h 
2 3 arte), = Morte 


v-1 


r+ 
2 5 „er+2 (—_ 1) = rt 2) DR (54 


v1 


23) Stange, a.a.0., 8.387, Gl. (11,). 


..236 Stange, Das Bildungsgesetz für die Fehlerformeln use (8 Jalıjaug. 1948 
FE EEE EEE SEIFE GER ELBE ERREGER EB RB EG ar TE TE ET 


EB j 2 2 
Vergleicht man die r Gleichungen (54,) miteinander, welche in A, bzw. (— 1yr er ) linear 
Rp / . r+1-+v 
und homogen sind, so folgt aus der Übereinstimmung ihrer Koeffizienten 


EHI RN wal:2;...3 lH Bee (55)- 


Im Sonderfall der Grenzlösung sind die Gewichte A, (abgesehen vom Vorzeichen) den Binomial- 
Tfaktoren Sr % ,) verhältnisgleich. Mit (55,) folgt weiter aus dem ersten. und letzten Glei- 
chungspaar (54,) und (54,) 


2r +2 
a CE RSS Bu (55,) 
und 
NM, n=c(-1f!(2r +23)! =Ayıl2r +2)!» 2... 2... (555). 
Der mathematische Fehler fy wird infolgedessen gemäß (7,) bei zwei überzähligen Punkten 
Ju (R; R +3) = Ar+i aüart8) art = Ar Fe ET ee (56), 


wobei A25+? die Differenz der (im allgemeiner nicht bekannten) meßfehlerfreien z- Werte 
bedeutet. Vergleicht man (56) und (53), so sieht man, daß die Herabsetzung der Streuung der 
x-Werte mit Hilfe von (53) ganz auf Kosten der Fälschung fy geht. Die letzte Gleichung ist 
auch anschaulich unmittelbar einleuchtend. Wendet man nämlich das Ausgleichverfahren (53) 
auf meßfehlerfreie Werte &;, an, so besteht der Erfolg nur in einer Fälschung dieser Werte um 
den Betrag Ti, — 24 = Artı A47t?. Diese Fälschung kann nach dem ANSERF BEN" mit 
Hilfe von (56) wieder vollständig "rückgängig gemacht werden. 

Zur Berechnung des noch unbekannten Faktors c setzen wir (55,) und (55,) in (8,) ein. 
‘Dann kommt mit (43) 


2r +2 
er ee ER 1-3... @2r+h) (57) 
+1 ( At6R) (2r +3) (2r +5) "(dr +3) R 3 
»=-4y \r +14+9 
Damit wird nun endgültig für die a: der geglättete Wert 30) 
2r +1) 
Az zn moee BEN 


der zugeordnete physikalische Fehler 


ru) 
fp(R;R+3) = YA, u mit A-1-74=1-(/ ) 1:37 (OrHD) 055, 


4r +4 r+1/(2r +43) --- (dr+3 
| 2) (2r +3) (4r+3) 
und der mathematische Fehler 
YATbts a ei Be ; 
Ju(R; R-+3) = (—1) Gr +3)@rf5)- (dr +3) 4 TR RN (58;). 


Blld 8. Verlauf des physikalischen Fehlers fp(R,R + 8) im Grenzfall des Ausgleichens mit nur zwei „über- 
zähligen‘‘ Punkten in Abhängigkeit von der BE 3 = ?2r der abschnittwelse benutzten Ersatzfunktion 
renzlösung) 


14) Gleichung (58,) erscheint bei Willers (U, 8. 171 unten) in etwas anderer Gestalt. 


„D r Verlauf von Er R Fe in eek i von der Ordnung ER Auszleichtunktion 
aus Bild 3 ersichtlich und bestätigt noch einmal die bereits in der früheren Arbeit ausgesprochene . | 


‚ e Tatsache, daß das Ausgleichen mit Hilfe ganzer rationaler Funktionen höheren Grades nur dann 


v e‘ 
w L 


a nl 2 u lea nn na. > 


Da 


die Streuung der Ausgangswerte merklich herabsetzt, wenn man „große“ | Punktzahlen N wählt. r 
Bei nur wenig überzähligen Punkten sind die Kurven ‚zu DIA und ah sich u 


dessen dem Verlauf der fehlerhaften Punktreihe zu gut an. | , N 
Die bekannten RE 2 E n Br 5 Be 
4 IR tr Ag bzw. = ug. x r ER (8) 


ana 
N Bi 


Bee Nr sr 
PSGdEREEEEN 
> 

are 

RE 


Bild4. Aufteilung der (n,r)-Ebene in drei Teilbereiche I ee besitzt Eelre Lösung), II Biel 
gebiet) und III (Gültigkeitsbereich der asymptotischen. Formeln für große n). Das ‚‚Zwischengebiet‘‘ II kaun mit 
Hilfe der bekannten Grenzlösung (P,) fürn =r + 1 und der asymptotischen Näherung (?P,) für n=4r +3 „über- 
brückt‘“ werden. Beispiele einer solchen Rechnung sind für den physikalischen Fehler fp (R, N) aus Bild5 und für 


den mathematischen Fehler fy(R,N) aus Bild6 und 1 ersichtlich 


beim Glätten mit Hilfe einer Geraden (r=0; R=0) durch N=3, bzw. mit einer Parabel (r=1 
R=2) durch N=5 Punkte erscheinen als die einfachsten Sonderfälle von (58,). 

Bei nur zwei überzähligen Meßpunkten nimmt die Lösung der Ausgleichaufgabe die einfache 
Form (58,) bis (58,) an. Es ist in dem Falle also nicht erforderlich, die wesentlich mühsamere 
Berechnung der Summen $,(n) und der aus ihnen gebildeten Determinanten D,,, durchzuführen. 

Nach Bild 4 wird der erste Qua- 08 
drant der n,r-Ebene durch die zwei Ge- 
raden n=r+1 und n=4r-+3 in drei „EIRN) 
Bereiche I, II und III geteilt. Im Be- 0; 
reich I besitzt die Ausgleichaufgabe keine 07 
Lösung. Längs der gestrichelten Ge- 
raden n=r entartet sie zur ‚‚Interpola- 
tion“, da in dem Falle die Zahl 06 
N=2n-H1=2r-+1 der Meßpunkte ge- 
rade hinreicht, die (2r +1) Faktoren a, 
... Agr der „‚Ausgleichfunktion“ (1) so 
festzulegen, daß sie durch alle Punkte 
hindurchgeht. Längs der Geraden n= 
r +1 gilt die Grenzlösung (58,) und (58,), 
deren Verlauf für fpaus Bild3ersichtlich 
ist. Im Bereich III gilt mit ausreichender 
Genauigkeit die Näherung (25) und (26). 

Im Zwischengebiet Il ist man jedoch hin- 
sichtlich der strengen Lösungauf dienu-r 7377, 757775773 


05 


UNE TOR IE TE 
merische Berechnung von Fallzu Fall an- Bild 5. Verlauf des physikalischen Fehlers fp (R, N) fürr=2 


(bzw. R=4 und 5) im Zwischengebiet II über der Punktzahl n. 


gewiesen, die in der eingangs genannten 
Arbeit wenigstens für das praktisch wich- 
tige Teilgebiet durchgeführt wurde, wel- 
ches in Bild 4 geschrafft dargestellt ist. 


15) Willers I,a.a.0., 8.253, Gl. (14). 


Die Kurve entspricht einer aus Grenzlösung (P’) und asympto- 
tischer Lösung (P}) berechneten Näherung; die zum Vergleich. 


eingetragenen Punkte sind das Ergebnis der strengen Rechnung. 
(Vgl. den Wortlaut zu Bild 4.) 


nu Bl Fehlerabschät 
Bild 4 mit Hilfe der en r 


j Be, die re "Näherung, Kr 
N ar en ge +5) --(4r +3) 
ar it | Ex 
durch das Zusatzglied je zu erweitern und ns so zu bestimmen, ER die verbesserte zweite L BE p 


Merz2 (n)=1 yet 4 


Be = Näherung Ei 
Bi ar + ar +5) „+3 


r+1 
TE EEE ne ste 


Z 7 
an der Grenze n=r +1 mit der bekannten Lösung (98,) übereinstimmt. Dann muß ze 08), 
und (58,) gelten 
RER (2r +2)! _@r+2)! 
ee zen 
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Diese Gleichung liefert für die einfachsten Fälle r=1 und r=2 die Werte „= —2 a, hr © 
= —)9. Es ist Bu in besserer ee als we 


—2!(r +12. 


9 


Der daraus berechnete Verlauf des ESTER Fehlers fy wird Sr r=1lin Bild 1 und für 
1= a in Bild 6 jeweils Auch die untere Kurve NSS DE Diese Kurven passen sich dem wirk- 


ep" 4 6 8 7) 12 4 n 
0. 


Bild6. Zur Berechnung des mathematischen Fehlers Iy® N)w M., +2 (n) für 


r=2 (bzw. R=4 und 5) im Zwischengebiet II. Die obere Kurve entspricht der 
asymptotischen Lösung; die Punkte sind das Ergebnis der strengen Rechnung und 
die untere Kurve gehört zu der mit Hilfe der Grenzlösung (P’*) berechneten Näherung 


für das Zwischengebiet II. (Vgl. den Wortlaut zu Bild 4.) 


se 
er 
) „ 


lichen Verlaufe gut an, der durch die Kreispunkte veranschaulicht wird. Die Abweichungen 
zwischen der strengen Lösung und der mit Hilfe der Grenzlösung verbesserten zweiten Näherung 
überschreiten im ganzen Zwischengebiet I für r=1 (bzw. r=2) die Beträge 2,5% (bzw. 3,8 yA' 


nicht und fallen infolgedessen bei einer Abschätzung des mathematischen Fehlers gar nicht 
ins Gewicht. 


Eingegangen: 13. 11. 1948. 
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Beitrag 2 zur Enstchung von Wa rellen au Wind. 
Von Walter Wuest in Göttingen 


. Die Entstehung vo von Wasserwellen durch Wind wird auf eine Instabilität der ebenen Wasserfläche gegen- | 
über kleinen Störungen zurückgeführt. Dabei werden laminare Grenzschichten in Wasser und Luft berück- 
sichtigt und für die Grundströmung lineare und gekrümmte Geschwindigkeitsprofile zugrundegelegt. Die 
Rechnung ergibt einen wesentlichen Einfluß der Luftgrenzschicht und die a: der zuerst Ban 


‚Wellen hängt von der Grenzschichtdicke ab. 


thickness of ihe boundary layer. 


La formation d’ondes sur l’eau par l’action du vent est deduite de Pinstabilite de 1a surface d’eau ini- 


tialement plane. sous l’influence de pelites perturbations. La calculation comprend des couches limites 
laminaires ei quant au courant de ‚base, des profiles de vitesse lineaires et nonlindaires. Il resulte que la 


couche limite de l’air influence beaucoup n a et la longueur des ondes initiales depend de I’ epaisseur 


de la couche limite. 
IlokassiBaeTcH, YTO BOSHUKHOBEHHE BONAHHX. BONIH IPu BeTpe MO3KHO IIPHIIHCATB HEYCTOH- 


YUBOCTH TIAAKOK BOAAHON TIOBePXHOCTH IO OTHOIIEHHW K MAJIEHBKUM IHOMeXaM. YA4HTB-. 


BAIWTCA JIAMHHAPHLHE IOTPAHHYHHIe CIIOH BOAH M BOSAYX&; ' IIPHHHMAeTCH B pacyeT IPAMO- 
JIHHeÜHHÜÄ U KPHBOIHHeÄHHH IPohmN cKOPpocTef OCHOBHOTO TeyeHun. Us pacyera cnenyer, 
YTO BIHAHHe HOTPAHHYHOIO CAHOA BOBAYXA 3HAYMTEIBHO; INUHA BOBEURADUIGE "DONE: BABNOHT 
OT TONAIIHHEI HOTPAHHYHOTO CIIOA. 


sı. Beobachtungen und Windkanalversuche. 


Wenn die Luft gleichmäßig über eine ebene Wasserfläche strömt, bewirkt die Flüssigkeits- 
reibung zunächst nur eine Mitnahme der obersten Wasserschichten (,‚Windtrift“). Es bilden 


sich also sowohl in der Luft als auch im Wasser Grenzschichten aus. Schon bei sehr kleiner 
Windstärke wird die Spiegelungsfähigkeit der Wasseroberfläche stark vermindert und von einer 
gewissen Windgeschwindigkeit an beginnen sich kurze Wellen über die Oberfläche des Wassers 
zu bewegen, die sich immer mehr vergrößern, je länger die Windbahn ist. Dies läßt darauf 
schließen, daß die ebene Wasserfläche von einer gewissen Windgeschwindigkeit an instabil wird. 
Es ist naheliegend, hierbei an einen Zusammenhang mit der bekannten Instabilität lamjnarer 
Grenzschichten zu denken. Auf diesen Zusammenhang scheint auch die Beobachtung hinzu- 
weisen, daß ebenso wie bei laminaren Wandgrenzschichten die Temperaturschichtung der Luft 
von Einfluß auf das Stabilitätsverhalten ist. Denn nach einer Beobachtung von H.Frank[1] 
werden nachts infolge der stabileren Temperaturschichtung die Wasserwellen weniger leicht, 
d.h. erst bei höheren Windgeschwindigkeiten angeregt als bei Tage. 

Immerhin unterscheidet sich die Bildung von Oberflächenwellen doch in einer Hinsicht vom 


Instabilwerden einer laminaren Grenzschicht. Während nämlich die laminare Wandgrenzschicht 


bei noch so kleiner Strömungsgeschwindigkeit instabil wird, wenn nur die Grenzschichtdicke 
genügend angewachsen ist, sprechen andererseits die Beobachtungen und Versuche für 'eine 
untere Grenze der Windgeschwindigkeit, welche Oberflächenwellen hervorzurufen vermag. 

- In qualitativer Übereinstimmung hiermit konnte Lord Kelvin (William Thomson) [2] 
theoretisch nachweisen, daß oberhalb einer gewissen, im Vergleich zu den Beobachtungen aller- 
dings viel zu hohen Windstärke Oberflächenwellen auch ohne Berücksichtigung der Flüssigkeits- 
reibung alsö lediglich infolge dynamischer Druckänderungen angefacht werden können. Bei 
Annahme von Potentialströmung ergibt sich nämlich für die Wellengeschwindigkeit c bei einer 
Windgeschwindigkeit-U die Formel: 


= Be ee EN ee NIE 


Hierbei ist r das Verhältnis der von Luft und Wasser. Für Wasser von 15° und Luft von 20° 
und 760 mm Hg ist r = 0,001218. c, ist die einer bestimmten Wellenlänge zugeordnete Wellen- 
geschwindigkeit bei Windstille. c, hat ein Minimum von Comin = 23,2 cm/sec bei einer Wellenlänge 
A=1,72cm. Esist dies die Grenzwelle zwischen den kurzen Kapillarwellen und den langen Schwere- 
wellen. Nach Formel (1,1) wird die Wellengeschwindigkeit c dann konjugiert komplex, wenn 


1 
Beet 
r 


offenbar durch den Kleinstwert von c, gegeben, so daß man Uyin = 666 cm/sec erhält. 


.Cg. Die niedrigste Windgeschwindigkeit, bei der diese Instabilität auftritt, ist 


r 


The formation of water waves due to the action of wind is deduced from the insiabälty of the plane water „ 
. surface under the influence of small disturbances. The compulation considers laminar boundary-layers in 
water and air with linear and monlinear velocity profiles of the base current. It resulis that Ihe problem 
is especially influenced by the boundary kagpi3 in the air, and the. length of the initial waves depends on the 
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Die von verschiedenen Beobachtern [3], [4] in der Natur festgestellten Zahlenwerte für die 


kritische Windgeschwindigkeit betragen etwa 85 bis 110 cm/sec, liegen also erheblich unter dem 


Wert der reibungslosen Theorie. Allerdings sind diese Beobachtungen nur teilweise zu einem 
quantitativen Vergleich geeignet, weil durchweg die Änderung der Windgeschwindigkeit mit der 
Höhe über dem Wasserspiegel nicht beachtet wurde und daher die Angaben unzuverlässig sind. 
Zuverlässiger sind Messungen, die an Versuchsanordnungen mit künstlichem Luftstrom gewon- 
nen wurden. Solche Messungen sind von Th. Stanton [5], G. Weinblum [6] und W. 
Wuest [7] durchgeführt worden. 

Stanton benutzte für seine 
Versuche einen etwa 15 m langen 
geschlossenen Wasserkanal, wobei 
er durch zwei Fenster beobachtete, 
die 4,57 bzw. 12,48 m vom Anfang 
des Kanals entfernt waren. Die 
Kelvinsche Formel (1,1) konnte 
durch diese Messungen nicht be- 
stätigt werden. Vielmehr ist die 
einer bestimmten Wellenlänge zu- 
geordnete Wellengeschwindigkeit, 
wie aus Bild1 zu ersehen ist, merk- 
lich größer als bei Windstille im 
Gegensatz zur Kelvinschen Theorie. 


40 


Wellengeschwindigkeit [em/sek] 


Bild1. Abhängigkeit der Wellenlänge von der Wellengeschwindigkeit bei verschie- 
denen Windgeschwindigkeiten nach der Kelvinschen Formel (1.1) und nach Merkwürdig erweise S cheinen alle 


Messungen von Th. Stanton [5]. 


Meßpunkte unabhängig von der 
Windstärke auf einer einzigen Kurve zu liegen. Auch bezüglich der unteren Windschwelle er- 
geben sich wesentliche Abweichungen gegenüber der Kelvinschen Theorie. Bei einer Windge- 
schwindigkeit von 2 m/sec bilden sich nach Stanton ‚„Rippeln‘“, d.h. Wellen kleiner Wellen- 
länge, während sich größere Wellen erst von 2,5 m/sec an bilden sollen. In Bild 2 sind die an 
beiden Beobachtungsstellen gewonnenen Meßergebnisse aufgetragen. Die Kurven scheinen sich 
in der Tat bei 2—2,5m/sec Windgeschwindigkeit zu treffen. Aus Bild 2 erkennt man gleichzeitig, 
daß’die Wellenlänge mit zunehmendem Abstand vom Anfang des Beckens um so schneller an- 
wächst, je mehr die Windgeschwindig- 


® Stroboskop 
©. Plattenkamera 
+ Filmkamera: 


Windgeschwindigkeit [mjsek] 
Bild 2. Abhängigkeit der Wellenlänge von Ya 
verschiedener Mage Her „Windbahhn "nach, Windkanaivorsuchen. von ag Arch Wind an Ki e s nt br 
G. Weinblum [6] erhielt bei Versuchen, die in der Preußischen Versuchsanstalt für 
Wasserbau und Schiffahrt durchgeführt wurden, für die kritische Windgeschwindigkeit den Wert 
2 m/sec. Im Bereich von 1—2 m/sec konnte Weinblum ‚‚nur ein leichtes Vibrieren der 
Wasseroberfläche beobachten, das sich grundsätzlich von den später auftretenden Wellen unter- 
scheidet. Dabei war bezeichnend, daß bei den genannten geringen Windstärken selbst durch 
elastische Schwingungen angefachte kurze Wellen bald wieder abklangen.“‘ Zur Ergänzun 
dieser Versuche wurden vom Verfasser 1940 Messungen an einem Wasserbecken mit künstliche 
Luftstrom ausgeführt, wobei insbesondere die Vorgänge am Anfang des Beckens untersucht 


keit den kritischen Wert überschreitet. 


DE 2 0 a u a Bd 


Eu - Kar Für die Durchführung der Versuche 4 a eine Sofern) von 2m ne Be 
BR: 0,8 m Breite und 0,15 m Tiefe hergestellt. Um eine Reflexion der Wellen am Ende des Beckens 
„zu verhindern, war dieses alsflach ansteigendes „Ufer“ ausgebildet und aufgenietete Blechstreifen 


sollten die kinetische Energie der Wellen vernichten (Bild 3). Das Wasserbecken wurde im 
„Rauhigkeitskanal‘“ des KWI für Strömungsforschung i in Göttingen aufgestellt. Es war aller- 


dings mit dieser Anlage, die an sich für größere Windgeschwindigkeiten bestimmt ist, nur durch 


improvisierte Maßnahmen möglich, zu kleineren Windgeschwindigkeiten zu gelangen. Bei einer 
Windgeschwindigkeit von 2,5 m/sec trat eine deutliche Veränderung der Wasseroberfläche ein. 


Jedoch sind stationäre Versuche bei kleineren "Windgeschwindigkeiten nicht durchgeführt wor- 
den. Besonders auffällig war diese Änderung der Wasseroberfläche, wenn bei Ingangsetzung des 


Gebläses durch eine Schwingung der Windkanalanlage bei niedrigen Drehzahlen Kapillar- 
wellen von etwa 1,2 cm Länge hervorgerufen wurden, die aber bei Überschreiten der kritischen 
Windgeschwindigkeit ganz plötzlich zu Windwellen umschlugen. Während die Kapillarwellen 
sich durch eine ganz gerade Front auszeichneten, waren die Wellenkämme der ersten Wind- 
wellen ganz unregelmäßig. Die Wellenfrequenz wurde stroboskopisch gemessen. Zur, Bestim- 
mung der Wellenlänge waren auf dem Boden des Beckens verschiedenfarbige Streifen mit Ab- 
ständen von 1,7 bis 2,0 cm. angebracht, und es wurde bei stroboskopischer Beobachtung versucht, 
die Reflexe der Wellen mit einem der Streifensysteme zur Deckung zu bringen. Zur Kontrolle 
wurde außerdem noch die Wellengeschwindigkeit zwischen zwei Marken gestoppt: In einem 
Abstand von 30 cm vom Anfang des Beckens entfernt wurden folgende Werte gemessen: 


Windgeschwindigkeit cem/sece — 300 400 
Wellenlänge cm 1,70—1,75 1,75—1,80 
Wellengeschwindigkeit em/sec 26,5 "2748 


Eine wesentliche Rolle spielt auch die Beschaffenheit der Wasseroberfläche. Wenn der‘ 


Windkanal einige Zeit in Betrieb war, hatte sich auf der ‘Wasseroberfläche feinster Staub ab- 


‚gelagert. Eine derartige Staubschicht wirkt wie eine feste Wand, so daß sich im Wasser keine 


Grenzschicht ausbildet und die Oberflächengeschwindigkeit des Wassers Null ist, während sie 
sonst etwa 3% der Windgeschwindigkeit beträgt. Die Wellenbildung tritt erst bei höheren Wind- 
geschwindigkeiten ein. Nach einem anderen Versuch des Verfassers, der 1941 mit einem kleinen 
Leybold-Gebläse durchgeführt wurde, wirkt sich eine Petroleumhaut auf der Wasseroberfläche 
dahingehend aus, daß sich die ersten Wellen erst etwa bei 5 m/sec Windgeschwindigkeit bilden 
und die Wellenlänge der ersten Wellen wurde aus photographischen Aufnahmen zu 2,5 cm be- 
stimmt. 

Für die Art der Grenzschicht im Was- 
ser ist die Oberflächengeschwindigkeit des 
Wassers kennzeichnend [7]. Bei konstan- 
ter Windgeschwindigkeit ist die Oberflä- 
chengeschwindigkeit ebenfalls konstant, wenn 
sowohl Luftgrenzschicht als auch Wasser- 
grenzschicht entweder laminar oder beide 
turbulent sind und zugleich mit verschwin- 
dender Dicke beginnen. Ist dagegen die 
Luftgrenzschicht turbulent und die Wasser- 
grenzschicht laminar, -so wächst die Ober- > Entfernung vom Anfang des Beckens 
flächengeschwindigkeit immer mehr an, Di nn nee der Wdbahn für verschiedene Windgechwindigkelen 
schließlich die Wassergrenzschicht ebenfalls (nach Messungen von W. Wuest [6]) 
turbulent wird und die Oberflächengeschwin- 
digkeit einem konstanten Wert zustrebt. Die Oberflächengeschwindigkeit des Wassers kann 
aber leicht durch geeignete Schwimmkörper bestimmt werden. Man muß nur dafür sorgen, 
daß diese nicht aus dem. Wasser herausragen und dadurch, daß sie vom Winde angetrieben 
werden, einen falschen Wert anzeigen. Bei Windkanalmessungen des Verfassers, die mit 
der früher beschriebenen Versuchseinrichtung durchgeführt wurden, erwiesen sich Tropfen 
von Paraffinöl, das durch Anilinfarbstoff rot gefärbt war, als besonders geeignet. Bild 4 
zeigt einige Meßergebnisse. Die Luftgrenzschicht war bei den Versuchen turbulent und hatte 
schon eine gewisse Dicke (etwa ZEN erreicht, ehe sie mit der Wasserfläche in Berührung 

16 


Ober flächengeschw.| Windgeschw. 
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trat. Die Wassergrenzschicht ist zunächst laminar und wird erst nach Durchlaufen einer ge- 
wissen Länge turbulent. Nach den Versuchen tritt übrigens die Wellenbildung bereits ein, ehe 
die Wassergrenzschicht turbulent geworden ist. Nach dem Umschlag zur Turbulenz strebt das 
Verhältnis der Oberflächengeschwindigkeit des Wassers zur ungestörten 

dem Grenzwert 0,032 zu. 


$2. Die Stabilität der Wasserfläche bei Wind t 

Wir beschränken uns auf laminare Grenzschichten in Wasser und Luft und lagern 
dieser Grundströmung Störungsbewegungen, deren Amplitude so klein ist, daß eine Lineari- 
sierung der Bewegungsgleichungen hinsichtlich dieser Größen gerechtfertigt ist. Die ungestörte 
Grenzschicht möge sich über mehrere Wellenlängen der Störung nicht oder allenfalls so wenig 
ändern, daß in der nachfolgenden Theorie für die Grundströmung ein nur vom Wandabstand y 
abhängiger Mittelwert U(y) angenommen werden kann. Dabei wird aber nicht vorausgesetzt, 
daß die Grenzschichtdicke klein gegen die Wellenlänge der Störung ist, so daß also der Störungs- 
druck auch vom Wandabstand abhängig ist. j : 

Zur Erfüllung der Kontinuitätsgleichung wird die Störungsbewegung von einer Strom- 
funktion W abgeleitet, die als eine jn «-Richtung (Hauptströmungsrichtung) fortschreitende 
Wellenbewegung angenommen wird, deren Amplitude @ nur von y abhängt: 


Van Er oa (2.1). 
2 : 
Hierbei ist A = = die Wellenlänge und 8 demnach immer reell, während ce komplexe Werte 


annehmen kann, deren reeller Teil die Phasengeschwindigkeit der Störung und deren imaginärer 
Teil die Anfachung oder Dämpfung ergibt. Im folgenden soll der Durchgang zur Instabilität 
untersucht werden, so daß c im allgemeinen, d.h. sofern nicht ausdrücklich anders vermerkt, 
als reell vorausgesetzt wird. Für die Störungsamplitude erhält man aus den Navier-Stokesschen 
Differentialgleichungen durch Elimination des Druckes: 


WR U ET 2 Hg... 23. 4 


Hierbei bezeichnet v die kinematische Zähigkeit. Die Störungsgeschwindigkeiten ergeben sich 
aus der Stromfunktion zu 


DOT 
IR dy —— o (y) eßtz-eh) ee ee (2.3) 
9% N iß (ı- 
Vers Ran IR ei re Pe =. 124): 
Der reibungslose Störungsdruck ergibt sich durch Integration der Bewegungsgleichung zu: 
p=olÜ!’py—(T—-)p—gi,jefteed) or... 2.5). 
Dabei ist DER Rs 
u A a (2.6) 


die Auslenkung eines Flüssigkeitsteilchens infolge der Störungsbewegung, g die Erdbeschleuni- 
gung und o die Dichte. Mit der Normalkomponente der Geschwindigkeit und damit nach (2.4) 
mit der Stromfunktion besteht folgender Zusammenhang: Ü 
dt Lern. 
Sr UT LP le=00 ee (2.7), 
so daß 
1 


= — pl) ee (2.8). 


Bei Berücksichtigung der Flüssigkeitsreibung ist die Normal- und T j 
> angentialspa 
nach H. Lamb [8] gegeben durch: BETH 


/ 0v 
De BT tn a Hr (2.9) 
ir 9v Au = 2 
Platz) Ro teren... RIO 


Hierbei ist 4 = ev die dynamische Zähigkeit. 


2 
2): ' Ordnung und hat daher vier EREEN Bi 
| Reibungsglieder auf der rechten Seite ist die Differential- 
a on, nur a zweiter Ordnung. Zwei Lösungen können demnach als „reibungslose“ ei [5 = 
Be Lösungen angesehen werden. Im Gegensatz zu den restlichen beiden ‚„‚Reibungslösungen“ sind a 
+. ale nur langsam ‚mit. dem Wandabstand y veränderlich. Die vollständige Differentialgleichung Bi, bi 
r führt offenbar auf je vier Integrationskonstante für die Luft- und Wassergrenzschicht. Insgesamt 2 
.- ee er sind also acht Randbedingungen zu in um diese a zu Bu: bi “ 
Diese N renneth lauten: We ER N KEE 


Pi N a es ee u=v=0 a Zn BA ERDE RENT MEN ee R 

a RR ee ON D, Pry ‚stetig EB EEE DONE Sr. Et BR NEE ar a 

E: TR REN Pyv unstetigum den Betrag der Oberflächenspannung, RR Wi 
Rare zwischen Wasser und Iuft ....... RR .@ 13). 


Wenn das der Hauptströmung an der Grenztläche zwischen Wasser 
- und Luft einen Sprung aufweist, nu statt der a von v die Sn von & füry= 0 . 


zu fordern. n | 

$ 2 ae N 3. Anrondink auf einen Sonderfall | 7 en 6: 
E- R Die Rechnung soll zunächst an einem besonders einfachen Sonderfall euer werden, | & 2 

3 der sich von der Kelvinschen Theorie nur dadurch unterscheidet, daß für die Störungsbewe- 
gungen die Flüssigkeitsreibung berücksichtigt wird, während von einer Grenzschicht der Haupt- 
B:=; = apesthnng gänzlich abgesehen wird. Die Grundströmung ist demnach gegeben durch: : 
5 y>o0 ZUR = const i B e 
A | ee | A 
E | Die Störungsgleichung @. 2) hat daher bei IR der Bedingungen für = E50 B ” 
2 in Luft und Wasser folgende Lösungen: $ Be. 
j BE re E 
3 | EBE ORRL ER RER (3,1). Be 
; Die drei eahiingläsen (2.13) ergeben Folgehde Beziehungen: = 
4 HER. — — (U,—o)le(C +D) | je 
3 A+B+[YR FR, — 1] B=—(C +D)— [I —iRw — 1]D Bi 
E A+B-+iR;/2 B = uw/un (C +D) — uwlun iRwl2D.  . 
4 Hieraus errechnet man: | e 
| fer RellfI HR, —1) Ur 2 +2 + mOe—lNCHD) 25, E 


MER DS RER.) ee N) 
| = (A, +iA)(C +D) 

| Die Oberflächenspannung T an der Grenzfläche zwischen Wasser und Luft erzeugt im 
B Wasser einen Überdruck von: 


IT N 


Eh Be 27 iBßtz—ct) 
| Ta” +PTop(0)(U,— 0) e i 
Setzt man 
9 er 70: Tpß we 
=, 


Bew+ter, ewter 
wo c, die Kelvinsche Wellengeschwindigkeit bei Abwesenheit von Wind ist, so erhält man mil 


(2.5), (2.9) und (2.13): 


(fe)? = |1 HU; Zee ra) ln Ga a ee (3.3) 
en a 2 A| 
Helft +Wle 27 |A, a 
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oder Näherungsweise Rn = Ser <1, wenn c reell sein soll: 
w 
(ei)? = : g (Dile Bit 7 eultew + 01) * . ” “ > D ” ” “ * (3.4). 
A, r0 


A,ist aus (3.2) zu berechnen. Näherungsweise kann man wegen u, <um das erste Glied 
im Nenner vernachlässigen und Yl —iRy —1 = y—iRy setzen. Jedoch bereitet es keine 
grundsätzliche Schwierigkeit, mit den vollständigen Formeln (3.3) und 3.2) zu rechnen. Die 
Nähernngsformel ergibt: 


= Be >razai — I u a SE a En a ee . . ” 3.5 ” 
YR,=y2 7 +Tıle ı) es (3.5) 


Man beachte, daß man für U, = 0 nicht genau die Kelvinsche Formel bei Windstille erhält- 
Dies kommt durch die Forderung, daß die Tangentialgeschwindigkeit der Störung an der Grenz- 
fläche stetig sein soll im Gegensatz zur Kelvinschen Theorie, wo Potentialströmung voraus- 
gesetzt wird. In Bild 5 sind die berechneten Werte aufgetragen. Es ergibt sich eine Mindest- 
geschwindigkeit von 160 cm/sec bei einer Wellenlänge A= 3 cm. Hinsichtlich der Windschwelle 
ist die Übereinstimmung mit den Messungen als gut zu bezeichnen, wenn man bedenkt, daß die 
bei Windgeschwindigkeiten über 2 m/sec beobachtete Wellenbildung ja bereits eine gewisse An- 
fachung voraussetzt, 
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Bild 5. Kritische Windgeschwindigkeit in Abhängigkeit von der Bild 6. Aus Geradenstücken zusam- 
Wellenlänge bei Vernachlässigung der Grenzschichten in Wasser und mengesetztes Geschwindigkeitsprofil 
Luft (jedoch Reibungseinfluß auf die Störungen berücksichtigt) als Grundlage der Näherungsrech- 


nung (schematisch) 


$ 4. Lineare Grenzschichten in Wasser und Luft 

Wir versuchen im folgenden, die Grenzschichten in Luft und Wasser näherungsweise da- 
durch zu berücksichtigen, daß wir die Grenzschichtprofile durch Geradensiücke entsprechend 
Bild 6ersetzen. Dabeisollein Geschwindigkeitssprung in der Trennungsfläche zugelassen werden. 
en ie Wahl der Geschwindigkeitsprofile verschwindet die Krümmung durchweg, so daß 

also U”’{y) = 0. Durch Einführung der Wirbelgröße ® — p" — i Ö lei 
en 1) ß° 9 geht die Störungsgleichung 
U—c—ivß)o=—iv/ßo" (y) A . (4.1). 


Wir betrachten schwach gedämpfte Schwingungen, für die ce = G@— ivß. (Wegen der 
Kleinheit von v ß handelt es sich praktisch um indifferente Schwingungen, so daß wir in den 
späteren Rechnungen c wieder als reell voraussetzen können. Statt y wird eine neue Veränder- 
liche eingeführt, die so gewählt ist, daß n= 0 an derkritischen Stelle U = € ist. Wir setzen daher 

U-4=U'y)en, 
wobei die Konstante & so gewählt wird, daß (4.1) übergeht in 


1a In) ee 
Daraus folgt umgekehrt: i 


ET TE TREE TEN, OR SEN ER TER 
kai) ji EN th > Kada) ae 


DRAN 
oh 
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Er 


ae 1000 E otrag zn | ‚ur tt u yon Waren durch Wind 


"An der Grenzfläche y= 0 zwischen Weber ind Luft ist enlsprecheh Bild 6. U= 0: 


in der Luft bzw. U=T, im Wasser. In der: Ba alen Be also Pan Wert ae 0 


folgender Wert von ie 
Br, U,—eo U,oe 


No Fe — Then app. % 2 Beet oo (4.4). 
- Eine a: Beziehung ergibt, sich in der Wassere nzschicht für y=0. | 


Ein Fundamentalsystem der Lösungen der DNS UDE (4. 2) ist in bekannter | 
| Weise [9] durch Besselsche Funktionen gegeben: 


zn Iulz qea en) — (1ja)1" au). “Fi Be ehe 


bzw. mad fe Pre ae) = (1B)AR (UT 2p8) + ne 


Die Funktionen F\ und F/ sind von H. Schlichting [10] als: Reihenentwicklung 


dargestellt und für Werte. von Sa n< 4 tabelliert. Für die. Luftgrenzschicht kommt nur 
dasjenige Lösungsaggregat 


o=GF +HO,F/ . PER ER an) 


in Frage, das für große reelle positive 7 gegen Null geht, für die een dagegen 


dasjenige, das für große reelle negative n verschwindet. Diese Bedingung läßt sich erfüllen, 


wenn man die Lösung für große ») durch Hankelsche Funktionen darstellt. Durch Anwendung 


- der Übergangsformeln zwischen der Potenzreihenentwicklung durch Besselsche Funktionen und 


der asymptotischen Entwicklung durch Hankelsche Funktionen kann in der Lösung (4.7). die 


' eine Integrationskonstante eliminiert werden. Für die Luftgrenzschicht ergibt sich dann: se 
o,=(;L(F} —e'” 3r@ 7 (a/3)]7° (2/3) Dur REIFEN "J 


oder zahlenmäßig: 
o,=(,L ir ll, 1882 — ‚0,6860 ;) F7]- 


Die Indizes L und W beziehen sich jeweils auf die Luft- und Wassergrenzschicht, während 
die Größen ohne diese Indizes für beide gemeinsam gelten. 
Man hat also weiterhin die Differentialgleichung 


ey) -Boe)=CF'M) :..::.2 2... 2. (49) 
zu behandeln mit der allgemeinen Lösung 


(y)= l4+5 ze 5). nom or" man] + wi => Al ven or" man] er (4.10). 


Durch partielle Integration kann man hier noch eine Heihenentwiriäung nach der en Größe. 
vornehmen: 


o=4Aer + Bew + OF +höhere Gliedrine. ....... (41l) 
wobei a 
F= Sr pr (Mandy: cs une . . (4.12). 
Nm ee 


Bei der Berechnung der Eufforensin ist zu beachten, daß die kritische Stelle. n=0 
im Integrationsbereich liegt. Es ist daher erforderlich, mit den Potenzreihenentwicklungen von 
n = 0 aus zu rechnen, während bei der Wassergrenzschicht die asymptotischen Entwicklungen 
im allgemeinen ausreichend sind. Aus diesem Grunde werden folgende Umformungen vorge- 
nommen: 


n Fe 7 ’ ’ 
m = [ nmdm=—[Frdn ee Br el) 
Mo ee 
Aus der für Zylinderfunktionen J, (t) für beliebiges » gültigen Formel (vgl. z. B. [11]) 
2 
[A Wdt=1 
ö 


folgt aberr ‘=, 
Fa Pe = Fam = —(d —e%) 3UR T (AP) = 0,64395 — 111543 . .. (414). 
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Durch weitere Integration ergibt sich: 


M/F, n)dn= IF: (m) am —[ Fr (m) dm 


. . (4.15). 
n Kenn d 
-] [Km dndn- Kın-Fer= MM —Kın-Fes 
n=0 m=0 ; 
Hierbei ist aber nach Schlichting [10] S. 56, Formel (19): _ Ko 
00. Na h x . 2 
Pı> [ [ F)dmdn=—i ER a ar ae EEE 


y=0 = 
Insgesamt können wir also die Lösungen [..: die Luft- und Wassergrenzschicht in folgender 
Weise anschreiben: . 
a 
2 ae ee . (4.17). 
yw= AweV + Bwer! +CywFw 
Die ersten beiden Glieder. sind jeweils die reibungslosen Lösungen, das dritte Glied ist da- 
gegen die Reibungslösung. Die Reibungsglieder klingen im Vergleich zum reibungslosen Lösungs- 
anteil nach außen rasch ab, so daß man am Rande der Grenzschicht nur noch die reibungslosen 
Lösungsanteile zu berücksichtigen braucht. Außerhalb der Grenzschichten lauten aber die 
Lösungen: . 


Gum Ace DR ER are (4.18). 
m= Bev 
Am Rande der Grenzschicht, der durch y= 6, bzw.. y= — öw dargestellt sei, muß aber 
p und o’ stetig sein. Bezeichnet man in (4.17) die reibungslosen Lösungsanteile mit 
p* = Aexp (ßy) + Bexp (— By), 
so folgt aus der Stetigkeit von p und 9’ mit (4.18) nach Eliminieren von A und B: 


Pi (Or) IB P% (dr) = —1 | 
we Be ee (4.19). 
ow (— öw)lPpw(— im) = 1 A 

Durch die Beziehungen (4.19) kann eine der beiden Konstanten der reibungslosen Lösungen 
eliminiert werden, so daß insgesamt noch vier Integrationskonstanten zu ermitteln sind, wofür 
die Randbedingungen (2.12) und (2.13) ausreichend sind. 

Aus den drei Randbedingungen (2.12) findet man ganz entsprechend zu (3.2): 
(Fow/B em — OwF im) Or = 

_ Alu [Ow + 9 (U: o)I(V,—0)] Br—2 [1 —ur (U,—0)/uw (Us—o)] 
Dw —urluw Dr ke 
Hierbei sind folgende Abkürzungen gebraucht: 


9, = — Bor (O)/pz (0) 
Or=  Yw (0) BY (0) 
wobei 9* wiederum den reibungslosen Lösungsanteil darstellt. 
o, — FOulFon Ben — Ber Fon/For 
re (4:22), 
Fow/Fiw Bew — Bew FowlFow. 
Dy = ng ag pa a oh ae 2 © (4.23). 
w Bew FowlFow 
‚ „Der Index:0 deutet dabei an, daß die Werte der Funktionen F7, und Fy und ihrer Ab- 
leitungen an der Stelle y=0 (n= nor bzw. 7 = 7ow) zu nehmen sind. Nach (4.4) ist aber 


——_ 


Ber = Y— N V£ T, x 


Bon = nn VOR. 
3 


Eu 


EVER Kat 


BB Ba a Een a a ai ai 


Tabelle in Abhängigkeit VON Noz, zusammengestellt. Für die a: Ieltet, man aus 
„den en asy. ptotischen Reihen folgende Formeln ab: 
Folie Vom = VZELL +0,75 VE nm) + LBL EL 1m) ‚da. a & 


Y en FowlRomw = VELL 1.25 YEl— nom)? — 1,5988 5 (— Mom)? vr . 2). 
In Tabelle 1 sind diese A für einige Near von Now ausgewertet. % 


Bu; Ta b ellel 
Luftgrenzschicht 
NoL ForlFor V—NoL | | Y—nL Foul Feb 
0 has 0 ‚9220 — 0,484 Y nz (— 0,602 + 0,3 397 i) m 
0,5 — 1,178 — 0,5217 — 0,545 + 0,269i 
' 1,0 — 0,7490 —0, 18501 i 1 —0,892 + 0,350i 
1,5 — 0,5451 + 0,0163i — 1,254 + 0,320i , 
2,0 — 0,4152 + 0,1775 i — 1,698 +0,191i 
2,5 — 0,3294 + 0,3330i — 2,147 —0,16i 
3,0 | — 0,2789 + 0,4991i — 2,420 — 0,899 i 
3,5 — 0,2972 + 0,698 i - — 2,130 — 1,915i 
4,0 — 0,438 + 0,852 i — 0,986 — 2,428i 
oo — 0,7071 + 0,7071i — 0,7071.— 0,707 i \ 
Wassergrenzschicht 
NW 2 a a 720 700 
6 0,6603 — 0,7029 i 0,7123 + 0,7869 i 
8 0,6757 — 0,7053i 0,7093 + 0,7602 i 
10 0,6843 — 0,7062 i 0,7082 + 0,7456 i 
oo 0,7071 — 0,7071i 0,7071 + 0,7071i 


Für die weitere Rechnung legen wir foıgende Werte für das Geschwindigkeitsprofil ZU- 


grunde, die sich nach einer früheren Untersuchung des Verfassers [7] für Jaminare Grenzschichten 
. in Luft und Wasser ergeben: 


0,10; = U,U,— 0,0, 
öwlör = 1,69. 
Mit diesen Werten wird aber 
Dow =, (@) a — (13,12 2,86 - 1045)? — 34,0, 
NoL ÖL U3 
Man ersieht BT daß der Einfluß der Wassergrenzschicht nicht wesentlich sein kann, 


solange nicht etwa die Wellengeschwindigkeit nahezu gleich der Oberflächengeschwindigkeit 
des Wassers ist. Dieser Fall liegt aber weit außerhalb der auf Grund der Beobachtungen zu er- 


_ wartenden Werte für die kritische Windgeschwindigkeit. 
Wir haben schließlich noch die vierte Grenzbedingung (2.13) einzuführen, die wiederum. 


wie (3.3) in einen reellen und imaginären Bestandteil zerfällt, wobei das Verschwinden des reellen 
Teiles die Phasengeschwindigkeit, das Verschwinden des imaginären Teiles dagegen eine Bedingung 
für. den Durchgang zur Instabilität liefert. 
Setzen wir an Stelle von (4.20): 
(Fow/ß ew — Ow Fow) Cw = Aı Hi As: 


so ergeben sich somit folgende beiden: Gleichungen: 
2 U,+U,—2c 
r) # 1 — BR ar ru ey 


Ow U; 1 2: Re 


 0w +05 U;—c Bow ow +0, U,—ec U,—e Pö, 
[—0, (U, — 0) + (U; — co] A = 2 (uw — ur) B(Ow + A): - » » (4.27). 


3° SER Now. sind imme = negativ, Dt ‚die kritische Stelle innerhalb der Eufigrenzschiehte 2 
liegt. Für die Luftgrenzschicht sind die Funktionen Fon —roz und Y—nor Fox/For in 


Om Re ne . (4.36), 
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Näherungsweise kann die rechte Seite in (4.27) vernachlässigt werden, so daß man als 


Bedingungsgleichung erhält: 
A= 0. TR 


Wenn man als Veränderliche e/U,, ßör und e/ßv; wählt, liefert (4.28) eine Beziehung 
zwischen diesen drei Größen, so daß demnach nur zwei davon frei wählbar sind. Mit Hilfe der 
Gleichung (4.26) bestimmt man ferner o,/ßv, wo c, wie früher die Kelvinsche Wellengeschwin- 
digkeit bei Windstille ist. Andererseits ist aber wegen der bekannten Formel für die Kelvinsche 
Wellengeschwindigkeit c, die Größe c,/ßvz nur eine- Funktion.von ß allein. Auf diese Weise 
gelingt es, zu einem gegebenen Wertepaar c/U, und ßö, in’einem gewissen Wertebereich die Grö- 
Ben c, U,, ß und ö, für sich getrennt zu bestimmen. 


$5. Berücksichtigung der Profilkrümmung 

Wir haben bisher der Rechnung Geschwindigkeitsprofile zugrundegelegt, die sich aus Ge- 
radenstücken zusammensetzen, so daß die Profilkrümmung durchgängig verschwindet. W.Toll- 
mien [9] hat aber gezeigt, daß die Profilkrämmung für die Stabilität laminarer Wandgrenz- 
schichten von ausschlaggebender Bedeutung ist. Es(ist also auch beim vorliegenden Problem 
ein derartiger Einfluß zu vermuten. Die Profilkrimmung beeinflußt insbesondere die „‚reibungs- 
losen“ Partikularlösungen, die ja bisher nach (4.11) lediglich durch Exponentialfunktionen dar- 
gestellt wurden. Man erhält diese Partikularlösungen auch bei allgemeinen Geschwindigkeits- 
profilen verhältnismäßig leicht durch Reihenentwicklungen aus der reibungslosen Störungs- 
gleichung, die aus der vollständigen Störungsgleichung (2.2) durch Vernachlässigung der rechten 
Seite hervorgeht. Die kritische Stelle U = «,, wo Phasengeschwindigkeit der Störungsbewegung 
und Geschwindigkeit der Grundströmung übereinstimmen, ist eine singuläre Stelle (Stelle der 
Bestimmtheit) der reibungslosen Differentialgleichung. Zwei linear unabhängige Lösungen 
lassen sich als Potenzreihen P, und P,in der Umgebung der kritischen Stelle y= y; mit y, 
= (y— %,)/Ys als Veränderliche in folgender Form. darstellen: 


9% = YıPı(yı) 2 ; 
[73 [2 > 
U MO [ir Ur Tip: “ <o 


Dieses Hinzutreten eines imaginären Bestandteiles beim Durchgang. durch die kritische 
Stelle ist von W. Tollmien nachgewiesen worden. 9, besitzt für y, = 0 eine Singularität, 
indem dort @, logarithmisch unendlich wird. Damit dieses singuläre Verhalten verschwindet, 
muß die Reibung in einer kleinen Umgebung der kritischen Stelle berücksichtigt werden. Hier- 
durch geht auch der Phasensprung auf einer endlichen kleinen Strecke vor sich. 

Die Geschwindigkeitsprofile, die sich bei laminarer Strömung in Luft und Wasser aus- 
bilden, entsprechen den ‚‚Blasiusprofilen‘“ der Plattenströmung. Nach dem Vorgang von W. 
Tollmien nähern wir dieses Profil durch eine lineare und quadratische Funktion an (Bild 7): 


0 <ylö< 0,175 UT = 108 ydr RRED a (5.2) 
0,1755< yJö< 1,015 Ur IB Re an (5.3) 
1,015< yJö< BR N N Er: 


Hierbei ist die Bezugslänge ö die halbe Breite der dem Blasiusprofil einbeschriebenen Pa- 

rabel. Mit der Verdrängungsdicke ö* der Grenzschicht besteht folgender Zusammenhang: 
ö* = 0,341 6. 

Die Geschwindigkeit an der Trennungsfläche von Luft und Wasser wurde dabei gleich Null 
gesetzt, d.h. der Vorgang wird von einem mit der Wasseroberfläche mitbewegten Koordinaten- 
system 'aus betrachtet (vgl. Bild 7). Außerdem beschränken wir uns darauf, den Krümmungs- 
einfluß in der Luftgrenzschicht zw berücksichtigen, da ja nach dem Bisherigen der Einfluß der 
Wassergrenzschicht von untergeordneter Bedeutung ist. 

Im Bereich der linearen Geschwindigkeitsverteilung ergeben sich natürlich wie früher Ex- 
ponentialfunktionen für die reibungslosen Lösungen. Wegen des Zusammenhanges mit dem 
Folgenden schreiben wir diese Lösungen in der Form: | 

ac 
M= Ein P (y—yı) | 


%= ip (y—yı) 


FE Pe rn pe 


WAL a et a 


| Dabei ist die kritische Stelle y% en durch ee =6 c/1,68 v, Im wandnahen Bereich r 
stellen wir zwar U durch eine lineare "Funktion dar, doch dürfen wir nicht etwa U'= = const, 


U” —=0O setzen. Es ist vielmehr entsprechend den ersten N 
Gliedern der Blasiusschen Potenzreihe zu setzen: | 6 | ee 


U/U, = 1,68 yJö — 1,520 (y/6)®, 
so daß an der kritischen Stelle yz gilt: Er 
| SU;/U; = — 3,88 (eJU,)? . ©. 6) 


Im Bereich der parabolischen ne 
teilung ist die kritische Stelle nach (5.4) gegeben durch: 


= (1,05 —Y1—c/V)). 7): ulu, 
An Stelle des Wandabstandes y führen wir ‚olende | | 
neue Veränderliche ein: a) lineare Geschw. verteilung 
« | | bj parabolische m 
A SAN ER S | (5 8) K t a 
3 A en RL a NER. He en 2 


und setzen ferner: 


 A=pöyfl-elt,.. DR 


Bild7. - Grenzschichtprofil bei Berücksichtigung 


Auch wenn die kritische Stelle in den Inearen Ge- .der Profilkrümmung in der Luftgrenzschicht 


schwindigkeitsbereich fällt, ist y£ durch (5.7) gegeben und has socı 


stelit dann lediglich eine Rechengröße dar. Die reibungslosen Lösungen werden i im parabolischen 
Bereich durch folgende Reihenentwicklungen dargestellt 


$ı 
yl-e/U, 


- Das unterstrichene Glied verschwindet außerhalb der kritischen Stelle, also für 2 0. 
Die Koeffizienten sind gegeben durch: 


Dr ae l;,=-12; ,= Pl: „= — Bi/18; 

a, = — ß?/720 + 81/120; ag = — B?/2400 — 23 8} /10800; 

=1l;b4,=0; = —1-+Pi/2; db, = 1/8 + P/18; 

b* = 1/48 — 61 87/432 + 84/24; b5 = 1/192 + 29 82/5760 + 84/225; 

bg = 1/640 + 623 87/3384 000 — 2051 81/324 000 + 8 /720. 

Die Reibungslösungen klingen bereits im linearen Teil des Geschwindigkeitsprofiles ab. 

Sie können daher aus dem vorhergehenden Abschnitt übernommen werden, wobei lediglich für. 
U’ der neue Wert U’ = 1,68 U,/ö zu setzen ist. Am äußeren Rand der Grenzschicht, d.h. für 
Yyı= 1 geht die reibungslose Lösung über in 9 = const :exp (—fy), so daß dort 9 = — ßp. 
Durch diese Bedingung kann in der Lösung des parabolischen Geschwindigkeitsbereiches 


9= Al) 1 — e/U, + BP 


eine Integrationskonstante eliminiert werden, Auch an der Übergangsstelle vom parabolischen 
in den linearen Bereich muß @ und g’ stetig sein. Dabei ist das Hinzutreten eines imaginären 
Gliedes zu berücksichtigen, wenn die kritische Stelle durchschritten wird. Man erhält schließ- 
lich genau analog zu (4.21) ©, als Funktion von c/U, und ßö, jedoch im Gegensatz zu früher als 
komplexe Größe. Die weitere Rechnung unterscheidet sich nicht von der früheren, da für die 
Wassergrenzschicht die Krümmung nicht berücksichtigt werden soll. 


$6. Ergebnisse der Stabilitätsrechnung 
In einer früheren Arbeit [7] führte der Verfasser die Stabilitätsrechnung ohne Berücksich- 
tigung der Flüssigkeitsreibung durch, wobei ebenfalls lineare Geschwindigkeitsprofile in Wasser 
und Luft angenommen wurden. Wegen der Vernachlässigung der Zähigkeit treten als wesentliche 
Veränderliche nur c/U, und 86, in Erscheinung. Eine (4.27) entsprechende Beziehung liefert 


=a4Yı +, y} +3yJ +... HER ka Tefih are . (5.10) . 


%»=b tb Yı +beyi+... +6U;/Uryp: (log Nm... 6). 
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erhält, fünf RER die drei Instabilitätsbereiche abgrenzen, wobei der m 
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Wassergrenzschicht, der obere di 
Luftgrenzschicht bedingt ist. 1 
obere Instabilitätsgebiet führt auf W: B 
von c/U,, die der Größenordnung nach 
mit beobachteten Werten von c/U, ver- 
gleichbar sind. Man wird aus alkan 
Grunde der Luftgrenzschicht einen we- 
sentlichen Anteil bei der Entstehung von 
Wasserwellen zuschreiben müssen. Eine 
mit :Gl. (4.26) vergleichbare Beziehung EN 
liefert im reibungslosen Fall 0,/U, wnın 
mit c,wieder die Kelvinssche Wellen- 
geschwindigkeit bei Windstille bezeich- 
net wird. Jedoch ist die Wellenlänge: 
und damit ß noch frei wählbar. Die 
geringste Windgeschwindigkeit ist dann 
offenbar durch den Kleinstwert von c, 
gegeben, was einer Wellenlänge von Hi 


- A=1,72cm entspricht. Jedoch erhält 
'man keine brauchbaren Werte für die i 


untere Grenze der Windgeschwindigkeit, 
weil für genügend große Werte von ßö, 
die Windgeschwindigkeit gleich der Wel- 
lengeschwindigkeit wird. 

Im Gegensatz hierzu tritt bei Be- 


> Hücksichtigung der Reibung neben c/U, und B8, noch eine dritte Veränderliche auf, als welche 
Be." man etwa die Reynoldsche Zahl c/ßv, wählen kann. Hiervon sind zwei Veränderliche in 
E, einem gewissen Wertebereich frei 
R; | wählbar. Für die numerische Rech- 
K nung empfiehlt es sich, als Ausgangs- 


 Wertepaar die davon abgeleiteten 
Größen 1/@, und nz, zu wählen, 
die durch (4.21) bzw. (4.4) und 
(4.3) gegeben sind. In das Diagramm 
Abb. 9 mit den Koordinaten c/T, 
und‘ 8, sind Linien 1/0, = const 
eingezeichnet. Ferner wurden die Lö- 
sungskurven der Gl. (4.28) für meh- 
rere Werte von No, (No = — 0,5; 
—1,0; —2,0) eingetragen. Die Kurve 
nos = —2 bildet nahezu die Grenze 
des Lösungsbereiches, weshalb keine 
Werte darüber hinaus berechnet wur- 
den. Auch die übrigen ungefähren 
Grenzen des Lösungsbereiches wur- 
den eingetragen. Eine genauere Fest- 
legung der Lösungsgrenzen hätte die 


BIIA9, Stabilitätsdiagramm mit Reibung bei linearen Grenzschichten m Derechnung einer größeren Anzahl 


Wasser und Luft 


von Punkten erfordert, was aber 
wegen der Umständlichkeit der Rech- 


nung unterlassen wurde. Beim Vergleich mit dem reibungslosen Stabilitätsdiagramm der 
Abb.8 fällt auf, daß der Lösungsbereich nur einen Teil des instabilen Bereiches von Abb. 8 aus- 


!) Als Abzisse ist dort ß.öw gewählt, wobei aber im Gegensatz zur vorliegenden Arbeit öy = 1,1967, 
entsprechend turbulenten Grenzschichten im Wasser und Luft gewählt wurde. - 


FR = in a.) - 
DR a ” 


sungsbereiches in Bild 9 um- 


schließen also nicht. wie in 
Bild 8 den instabilen Bereich, 


sondern lediglich den Bereich, 


' in dem die Instabilitätskurven 
liegen. Im Gegensatz zum rei- 
bungslosen ‘ Stabilitätsproblem. 


erhält man bei weiterer Aus- 
rechnung die Größen U,, 6, 6, 


und ß getrennt. In Bild 10 und 
11 ist U, und d, über A=2n/ß 


aufgetragen. Als wesentliches 


Ergebnis erkennt man insbe- 


sondere aus Bild 11, daß die 
Länge der zuerst entstehenden 
Wellen wesentlich von der Grenz- 
schichtdicke beeinflußt wird, in 
dem Sinne, daß :um so längere 


‚Wellen entstehen, je größer die 
 Grenzschichtdicke ist. Bei ge- 


gebener Grenzschichtdicke wird 
nur ein bestimmter Wellen- 


längenbereich angefacht und 


zwar am leichtesten die kür- 
zesten Wellen dieses Bereiches, 
die in Bild 10 durch die aus- 
gezogene Linie dargestellt sind. 
Am leichtesten werden Wellen 
von etwa 3 cm Länge angefacht, 


was einer Verdrängungsdicke 


von etwa 1,5 mm entspricht. 
In Bild 12 und 13 sind die 
Ergebnisse der: Stabilitätsrech- 
nung bei Berücksichtigung der 
Profilkrümmung aufgetragen. 
Die Rechnung wurde durchge- 
führt für (c— U,)/(U,—UD;) 
= 0,234 und 0,2944. Bei Auf- 
tragung der Grenzschichtdick, 
über der Wellenlänge fallen 
beide Kurven nahezu zusam- 


men. Die Berücksichtigung der 


Profilkrümmung bedingt keine 
wesentlichen Änderungen des 
Lösungscharakters, weshalb von 
der. Berechnung weiterer Werte 
Abstand genommen wurde. 


$?7. Zusammenfassung 
Für das Stabilitätsverhal- 
ten.einer Wasserfläche bei An- 


Beitrag zur Entstehung von Wasserwellen durch Wind 


Air a! re ES Din 
Ale a ” E 
4 e ea 


. füllt. Lösungen existieren nur für kleine Werte von B ör- Während in Abb.8 zwei Lösunges 
'„ kurven den instabilen Bereich abgrenzen, ist in Abb. 9 der Lösungsbereich von.einer Kurven- 
schar. überdeckt, von der jedoch nur einzelne Punkte berechnet wurden. Die ‚Grenzen des Lö- ° 
Windgeschwindigkeit [m/sec] e Bun. ea, 


‚Bild 10. Kritische Windgeschwindigkeit in Abhängigkeit von der Wellenlänge bei 


linearen Grenzschichten‘ 


.d.Imm] ' 307 ey} 


Bild 11, Grenzschichtdicke 6 und ‘Verdrängungsdicke 8* in Abhängigkeit ‚von 
i der Wellenlänge ; 


 Windgeschwindigkeit [m/sec] 


. 1 RER | 7 ae ee [em] 


Bild12. Kritische Windgeschwindigkeit in Abhängigkeit von der Wellenlänge 
bei Berücksichtigung der Profilkrümmung : 


wesenheit von Wind sind die Grenzschichten (insbesondere diejenige in der Luft) von erheb- 
licher Bedeutung. Dabei bringt die Berücksichtigung der Flüssigkeitsreibung eine wesentliche 


A 


ie 


EN et 
a 


men, “7 > 5 
ya en cn Zee 
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Änderung der Rechenerg 

_ über der reibungslosen Theor To 
einer früheren Arbeit des Verfas: ers. 2 
handelt wurde. Die Wellenlänge der zu- 


ed 


von der Grenzschichtdicke ab und wächst 


der Profilkrümmung gegenüber der ver- 
einfachenden Annahme von linearen Ge- 
schwindigkeitsprofilen in Luft und Wasser 
keine grundlegende Änderung der Er- 
gebnisse. Die am leichtesten angeregten 
Wellen haben eine Länge von etwa 2 bis 
4cm. Windgeschwindigkeiten von etwa 
0,7 m/sec genügen bereits zuihrer Anfachung. Ein Vergleich mit beobachteten Werten ist deshalb 
nur unvollkommen, weil bei den Beobachtungen und Messungen wohl stets turbulente Grenz- 
schichten vorgelegen haben. Während nach Beobachtungen in der Natur, bei denen aber die Än- 
„derung der Windgeschwindigkeit mit der Höhe über dem Wasser nicht in Betracht gezogen wurde, 
kleinste Windgeschwindigkeiten von 0,85 bis 1,1 m/sec gefunden wurden, sollten nach Messungen 
am Windkanal erst Windgeschwindigkeiten von 2 bis 2,5 m/sec imstande sein, Wellen zu erregen. 
Wenn man annimmt, daß die Luftgeschwindigkeit ohne allmählichen Anstieg in einer Grenz- 
- sthicht über dem Wasser sofort den vollen Wert annimmt, wobei aber der Reibungseinfluß auf 
die Störungen im Gegensatz zur Kelvinschen Theorie der Potentialströmung berücksichtigt 
wird, ergibt sich merkwürdigerweise eine kritische Windgeschwindigkeit von 1,6 m/sec, also ein 
Wert, der einigermaßen vergleichbar mit den bei turbulenten Grenzschichten am Windkanal 
‚beobachteten Werten ist. Ein Zusammenhang könnte insofern vermutet werden, weil auch bei 
turbulenten Geschwindigkeitsprofilen der Anstieg der Geschwindigkeit auf einen endlichen Wert 
in einer außerordentlich dünnen Grenzschicht erfolgt. 
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‘ Bild 13. Grenzschichtdicke ö und Verdrängungsdicke ö* bei Berücksichti- 
: gung der Profilkrümmung 


KLEINE MITTEILUNGEN 


u 773 


erst gebildeten Wellen hängt wesentlich \ 


mit zunehmender Grenzschichtdicke n. 
Dagegen bringt eine Berücksichtigung 


Ein Kurvenblatt zur Berechnung von a cos? 
und — asin2a. 


Beim tachymetrischen Messen werden der zwischen 
den Meßfäden des Fernrohres erscheinende Abschnitt a 
der Meßlatte und der Höhenwinkel der Fernrohrachse a 
gemessen. Aus ihnen ergibt sieh die gesuchte Entfer- 
nung der Latte vom Fernrohr als L= acos?a, der 
Höhenunterschied vom Fernrohr zum Zielpunkt auf 


1 
der Latte als k = 7 asin2al). 


‘) Vgl. Jordan-Eggert, Handbuch der Vermessungskunde, 
2.Bd., 2. Halbband, 1933, 63, insbesondere 8,273, 


In der Regel reicht für die Berechnung von Z und A 
die Genauigkeit des Rechenschiebers aus. Verwendet 
man einen gewöhnlichen Rechenschieber, dann muß 
man dabei a und « je zweimal einstellen. Sonderschie- 
ber ersparen von diesen vier Einstellungen eine, näm- 
lich eine Einstellung von a. & muß man immer noch 
auf zwei recht verschiedenen Skalen einstellen. Weil 
zwischen diesen beiden Einstellungen Z abgelesen und 
notiert wird, muß man sich & merken, was ermüdet 
und eine Quelle von Fehlern ist. 

Bei Verwendung des auf dem nachstehenden Bild 
schematisch dargestellten Kurvenblattes stellt man 
a und & nur je einmalein und kann dann Z und A sofort 


» N 4% 
ath. Mech. Key 1% % 
8 Juli/Aug. 1949 


ablesen. Das Kurvenblatt enthält in ‚kartesischen 


Koordinaten (z,y) die Kurve L mit der Gleichung 
y=logceos’x und die Kurve mit der Gleichung 


y=log (= sin 2a) . Dabei wird a auf der x-Achse 
aufgetragen. Das Blatt enthält nur Punkte (a, y), für 


die —1<y=0 ist. Wegen der Periodizität der log- . 


arithmischen Skala ist die Kurve % dort, wo 


' Zeinda= 0,1 ist, d.i. bei a = 5°46’, abgesetzt. 


Boy. 


Kurvenblatt zur Berechnung von acos?& und a sin? 2% 


Aus diesem Kurvenblatt mißt man mit einem loga- 
rithmischen Maßstab, dessen Einheit die Einheit der 
y-Achse ist, die Größen Z und h auf folgende Weise 
heraus. Manlegtihn so auf das Kurvenblatt, daß seine 
Oberkante die Teilungen links und rechts an der Stelle 
& schneidet und die rechte Teilung auf ihm die Stelle « 
trifft. Dann liest man auf ihm beim Schnittpunkt mit 
der Kurve Z den Wert L = a cos? a, beim Schnittpunkt 


1 
mit der Kurve k den Werth = 70 sin2& ab. Das Bild 


zeigt den falla = 123, = 12°20’. Auf einer genauen 
Ausführung des Kurvenblattes mit einer Einheit von 
25 cm auf der y-Achse liest man mit dem zugehörigen 
Maßstab L = 117,2, h= 25,7 ab, was mit BRechen- 
schiebergenauigkeit stimmt. 

Es ist vorteilhaft, ein Lineal, am bequemsten eine 
Reißschiene zu verwenden, auf der man den logarith- 
mischen Maßstab in sich verschieben kann. Dazu muß 
das Blatt natürlich auf einem Reißbrett befestigt sein. 


Zur Erleichterung der Bestimmung des Stellen- 
wertes von h sind an der Kurve h die Stellen, wo 


1 
3 sin2& = 0,1, bzw. = 0,01 ist, bezeichnet. 


Für die meisten Zwecke wird eine Einheit von 12,5 cm 
auf der y-Achse und eine Darstellung des Winkelgrades 
durch 12 mm ausreichende Genauigkeit ermöglichen. 
Es genügt, das Blatt für Winkel bis zu 25° anzulegen, 
weil größere Winkel beim Tachymetrieren so selten 
vorkommen, daß es nichtins Gewicht fällt, wenn man 

«für sie mit dem gewöhnlichen Rechenschieber rechnet. 

Die Handhabung des Kurvenblattes wird auch von 
ungeschulten Kräften sofort begriffen und richtig aus- 
geführt. Auch für den mit dem Rechenschieber Ver- 
trauten ist sie übersichtlicher, schneller, bequemer und 
sicherer als der Sonderschieber. 


Innsbruck. L. Vietoris. 


Transformation der Grenzschichtgleichung 
bei dem Problem des schräg angeblasenen 


Zylinders durch Anwendung von Funktional- 


determinanten. Ye 
L.Prandtl!)hat bekanntlich im Falle des ebenen 


Grenzschichtproblems gezeigt, daß durch Einführung | 
5 ae 

des „Gesamtdruckes“;:@=p+ 0 7 die nicht-lineare 

Differentialgleichung des Problems auf eine Gleichung _ 


vom Typus der Wärmeleitungsgleichung für die Funk- 


tion G, nämlich auf: @x = const. + u -G@yy?) zurück- 


geführt werden kann. In dieser Gleichung sind alle 
auftretenden variablen Größen als Funktionen von r 


(i.e. die Wandbogenlänge des Profils vom Stand- 


punkt ab gerechnet) und % (i.e. die Stromfunktion 
des ebenen Problems) zu betrachten. o bedeutet die 
Dichte, p die Druckfunktion und u die «-Komponente 
des Geschwindigkeitsvektors m = {u, v}. 54 
Bedeuten nun x und y die üblichen Grenzschicht- 
koordinaten, x also die Wandbogenlänge und y den 


‘senkrechten Abstand von der Profilwand, so reduziert 


sich das Problem für den schräg angeblasenen Zylinder 
von unendlicher Länge auf folgende Differentialglei- 
chung für die dritte Komponente w (x, y) in Richtung 
der Zylindererzeugenden des Geschwindigkeitsvek- 
tors w; \ 

u(2,y) we (2, y)+v(a,y)wy(,y)=v'wyy(z,y) (1) 
(v = const. = kinematische Zähigkeit). Die' übrigen 
Grenzschichtgleichungen verhalten sich wie im zwei- 


‘dimensionalen Fall. Denn zu den Gleichungen des 


Problems der ebenen laminaren Grenzschicht um einen 
Zylinder: wur +v°vy= U(z)- U’(2)+v-uy, und 
ur + vy=0, wobei U(x) die Potentialströmungsge- 
schwindigkeit am äußeren Rande der Grenzschicht be- 
deutet, tritt noch eine dritte Geschwindigkeitskompo- 
nente win Richtung der 2-Achse (parallel der Zylinder- 
achse) hinzu. Die Komponente W der äußeren Poten- 
tialströmung, die entsprechend hinzukommt, ist 
konstant. Nach wie vor ist U= U(x). Daher ist 
w=i(z, y)‚,d.h.vonzunabhängig. Für (1) gilt somit: 

= {u(z, y),v(2,y), w(z,y}}. Man kann also?), 
nach Ermittlung von u und v aus den unverändert ge- 
bliebenen Gleichungen des ebenen Problems, w(x, y) 
aus der linearen Differentialgleichung (1) berechnen ®). 
Wenn dann nach wie vor y(%, y) die Stromfunktion 
des ebenen Problems bedeutet, so gilt: u= y,,v= 
—yz. Es wird also aus (]): 


YyUu—yUyeylyyıı 2... (2). 
Diese Gleichung schreiben wir mit Hilfe von Funktio- 
naldeterminanten und erhalten: 


D(w, y) wur D(x,v:wy) 

Dia,y) Dey) 
Auf Grund der bekannten Rechenregeln über Funktio- 
naldeterminanten ergibt sich dann weiter: 


Dwy) _ 
D (z, v wy) 
und schließlich, wenn wir x und yals neue unabhängige 
Variable einführen: 


Diwy) _ Di, vw) 
Dia,y) Dis,Yy) 


ı)L.Prandtl:,,Zur Berechnung der Grenzschichten‘“, 
Z. angew. .Math. Mech. Bd.18 (1938), 8. 77 —82. 

2) Die genaue Behandlung dieser Differentialgleichung und 
eine Anwendung findet man bei:H. J. Luckert, ‚„Überdie 
Integration der Differentialgleichungeiner Gleitschiehtin zäher 
Flüssigkeit“. Schriften des mathematischen Seminars und des 
Institutsfürangewandte Mathematikander Universität Berlin. 
Bd.1 (1933), 8. 245 —274. 

8) Dies r merkte bereits L. Prandtlin seiner Arbeit „Über die 
Reibungsschichten bei dreidimensionalen Strömungen“, enthalten in 
der nichtim Druck erschienenen Festschrift,, Albert Betz zum 60. Ge- 
burtstag, 25. Dezember 1945‘, S.134—141. 

) Bei den bekannten ‚‚ähnlichen‘‘ Lösungen des ebenen Problems 
(sieheetwabeiW.Mangler:,,Die ähnlichen‘ Lösungender Prandtl- 
schen Grenzschichtgleichungen'‘, diese Zeitschrift, Bd.23 (1943), 8.241 
bis 251) führt diese Berechnung auf einfache Quadraturen, Die nicht 
schwierige Rechnung möge hier übergangen werden. 


Ausgeschrieben bedeutet ice aber: ) Lin, 


so gibt es eine Funktion H(z,y) derart, daß 
und h al ri ” ” - ” ” ” Ss 
R ; \ Hy=w. FR a WE a er. sl Ta)r 
Wegen: y=y(z,y) ist y=p(a,y), wobei 9, #0 
vorausgesetzt wird. } 
‚Aus (6), (7) und (7’) folgt'nun sofort: 
v- (Hy)y = H; Beet ae (8). 


und wegen a =0 schließlich die Gleichung 
A 
YvYyy Hyy=H, RER  ! 4 .. (9). 


Hierbei ist 9y= Yyla, y) = yyl2,9 (z,y)) eine be- 
kannte Funktion in x und y, während H(z,y) aus (9) 


h zu berechnen ist. 


Damit ist auch unser Problem auf einen Ausdruck 


. vomTypus ah ag re zurückgeführt, 


wobei naturgemäß die Randbedingungen ganz andere 
sind als in dem — eingangs zitierten — Prandtl- 
schen Ansatz für das ebene Problem, 


Karlsruhe/Baden. Hans Schubart. 


Ein Verfahren zur näherungsweisen Längen- 
bestimmung unregelmäßiger Kurvenzüge. 

E. Barbier hat 1860 einen besonderen Weg zur 
Lösung des Buffonschen N Sen gefunden, mit 
Hilfe dessen es aber auch möglich ist, die Länge eines 
gegebenen, in beliebiger Weise verschlungenen Kurven- 
zuges zu bestimmen. Bar bier formuliert zunächst 
folgende Aufgabe: eine beliebige polygone, geschlossene 


' oder offene Linie werde willkürlich auf ein System 


äquidistanter paralleler Geraden geworfen. Welches 
ist die mathematische Erwartung E eines Spielers, der 
für jeden Schnittpunkt einen Franken erhalten soll ? 

Die mathematische Erwartung setzt sich additiv aus 
denjenigen für die einzelnen Teilstrecken zusammen 
und wird damit proportional der Gesamtlänge L des 
Linienzuges, d.h. E=kL. Im Grenzfall unendlich 
kleiner Teilstrecken gilt das gleiche auch für eine belie- 
big gekrümmte Kurve. Nimmt man nun an, diese sei 
ein Kreis vom Umfang ar, dessen Durchmesser gleich 
dem Parallelenabstand « ist, so werden bei jedem Wurf 
gerade 2 Schnittpunkte entstehen, womit die mathe- 


. matische Erwartung 


Do na Hash DER 
an 


wird. Für eine beliebige Linie erhältman damit = = 
E ist aber die Zahl der zu erwartenden 'Schnittpunkte, 
womit man die Kurvenlänge L = a bekommt, wenn 


man die Kurve mit einem System, paralleler Geraden 
überdeckt und die Schnittpunkte auszählt. Man hat 
hier also eine einfache Methode, die Länge besonders 


er A a U=VUyyı.. .. Kr 
Ist nun w (x, y) stetig und differenzierbar in x und y, 


Bee tems ermittelt. _ j 
as Ve | 
sein, wo die auszumessende Kurve nicht auf dem Papier 
gegeben ist, sondern etwa in einem o) Gerät 
erscheint. Ein ins Okular gelegtes Netz gestattet die 
Auszählung dann in ee Weise. 
Mittweida. Helmut Lindner. 
Eine einfache Näherungskonstruktion für die 
Zahl n. x 
DievonKochanski1685 beneN 
konstruktion für die Zahl m stellt mittels einer ein- 
fachen geometrischen Konstruktion den Zahlenwert 


ER = s108,... die. ler ih as 


Fehler von 6 + 105, der weit unterhalb der Schranken 


der Zeichengenauigkeit liegt. (Die Konstruktion ist 
z. B. dargestellt in Reutter, Darst. Geom. Bd. 1, 
8. 64.) 

Hier soll eine andere Näherungskonstruktion mit- 
Be werden, die zwar nicht so genau ist wie die von 

ochanski, dafür aber konstruktiv noch einfacher. 
Wir gehen dazu aus von dem Näherungswert für =/2, 
nämlich «= 1/8- (7 + Y31) = 1,5710... mit einem 
Fehler von 2-10*, also ebenfalls weit unter der 
Schranke der Zeichengenauigkeit. 


Diese Zahl läßt sich konstruktiv sehr einfach dar- 
stellen: Einen Kreis vom Radius 1 schneidet man mit 
einer Geraden, deren Neigung 45°, die den horizontalen 
Durchmesser links vom Mittelpunkt im Abstand 1/4 
trifft. Lotet man den Schnittpunkt der Geraden mit 
dem Kreis rechts vom Mittelpunkt (?) auf den hori- 
zontalen Durchmesser herab (Z), so stellt der linke 
Abschnitt des Durchmessers die Strecke ‚„‚a‘“ dar. 


Heppenheim a.d.B. Dipl.-Ing. Heinrich Müller. 


BUCHBESPRECHUNGEN 


Karl Reicheneder, Fehlertheorie und Aus- 
gleichung von Rautenketten in der 
Nadirtriangulation, (Veröffentlichungen des 
Geodätischen Instituts in Potsdam.) 98 Seiten, mit 
10 Abb., Berlin 1949, Akademie-Verlag. ‘Preis brosch. 
11,— DM. 

In dieser neuesten VeröffentlichungvonReichen- 
eder wird die Fehlertheorie und die strenge Ausglei- 
chung von Rautenketten in neuartiger, sowohl den 
Photogrammeter und Geodäten, als auch den Mathe- 
matiker interessierender Darstellungsweise behandelt. 

Mathematisch interessant dürfte u. a. die Darstellung 
der Querfehler in der ausgeglichenen Bautenkette durch 


algebraische Funktionen, und vor allem die weitgehende 
Benutzung des von geeigneten abkürzenden Bezeich- 
nungen Gebrauch machenden ‚‚Indizes-Kalküls“ sein, 

Der Photogrammeter wird es besonders zu schätzen 
wissen, daß durch Entwicklung eines auf die strenge 
Ausgleichung in Rautenketten gemessener Rich- 
tungen bezüglichen Verfahrens eine auf dem Gebiet 
der Radialtriangulation bisher vorhandene Lücke ge- 
schlossen wird, und daß durch schematische Abwick- 
lung des betreffenden Rechnungsgangs und angegebene 
Tabellen dabei gebrauchter numerischer Werte sich we- 
sentliche Erleichterungen für die praktische Anwen- 
dung des Verfahrens ergeben. Beim gegenwärtigen 


“Dem Sr ED Eee 


‚schule München), Die strenge Lösung 


Kr‘ Fü iR 


Stand der Radialtriangulation bleibt es immerhin frag- 


lich, ob diese Vorzüge genügen, um die praktische An- 


wendung des Verfahrens im Hinblick auf die recht mä- 


Bige Genauigkeit der betreffenden Beobachtungen ge- 


rechtfertigt erscheinen zulassen. Auchist zu bedenken, 
daß die in Betracht’ kommenden Beobachtungsfehler 
nicht ganz unbedeutende systematische Komponenten 
enthalten, die grundsätzlichen Voraussetzungen für 
strenge Ausgleichung nach der M.d.kl. Q. also nicht 
recht erfüllt sein mögen. 


Dresden. 
Dr.-Ing. Ludwig Föppl (o. Prof. a. d. Techn. Hoch- 
für 
die rollende Reibung. 44. mit 13 bildern. 


‚München. 1947. Leibniz-Verlag, bisher Verlag R. Ol- 


denbourg. Preis brosch. 6,— uM. 

Der Verf. greift das Problem der Verteilung von 

Druck- und Schubspannungen in der Berührungsfläche 
zwischen Rad und Schiene (die aus gleichem Material 
bestehend vorausgesetzt werden) von neuem auf und 
gibt eine rechnerisch befriedigende Darstellung. 
‘ Zunächst wird die bei Anwesenheit einer tangentialen 
Kraft sich ergebende Schubspannungsverteilung errech- 
net unter der Annahme, daß über die gesamte Berüh- 
rungsfläche hinweg Hatten eintritt, und gezeigt, daß 
unter dieser Annahme die Spannungen an der Ablaut- 
seite unbeschränkt werden; ein Haften über die ganze 
Fläche kann also nicht eintreten. .Sodann wird die 
Schubspannungsverteilung errechnet unter der An- 
nahme, daß über die ganze Berührungstläche hinweg 
Gleiten eintritt. 

Nach Erledigung dieser beiden Grenzfälle wird vor- 
ausgesetzt, daß die Berührungsfläche sich auf- 
teilt in einen vorderen Anteil, in dem Hatten, und in 
einen hinteren Anteil, in dem Gleiten (,‚Schlüpfen“) 
eintritt. Um eine Spannungsverteilung zu erhalten, 
die alle erforderlichen Randbedingungen erfüllt, werden 
(mit insgesamt sieben offenen Parametern) drei Span- 
nungsverteilungen von der Art, wie sie sich aus den 
vorher betrachteten Grenzfällen ergeben hatten, über- 
lagert. Aus den Betrachtungen über die Schubspannun- 
gen können sechs der eingeführten Parameter bestimmt 


werden; zur Bestimmung des letzten muß die Druck-: 


verteilung mit herangezogen werden. 

Aus der Druckverteilung läßt sich auch das Moment 
gewinnen. Man erhält für dieses Moment den folgenden 
Ausdruck, die ‚‚Hauptgleichung‘‘ der Tiheorie 
M=4AoN+ 

+47 %+2a—(l + 0) [Ya +(a/2n)] a), 

Ho 3/a+ (a/2n) —2 

darin bezeichnet N die Normalkraft, T.die Tangential- 
kraft, M das Moment, A die Länge der gesamten Be- 
rührungsfläche, & das Verhältnis von ‚‚Schlüpfantetl‘“ 
zu „„Häftanteil‘‘ der Berührungsfläche, n die Abkürzung 
Fe ee und 4. den "Reibungskoeffizienten 
beim Gleiten. Bei gegebenen Werten von N, T und M 
ist (1) eine Bestimmungsgleichung für a, so daß die 
Aufteilung der Berührungsfläche angegeben wird. 

Bei der Diskussion der',,Hauptgleichung‘“ zeigt es 
sich, daß & = 0,64 ein ausgezeichneter Wert ist, und 
daß für alle von Null verschiedenen Tangentialkräfte 
T dieser Wert sich einstellen muß. Das Verhältnis 
4#= T/N, bis zu welchem die Tangentialkraft ansteigen 
kann, wird zu u/u, = 1,136 ermittelt. 

Der Verf. zieht noch einen Vergleich des Schlupfes, 
der sich aus & = 0,64 errechnet, mit Messungen von 
G.Sachs(ds. Zs. 4 [1924] 8. 1) und findet gute Über- 
einstimmung. 

Schließlich wird auch noch der Fall untersucht, daß 
beim gebremsten Rad der Bremsklotz auf der Ablauf- 
seite so tief sitzt, daß die im Rad erzeugten Schub- 
spannungen auf der Ablaufseite (dem, Betrage nach) 
erheblich über jenen der auflaufenden Seite liegen. 

Von seinen Vorgängern hat der Verfasser vor allem 
die Tatsache übernommen, daß die Berührfläche sich 


 Buchb prechungen 


Bryan 


aufteilt in ein Haft- und ein Gleitgebiet; er übertrifft 
andererseits seine Vorgänger durch die Geschlossenheit 


einer neuen Darstellung und die Prägnanz der Ergeb-. 


nisse, BETEN 
Karlsruhe. Karl Klotter. 

* Dipl.-Ing. Dr. phil, Albert Betz(Dir. des Max-Planck- 

Instituts für Stromungsforschung und Prot. a. d. Univ. 

Göttingen, Konforme Abovildung. VIIlund 

359 8. mit 276 Abb. Berlin-Göttingen-Heidelberg 1948. 

Springer-Verlag. Preis brosch. 36,— DM. 


Dieses schon seit langem erwartete Buch ist aus Vor- 
lesungen entstanden, die 
Hörern hielt, die mit den Grundlagen der konformen 
Abbildung bekannt gemacht werden sollten, bei denen 


. aber keine weitergehenden Kenntnisse der Funktionen- 
theorie vorauszusetzen waren. Es geht von der geome- 
‚trischen Anschauung aus, ohne zunächst den Begriff der . 


komplexen Zahlen zu verwenden. Um die Zuordnung 


zweier Ebenen anschaulicher zu machen, benutzt esin 


der Hauptsache die flächenhafle Verteilung elektrischer 
Strome. Durch Überlagerung solcher Verteilungen und 
aurch Vergleich zweier Verteilungen gewiunt Ve 
verhältnismäßig Kompliziert& Abbildungen. Nachdem 
dabei gleichzeitig eine Reihe von Beginien der Funk- 
tionentheorie eingeführt sind, werden erst jetzı mit 
komplexen Funktionen eine Reihe von Abbildungen be- 
handeıt, die durch einfache Funktionen gewonnen wer- 
den. Dabei wird immer das Anschauliche und die An- 
wendung hervorgehoben. So werden die wichtigsten 
Aufgaben der Stromungstheorie untersucht, z. B. Strö- 
mungen um 'ragflügel, Abbildung kreisähnlicher Fi- 
guren, Zusammenhang zwischen Profiliorm und Ge- 
schwindigkeitsverteilung usw. U.a. wird hier ein bis- 
her nichv veröffentlichtes Verfahren von Wittich-Rie- 
gels zur Abbildung von Tragtlügelprofilen auf Kreise 
angegeben. Die beiden letzten Abschnitte sind der 
Abbildung durch doppelperiodische Funktionen und 
den freien Strahlen gewidinet. : 
So ist ein Buch entstanden, das einen Überblick über 
die wichtigsten Verfahren der konformen Abbildung 
und ihre hauptsächlichsten Anwendungen gibt, aus- 
gehend von leicht verstänalichen, anschaulichen Me- 
thoden, fortschreitend zu. recht weitgehenden Pro- 
blemen. In hervorragender Weise erfüllt das Buch sei- 
nen Zweck, allen, die konforme Abbildung gebrauchen 
eine Anleitung zu ihrer praktischen Verwendung zu 
geben. Darüber hinaus vermittelt es eine sehr an- 
schauliche Einführung in entsprechende Teilgebiete der 
Funktionentheorie. 
Dresden. 


Dr.R. Rothe, ehem. Prof. a.d. Techn. Hochschule 
Berlin, Höhere Mathematik für Mathe- 
matiker, Physiker, Ingenieure, 

Teil 1. Differentialrechnung. und 
Grundformeln der Integralrechnung 
nebst Anwendungen (Teubners Mathema- 
tische Leitfäden, Band 21). 8. Aufl. 2088. mit 161 
Abb. Leipzig 1948. Verlagsgesellschaft B. G. Teubner. 
Preis: kart. 6,20 DM. 

Teil UI. Integralrechnung, Unend- 
liche Reihen, Vektorrechnung nebst 
Anwendungen. (Teubners Mathematische Leit- 
täden, Band 22.) 6. Aufl. 208 S. mit 98 Abb. Leipzig 
1949. Verlagsgesellschaft B. G. Teubner. Preis: kart. 
6,20 DM. 

Teil IV, Heft 1/2. Unter Mitwirkung von Studienrat 
0.Degosang, Bielefeld: Übungsaufgaben 
und Lösungen, (Teubners Mathematische Leit- 
fäden, Band 33/34.) 5. Aufl. 109 8. mit-97 Abb. Leip- 
zig 1949. Verlagsgesellschaft B. G. Teubner. Preis: 
kart. 3,50 DM. 

Teil IV, Heft 3/4. Unter Mitwirkung von Studienrat 
O0. Degosang, Bielefeld: Übungsaufgaben 


der Verfasser 1932/33 vor 


erlasser . 


Willers. 
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mit Lösungen. (Teubners Mathematische Leit- 


fäden, Band 35/36.) 4. Aufl. 108 8. mit 51 Abb. Leip- 
zig 1949. u sgesellschaft B. G. Teubner. Preis: 
kart. 3,40 A 

Es wird allgemein sehr begrüßt werden, daß die aus- 
gezeichneten Lehrbücher von Rothe wieder vorliegen. 
Rothe hat es verstanden, den Stoff von Anfang an so 
auszuwählen und darzustelten, daß sich seit der ersten 
Auflage wesentliche Änderungen nicht erforderlich ge- 
macht haben. Diese wurden von Trefftz (Teill, 
1. Aufl. Z. angew. Math. Mech. VIII 336/337, 4. Aufl., 
Z. angew. Math. Mech. XI, 391, Teil Il, 1. Aufl., 
Z. angew. Math. Mech. IX, 165) und von Grüß (Teil 
IV, 3, Z. angew. Math. Mech. XVIL, 309) rear 
Man, kann nur wiederholen, was dort gesagt es 
Die Bändchen bieten einen im Verhältnis zum Umfang 
fast überreichen Inhalt an Theorie und vielseitigen An- 
wendungen in ausgezeichneter Darstellung. Allerdings 
machen sie es dem Benutzer nicht immer leicht, eignen 


sich daher weniger zum Selbststudium als zum Ge- zeigten Au: £ Du 

brauch neben Vorlesungen. Auf jeden Fall fordert ihr Mathematiker, ern auch bei open eine Inge- 

Studium intensivste Mitarbeit. Erfreulicherweise wird, Fang ehier we. geschätzt, Emp- 

wie der Verlag mitteilt, die Neuauflage in nächster Zeit !ehfung brigt. 22 

durch Wiedererscheinen des Bandes III, der fehlenden Dresden, > illers. 
Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen 


„Konrad K Prof. der Mathematik a 
ri Be) TE und Anw nd 


A 

zZ 

g a 

(1934) Seite 190 dieser 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


Ernst Schörner, Mathematische Plasta- 
skap-Modelle für den Unterricht in der räum- 
lichen und darstellenden Geometrie und in den Grund- 
lagen des technischen Zeichnens. (Herausgeg. von Max 
Wendt.) TeilIII: Flächen, die Geraden enthalten. 
14 Abb. Berlin und Rangsdorf. 1947. Verlagsgesell- 
‚schaft Dreyer & Ca. 


Dr.-Ing. Max Schuler (ehem. ao. Prof. a.d. Uniy. 
Göttingen, Mechanische Schwingungs- 
lehre. TeilI: Einfache Schwinger (Bücher der Ma- 
thematik und Naturwissenschaften, herausgeg. von Dr, 
Henry Poltz). 168 S. mit 123 Abb. Wolfenbüttel und 
Hannover. 1949. Wolfenbütteler Verlagsanstalt. Preis 
kart. 6,-- DM-West. 


Wolfgang Gröhbner (o. Prof, a.d. Univ. Innsbruck) 
und Nikolaus Hofreiter (ao. Prof. a. d. Univ. Wien), 
Integraltafel(ersterTeil):Unbestimm- 
te Integrale. VIII--1668. Wien und Inns- 
Er Springer-Verlag. Preis brosch. 18,— DM, 
6,40 $. 


Heinrich W, E, Jung, Mathematische Ab- 
handlungen. 88. (Hallische Monographien, her- 


ausgegeben von Otto Eisfeld.) Halle 1948. Verlag 
Max Niemeyer. Preis geh.: 2,20 DM. 

Dr.-Ing. W. Kaufmann (o. Prof. a. d. Techn. Hoch- 
schule München), Statik der Tragwerke 
(Handbibliothek für Baui jeure, von 

-- 314 S. mit 364 Abb. 


Robert Otzen). 3. Auflage. 
Berlin-Göttingen-Heide 1949. Springer-Verlag. 


Preis brosch. 25,50 DM. 


Berichtigungen: 


Bei der Besprechung der Tafelwerke von Emde in 
Band 29 Seite 128 dieser Zeitschrift ist versehentlich 
der Titel des ersten Bandes: 


Dr. Fritz Emde (Prof. für Elektrotechnik i.R.a.d. 
TH, Stuttgart), Tafelnelementarer Funk- 


tionen (2. Auflage), XIT+181S. mit83 Abb, Lei 
zig 1948. Verlag B. G. Teubner. Preis geb. 11,60 D 
fortgeblieben. 

H. Heinrich, Genauigkeit für die Halbschrittver- 
fahren der graphischen Integration. Z. angew. Math. 
Mech. 29 (1949) 51—52. 

In dem Aufsatz soll es überall Halbschritt-, nicht 
Halbschnittverfahren heißen. 


ZUSCHRIFTEN AN DEN HERAUSGEBER 


Berichtigung zurkleinen Mitteilung: C.Weber, 
Zur nichtlinearen Elastizitätstheorie. Z. angew. 
Math. Mech. 28 (1948), S. 189—199. 


Bei den Ausführungen hat sich eine Reihe von 
Schreibfehlern eingeschlichen, so daß die Ergebnisse 
ebenfalls zu berichtigen sind. Herr Dr. Haus Kau- 
derer hat mich auf diese Fehler aufmerksam ge- 
maclıt und ich spreche ihm hiermit meinen Dank aus. 
Die von Herrn Dr. Hans Kauderer berichtigten 
Gleichungen lauten wie folgt: 


Ä 0A 0A aA 
4.) = 87, + 2 (2, — Em) 3, 4 &p8g 37, 


2A, (ll 
Gl.(17) 1 — om = 2 — Em) Kun 2) 
8 Al). Eee 
Gl. (18) om = a 


al. (19) om = Im(&i +&,+ &) 


01—0 
61.(20) ZZ = Flle — em)? + (e,— em)? 
2. +(&— em)? - 

Die Annahmen Va und Vb lauten dann wie die ab- 
geänderten Gl. (19) und (20). 
. Der Hauptgedanke der kleinen Mitteilung. durch 
immer weitergehende einschränkende Annahmen über 
die Zerrungsenergie, das Gesetz für den Zusammenhang 
zwischen Zerrungstensor und Spannungstensor abzu- 
leiten, bleibt voll erhalten. . 


Schlewecke (Nette). C. Weber. 
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